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Pr´ologo 


La asignatura Fundamentos Matem´aticos de la Ingenier´ıa de la titulaci´on de 
Ingenier´ıa T´ecnica en Dise˜no Industrial consta de 4.5 cr´editos te´oricos, 1.5 cr´editos 
de pr´acticos y 1.5 cr´editos de laboratorio. Tiene car´acter troncal, anual y se imparte 
en primer curso. 
La asignatura est´a dividida en tres partes bien diferenciadas: ´Algebra Lineal, 
que se imparte durante el primer semestre y objeto del presente material, C´alculo 
Diferencial e Integral, que se imparte durante el segundo semestre, y las pr´acticas 
de laboratorio que se imparten en varios grupos, unos en el primer semestre y otros 
en el segundo semestre. 
Los estudiantes pueden acceder a la Titulaci´on de Ingenier´ıa T´ecnica en Dise˜no 
Industrial desde cualquiera de las siguientes opciones: 


COU. Opci´on A y C. 


Bachillerato LOGSE. Opci´on cient´ıﬁco-t´ecnica, artes y ciencias sociales. 


Acceso mayores de 25 a˜nos. Opci´on cient´ıﬁco-t´ecnica, artes y ciencias socia- 
les. 


FP II. Diferentes opciones. 


Ciclos formativos de grado superior. Diferentes opciones. 


M´odulos profesionales de nivel III. Diferentes opciones. 


aunque la v´ıa de acceso predominante es el bachillerato LOGSE. 


El material que aqu´ı presentamos cubre los aspectos te´oricos/pr´acticos de la pri- 
mera parte de la asignatura ( ´Algebra Lineal), sin pretender ser un manual exhaustivo 
de los contenidos de la misma, ya que los cr´editos te´oricos se imparten una vez por 
semana, durante el primer semestre, en sesiones de una hora y media, y los cr´edi- 
tos pr´acticos se imparten una vez cada dos semanas, durante el primer semestre, en 
sesiones de una hora. 


Este manual est´a dividido en cuatro temas y cubre los aspectos fundamentales 
del ´algebra lineal. En los dos primeros temas se hace un repaso de aspectos ya es- 
tudiados en el bachillerato cient´ıﬁco-t´ecnico, como son los sistemas de ecuaciones 
lineales, estudi´andolos desde el punto matricial, y la geometr´ıa del espacio tridi- 
mensional, ya que la v´ıa de acceso del bachillerato LOGSE no garantiza que los 
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estudiantes hayan estudiado matem´aticas en el bachillerato. Si bien, desde un pun- 
to de vista matem´atico, los dos primeros temas se deber´ıan estudiar despu´es del 
tercero (espacios vectoriales), hemos cre´ıdo conveniente seguir el mismo esquema 
del bachillerato debido al poco tiempo disponible para impartir la asignatura. En 
el tercer tema estudiamos los espacios vectoriales reales, centr´andonos en los espa- 
cios vectoriales R2, R3 y R4. El cuarto tema lo dedicamos a la diagonalizaci´on de 
matrices reales. 
En cada uno de los temas, adem´as de exponer el contenido te´orico, se muestran 
ejemplos sencillos que ayudan al estudiante a comprender los aspectos te´oricos, y 
ejercicios para profundizar en los mismos. 
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Notaci´on 


La notaci´on que emplearemos en este material es la est´andar en matem´aticas. 
Designaremos por N al conjunto de los n´umeros naturales 


N = {1, 2, 3, 4, . . .} 


por Z al conjunto de los n´umeros enteros 


Z = {. . . , −4, −3, −2, −1, 0, 1, 2, 3, 4, . . .} 


por Q al conjunto de los n´umeros racionales 


Q = {m/n : m ∈ Z, n ∈ N} 


por R al conjunto de los n´umeros reales y por C al conjunto de los n´umeros com- 
plejos. 
El resto de la notaci´on empleada se ir´a introduciendo en cada uno de los temas. 
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Tema 1 


Sistemas de ecuaciones lineales 


1.1. 
Introducci´on 


En este tema se aborda el estudio de los sistemas de ecuaciones lineales con 
coeﬁcientes reales. El ejemplo m´as sencillo es el de una ecuaci´on lineal con una 
inc´ognita (ecuaci´on de primer grado): 


a · x = b 
(1.1) 


donde a, b ∈ R, con a = 0, son conocidos. Evidentemente, la ´unica soluci´on de la 
ecuaci´on (1.1) es 


x = b 


a 
El estudio matricial de los sistemas de ecuaciones lineales permite escribir for- 
malmente estos sistemas como en la ecuaci´on (1.1). Adem´as, y m´as importante, 
permite estudiar su compatibilidad mediante el teorema de Rouch´e-Fr¨obenius y re- 
solverlos mediante la regla de Cramer. 


1.2. 
Matrices 


Como se ha indicado en la introducci´on nos vamos a limitar a las matrices con 
coeﬁcientes reales. En esta secci´on vamos a introducir las deﬁniciones m´as usuales 
sobre matrices y las operaciones que se pueden realizar con ellas. 


1.2.1. 
Deﬁniciones 


Dados m, n ∈ N, se llama matriz real de tama˜no m × n a una colecci´on de 
m · n n´umeros reales, que representamos en una tabla rectangular de m ﬁlas y n 
columnas 


A = 


 


 
 
 
 


a11 
a12 
· · · 
a1n 
a21 
a22 
· · · 
a2n 
... 
... 
... 
... 
am1 
am2 
· · · 
amn 


 


 
 
 
 
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donde aij ∈ R, i = 1, 2, . . . , m, j = 1, 2, . . . , n. Los n´umeros reales aij se llaman 
elementos de la matriz A. Las matrices las denotaremos con las letras may´uscu- 
las: A, B, C, . . . , y sus elementos con la misma letra en min´usculas y sub´ındices 
indicando la ﬁla y la columna que ocupan en la matriz. 
Se llama l´ınea de una matriz a una ﬁla o a una columna de la misma. 
El conjunto de todas las matrices reales de tama˜no m × n lo representamos por 
Mm×n(R). 


Ejemplo 1.2.1 
La siguiente matriz A ∈ M4×3(R) 


A = 


 


 
 
 


−1 
0 
1 
2 
−3 
4 
7 
5 
9 
4 
2 
−3 


 


 
 
 


El elemento que ocupa la tercera ﬁla y la segunda columna de la matriz A es: 
a32 = 5. 


La matriz nula es aquella matriz que tiene todos sus elementos iguales a cero, 
y se representa por O. 


Ejemplo 1.2.2 
La siguiente es la matriz nula de tama˜no 4 × 3: 


O = 


 


 
 
 


0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


 


 
 
 



Si m = 1, decimos que la matriz A es una matriz ﬁla, en cuyo caso adquiere la 
forma 
A = 
a11 
a12 
· · · 
a1n 


Ejemplo 1.2.3 
La siguiente matriz A es una matriz ﬁla: 


A = (−1, 2, 0, 0) 


En la matriz A se han separado sus elementos con comas para que no haya confu- 
si´on. 


Si n = 1, decimos que la matriz A es una matriz columna, en cuyo caso ad- 
quiere la forma 


A = 


 


 
 
 
 


a11 
a21 
... 
am1 


 


 
 
 
 
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Ejemplo 1.2.4 
La siguiente matriz A es una matriz columna: 


A = 


 


 
 
 


2 
1 
3 
4 


 


 
 
 



Si m = n, decimos que la matriz A es una matriz cuadrada, en cuyo caso adquiere 
la forma 


A = 


 


 
 
 
 


a11 
a12 
· · · 
a1n 
a21 
a22 
· · · 
a2n 
... 
... 
... 
... 
an1 
an2 
· · · 
ann 


 


 
 
 
 


En este caso decimos que A es una matriz de tama˜no n. 
El conjunto de las matrices reales de tama˜no n lo representamos por Mn(R). 


Ejemplo 1.2.5 
La siguiente matriz A ∈ M3(R) 


A = 


 


 
2 
−1 
3 
0 
4 
2 
7 
−4 
9 


 


 



En una matriz cuadrada 


A = 


 


 
 
 
 


a11 
a12 
· · · 
a1n 
a21 
a22 
· · · 
a2n 
... 
... 
... 
... 
an1 
an2 
· · · 
ann 


 


 
 
 
 


decimos que los elementos a11, a22, . . . , ann forman la diagonal principal de A. 


Ejemplo 1.2.6 
Dada la matriz 


A = 


 


 
5 
7 
9 
3 
4 
0 
−6 
2 
−1 


 


 


los elementos 5, 4 y −1 forman la diagonal principal de A. 


Una matriz diagonal es toda matriz cuadrada cuyos elementos no situados en la 
diagonal principal son todos iguales a cero. Si D ∈ Mn(R) es una matriz diagonal, 
entonces, en ocasiones, la representaremos as´ı: 


D = diag [d11, d22, . . . , dnn] 
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Ejemplo 1.2.7 
La matriz 


D = 


 


 
1 
0 
0 
0 
2 
0 
0 
0 
0 


 


 = diag[1, 2, 0] 


es una matriz diagonal de tama˜no 3. 


La matriz identidad I es una matriz diagonal cuya diagonal principal est´a for- 
mada por unos 


I = 


 


 
 
 
 


1 
0 
· · · 
0 
0 
1 
· · · 
0 
... 
... 
... 
... 
0 
0 
· · · 
1 


 


 
 
 
 


Ejemplo 1.2.8 
La matriz 


I = 


 


 
1 
0 
0 
0 
1 
0 
0 
0 
1 


 


 = diag[1, 1, 1] 


es la matriz identidad de tama˜no 3. 


Ejercicio 1.2.1 
Obt´en las siguientes matrices de M4x3(R): 


(a) aij = 


1 
si i = j 
0 
si i = j 


(b) aij = 


i + j 
si i = j 
i − j 
si i = j 



Dada una matriz de tama˜no m × n 


A = 


 


 
 
 
 


a11 
a12 
· · · 
a1n 
a21 
a22 
· · · 
a2n 
... 
... 
... 
... 
am1 
am2 
· · · 
amn 


 


 
 
 
 


consideramos ciertos ´ındices de ﬁlas: i1, i2, . . . , ip, siendo 


1 i1 < i2 < · · · < ip m , 
1 p m 
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y ciertos ´ındices de columnas: j1, j2, . . . , jq, con 


1 j1 < j2 < · · · < jq n , 
1 q n 


La matriz de tama˜no p × q obtenida a partir de A suprimiendo las ﬁlas y las 
columnas distintas a las consideradas, se llama submatriz de tama˜no p × q de la 
matriz A. 


Ejemplo 1.2.9 
Sea la matriz 


A = 


 


 
 
 
 
 


3 
−1 
2 
5 
4 
7 
−5 
0 
6 
9 
1 
−2 
−4 
8 
4 
−7 
5 
4 
2 
−1 


 


 
 
 
 
 
∈ M5×4(R) 


Considerando la segunda y cuarta ﬁlas, y la primera, segunda y cuarta columnas, 
obtenemos la matriz 


B = 
4 
7 
0 
−4 
8 
−7 


∈ M2×3(R) 


que es una submatriz de A. 


Dada A ∈ Mm×n(R), se llama matriz traspuesta de la matriz A a una matriz 
B ∈ Mn×m(R) cuyos elementos bji son: 


bji = aij , 
i = 1, 2, . . . , m 
j = 1, 2, . . . , n 


La matriz traspuesta de A la representamos por At. 


Ejemplo 1.2.10 
La matriz traspuesta de la matriz 


A = 


 


 
2 
−1 
3 
9 
4 
5 
7 
8 
−1 
0 
2 
6 


 


 ∈ M3×4(R) 


es 


At = 


 


 
 
 


2 
4 
−1 
−1 
5 
0 
3 
7 
2 
9 
8 
6 


 


 
 
 ∈ M4×3(R) 



Diremos que una matriz cuadrada A es una matriz sim´etrica, si At = A. 
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Ejemplo 1.2.11 
La matriz 


A = 


 


 
1 
2 
3 
2 
4 
5 
3 
5 
6 


 


 


es una matriz sim´etrica. 


1.2.2. 
Operaciones con matrices 


Suma de matrices 
Dadas A, B ∈ Mm×n(R), se llama suma de A y B a una matriz C ∈ Mm×n(R), 
tal que 


cij = aij + bij , 
i = 1, 2, . . . , m 
j = 1, 2, . . . , n 


La matriz suma de A y B la representamos por A + B. 


Ejemplo 1.2.12 
Dadas las matrices 


A = 
−2 
1 
3 
0 
5 
6 


, 
B = 
3 
−2 
4 
−5 
−1 
3 



se tiene que 


A + B = 
1 
−1 
7 
−5 
4 
9 




Producto de un n´umero real y una matriz 
Dado α ∈ R y dada A ∈ Mm×n(R), se llama producto de α y A a una matriz 
B ∈ Mm×n(R), tal que 


bij = α aij , 
i = 1, 2, . . . , m 
j = 1, 2, . . . , n 


La matriz producto de α y A la representamos por α · A. 


Ejemplo 1.2.13 
Dados α = 2 y 


A = 
−1 
2 
0 
3 
5 
−6 



se tiene que 


α · A = 
−2 
4 
0 
6 
10 
−12 
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Dadas dos matrices A, B ∈ Mm×n(R) se deﬁne la matriz A − B como: 


A − B = A + (−1) · B 


Ejercicio 1.2.2 
Considera las matrices: 


A = 
1 
−2 
3 
0 
4 
−6 


, 
B = 
−2 
0 
1 
−5 
3 
7 



(a) Calcula 2A − 3B. 


(b) Encuentra una matriz C tal que A + 2B − C = 0. 



Ejercicio 1.2.3 
Calcula x, y, z y w si se veriﬁca que: 


3 
x 
z 
y 
w 


= 
x 
−1 
6 
2w 


+ 
4 
x + y 
z + w 
3 




Producto de matrices 
Dadas A ∈ Mm×p(R) y B ∈ Mp×n(R), se llama matriz producto de A y B a una 
matriz C ∈ Mm×n(R), tal que 


cij = 


p 


k=1 
aik bkj , 
i = 1, 2, . . . , m 


j = 1, 2, . . . , n 


La matriz producto de A y B la representamos por A · B. 


Ejemplo 1.2.14 
Dadas las matrices 


A = 


 


 
 
 


2 
1 
2 
0 
−3 
1 
4 
−2 
3 
7 
0 
2 


 


 
 
 ∈ M4×3(R) , 
B = 


 


 
1 
−1 
0 
2 
2 
−3 


 


 ∈ M3×2(R) 


se tiene que 


A · B = 


 


 
 
 


2 
1 
2 
0 
−3 
1 


4 
−2 
3 


7 
0 
2 


 


 
 
 · 


 


 
1 
−1 


0 
2 


2 
−3 


 


 = 


 


 
 
 


6 
−6 
2 
−9 
10 
−17 


11 
−13 


 


 
 
 ∈ M4×2(R) 


Llamando C = A · B, el elemento c32 = −17 se obtiene multiplicando elemento 
a elemento, los elementos de la tercera ﬁla de la matriz A por los elementos de la 
segunda columna de la matriz B y luego, sumando: 


c32 = 4 · (−1) + (−2) · 2 + 3 · (−3) = −4 − 4 − 9 = −17 
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1.3. 
Determinantes 


En esta secci´on no vamos a dar la deﬁnici´on de determinante de una matriz 
cuadrada, sino que nos conformamos con saber que a cada matriz A ∈ Mn(R) se 
le puede asociar un n´umero real llamado determinante de A y que representamos 
por |A|. 


1.3.1. 
C´alculo de determinantes 


El determinante de una matriz de tama˜no 1, A = (a), es 


|A| = a 


es decir, es igual al ´unico elemento que tiene la matriz A. 


Determinantes de matrices de tama˜no 2 
Sea A ∈ M2(R) 


A = 
a11 
a12 
a21 
a22 



Entonces 
|A| = a11a22 − a12a21 


Ejemplo 1.3.1 
Dada la matriz 


A = 
2 
−1 
3 
4 



se tiene que |A| = 11. 


Determinantes de matrices de tama˜no 3. Regla de Sarrus 
Sea A ∈ M3(R) 


A = 


 


 
a11 
a12 
a13 
a21 
a22 
a23 
a31 
a32 
a33 


 


 


Entonces 


|A| = a11a22a33 + a12a23a31 + a13a21a32 − (a13a22a31 + a11a23a32 + a12a21a33) 


Ejemplo 1.3.2 
Dada la matriz 


A = 


 


 
1 
−1 
0 
2 
3 
−4 
−1 
6 
5 


 


 


se tiene que |A| = 45. 
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Determinantes de matrices de tama˜no superior a 3 
Dada una matriz A ∈ Mn(R), llamamos menor complementario del elemento 
aij, con 1 i, j n, al determinante de la submatriz de A que se obtiene al quitar, 
en la matriz A, la ﬁla i y la columna j. 


Ejemplo 1.3.3 
Dada la matriz 


A = 


 


 
 
 


3 
−2 
0 
1 
1 
3 
4 
−6 
7 
1 
−2 
−2 
1 
0 
3 
0 


 


 
 
 


el menor complementario del elemento a41 = 1 es 30. El menor complementario 
del elemento a43 = 3 es 60. 


Llamamos adjunto del elemento aij, al resultado de multiplicar (−1)i+j y el 
menor complementario de aij. 
El adjunto del elemento aij lo representamos por Aij. 


Ejemplo 1.3.4 
Dada la matriz A del Ejemplo 1.3.3, el adjunto del elemento a41 = 1 es 


A41 = (−1)4+1 · 30 = −30 


El adjunto del elemento a43 = 3 es 


A43 = (−1)4+3 · 60 = −60 



Eligiendo una ﬁla de la matriz A, por ejemplo la ﬁla i, se tiene 


|A| = 


n 


j=1 
aijAij 


El determinante de una matriz es independiente de la ﬁla elegida para calcularlo. 
Igualmente, eligiendo una columna de la matriz A, por ejemplo la columna j, 
se tiene 


|A| = 


n 


i=1 
aijAij 


El determinante de una matriz es independiente de la columna elegida para calcu- 
larlo. 


Ejemplo 1.3.5 
Dada la matriz A del Ejemplo 1.3.3, se tiene que 


|A| = a41A41 + a43A43 = 1 · (−30) + 3 · (−60) = −210. 
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1.3.2. 
Propiedades de los determinantes 


El determinante de una matriz cuadrada tiene, entre otras, las siguientes propie- 
dades: 


1. Si A ∈ Mn(R), entonces 
|At| = |A| 


Ejemplo 1.3.6 
Dada la matriz 


A = 
1 
−1 
3 
2 



se tiene que |A| = 5 = |At|. 


2. Si en una matriz A ∈ Mn(R) se permutan entre s´ı dos l´ıneas, se obtiene una 
matriz B ∈ Mn(R) tal que 
|B| = −|A| 


Ejemplo 1.3.7 
Dada la matriz A del Ejemplo 1.3.6, consideramos la matriz 


B = 
3 
2 
1 
−1 



que se obtiene permutando las dos ﬁlas de la matriz A. Entonces, se tiene que 


|B| = −5 = −|A| 



3. Si en una matriz A ∈ Mn(R) una de sus l´ıneas se expresa como suma de 
dos nuevas l´ıneas, entonces su determinante es igual a la suma de los dos 
determinantes de las matrices B y C, que se obtienen al sustituir, en la matriz 
dada, dicha l´ınea por cada una de las l´ıneas sumandos. 


Ejemplo 1.3.8 
Dadas las matrices 


A = 
5 
−1 
4 
2 


= 
3 + 2 
−1 
1 + 3 
2 


, 
B = 
3 
−1 
1 
2 


, 
C = 
2 
−1 
3 
2 



se tiene que 
|A| = 14 = 7 + 7 = |B| + |C| 


Notar que la primera columna de la matriz A es la suma de la primera columna de 
la matriz B y de la primera columna de la matriz C. 
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4. Si en una matriz A ∈ Mn(R) se multiplica una de sus l´ıneas por α ∈ R, 
entonces el determinante de la matriz B as´ı obtenida es 


|B| = α |A| 


Ejemplo 1.3.9 
Sea α = 2 y sean las matrices 


A = 
1 
−1 
3 
2 


, 
B = 
2 
−1 
6 
4 


= 
α · 1 
−1 
α · 3 
2 



Entonces 
|B| = 10 = 2 · 5 = 2|A| = α |A| 



5. Si en una matriz cuadrada una l´ınea es combinaci´on lineal de otras l´ıneas, 
entonces su determinantes es cero. Como consecuencia, si una matriz cua- 
drada tiene todos los elementos de una l´ınea nulos o tiene dos l´ıneas iguales, 
entonces su determinantes es cero. 


Ejemplo 1.3.10 
Dada la matriz 


A = 
0 
0 
1 
−1 



se tiene que |A| = 0. 
Dada la matriz 


B = 
1 
2 
1 
2 



se tiene que |B| = 0. 
Dada la matriz 


C = 


 


 
1 
−1 
2 
1 
3 
−1 
1 
−9 
8 


 


 


se tiene que |C| = 0. 
Notar que la tercera ﬁla de la matriz C es la suma de la primera ﬁla multiplicada 
por 3 y de la segunda ﬁla multiplicada por −2. Se dice, entonces, que la tercera ﬁla 
es combinaci´on lineal de las dos primeras ﬁlas. El concepto de combinaci´on lineal 
se formalizar´a al estudiar los espacios vectoriales en el Tema 3. 


6. Si a una l´ınea de una matriz cuadrada se le suma una combinaci´on lineal de 
otras l´ıneas, entonces el determinante de la matriz as´ı obtenida es igual al 
determinante de la matriz de partida. 
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Ejemplo 1.3.11 
Dada la matriz 


A = 


 


 
1 
−2 
3 
4 
−1 
2 
3 
−5 
4 


 


 


se tiene que 


|A| = 



1 
−2 
3 
4 
−1 
2 
3 
−5 
4 


= 



1 
−2 
3 
0 
7 
−10 
3 
−5 
4 


= 



1 
−2 
3 
0 
7 
−10 
0 
1 
−5 



= − 



1 
−2 
3 
0 
1 
−5 
0 
7 
−10 


= − 



1 
−2 
3 
0 
1 
−5 
0 
0 
25 


= −25 


El ´ultimo determinante se obtiene f´acilmente aplicando la regla de Sarrus. 


1.4. 
Matriz inversa. Rango de una matriz 


El concepto de inversa de una matriz es similar, en cierto sentido, al de inverso 
de un n´umero real, aunque, como veremos, no todas las matrices tienen inversa. La 
noci´on de inversa de una matriz la utilizaremos en el Tema 4 cuando introduzcamos 
el concepto de matriz diagonalizable. Adem´as, es una herramienta ´util a la hora de 
resolver sistemas de ecuaciones lineales y ecuaciones matriciales. 
El rango de una matriz veremos que es una pieza fundamental a la hora de deter- 
minar la compatibilidad de un sistema de ecuaciones lineales (Teorema de Rouch´e- 
Fr¨obenius) y en problemas sobre variedades lineales aﬁnes (rectas y planos) que 
estudiaremos en el Tema 2. 


1.4.1. 
Matriz inversa 


Una matriz cuadrada A ∈ Mn(R) es una matriz regular si existe una matriz 
B ∈ Mn(R) tal que 
A · B = I = B · A 


La matriz B es ´unica y se llama matriz inversa de A; se representa por A−1. Se 
puede demostrar que (A−1)−1 = A. 


Teorema 1.4.1 (Matriz regular) 
Una matriz A ∈ Mn(R) es regular si, y s´olo si, |A| = 0. 


Ejemplo 1.4.1 
La matriz 


A = 
1 
2 
0 
1 



Fundamentos Matem´aticos de la Ingenier´ıa – 503 – 1r. semestre 
2009/2010 
12 
c 
UJI 


2 
M. Barreda Rochera / J. A. López Ortí - ISBN: 978-84-692-9833-6 
Fundamentos Matemáticos de la Ingeniería. Parte I: Álgebra Lineal - UJI 


es una matriz regular porque |A| = 1 = 0. Adem´as, 


A−1 = 
1 
−2 
0 
1 



ya que A · A−1 = I = A−1 · A. 


Dada una matriz regular A ∈ Mn(R), se tiene que 


A−1 = 1 


|A|adj(A)t 


donde 


adj(A) = 


 


 
 
 
 


A11 
A12 
· · · 
A1n 
A21 
A22 
· · · 
A2n 
... 
... 
... 
... 
An1 
An2 
· · · 
Ann 


 


 
 
 
 


Ejercicio 1.4.1 
Calcula la matriz inversa de las siguientes matrices: 


A = 
2 
3 
−4 
1 


B = 


 


 
1 
2 
3 
0 
−1 
4 
−2 
5 
0 


 


 
C = 


 


 
1 
−3 
2 
−3 
3 
−1 
2 
−1 
0 


 


 


D = 


 


 
2 
3 
1 
1 
2 
3 
3 
1 
2 


 


 
E = 


 


 
−1 
0 
0 
0 
−1 
0 
0 
0 
−1 


 


 
F = 


 


 
−1 
2 
0 
4 
−3 
5 
0 
1 
1 


 


 



1.4.2. 
Rango de una matriz 


Dada una matriz A ∈ Mm×n(R), llamamos menor de orden p de la matriz A, 
siendo 1 p m´ın(m, n), al determinante de una submatriz de tama˜no p de la 
matriz A. 


Ejemplo 1.4.2 
Sea la matriz 


A = 


 


 
 
 


3 
−2 
0 
1 
3 
4 
4 
1 
4 
2 
−5 
−4 


 


 
 
 


Cada uno de los elementos de la matriz es un menor de orden 1. 
Eligiendo, por ejemplo, las ﬁlas 2a y 3a, y las columnas 1a y 2a, tenemos el menor 
de orden 2: 
1 
3 
4 
1 


= −11 
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Eligiendo ahora, por ejemplo, las ﬁlas 1a, 2a y 4a, y las 3 columnas, tenemos el 
menor de orden 3: 


3 
−2 
0 
1 
3 
4 
2 
−5 
−4 


= 0 



Dada una matriz A ∈ Mm×n(R), A = O, llamamos rango de la matriz A, y 
lo denotamos por rang(A), al mayor de los ´ordenes de los menores no nulos. Por 
tanto 
1 rang(A) m´ın(m, n) 


Se conviene en decir que el rango de la matriz nula es cero: 


rang(O) = 0 


Dada una matriz A ∈ Mm×n(R) y un menor de orden p de A, se llama orlar el 
menor de orden p de la matriz A a obtener un menor de orden p + 1 de A a˜nadiendo 
una nueva ﬁla y una nueva columna al menor de orden p de la matriz A. 


Dada una matriz A ∈ Mm×n(R), A = O, el c´alculo de su rango se puede 
simpliﬁcar teniendo en cuenta las siguientes propiedades: 


1. Si todos los menores de orden p de una matriz son nulos, entonces tambi´en lo 
son los de orden p + 1. 


2. Como la matriz A es no nula, deber´a tener un elemento (menor de orden 1) 
distinto de cero. Se consideran, entonces, los menores de orden 2 orlados a 
partir de este menor de orden 1 distinto de cero; si todos estos menores de 
orden 2 son iguales a cero, entonces tambi´en lo son el resto de menores de 
orden 2 y, por tanto, rang(A) = 1. 


Si hay alg´un menor de orden 2 distinto de cero, entonces se consideran los 
menores de orden 3 orlados a partir del menor de orden 2 distinto de cero; si 
todos estos menores de orden 3 son iguales a cero, entonces tambi´en lo son el 
resto de menores de orden 3 y, por tanto, rang(A) = 2. 


Si hay alg´un menor de orden 3 distinto de cero, entonces se consideran los 
menores de orden 4 orlados a partir del menor de orden 3 distinto de cero . . . 


Ejemplo 1.4.3 
El rango de la matriz del Ejemplo 1.4.2 


A = 


 


 
 
 


3 
−2 
0 
1 
3 
4 
4 
1 
4 
2 
−5 
−4 


 


 
 
 


es igual a 2, ya que tiene un menor de orden 2 distinto de cero 
1 
3 
4 
1 


= −11 = 0 
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y los menores de orden 3 orlados a partir de este menor de orden 2 son nulos: 


3 
−2 
0 
1 
3 
4 
4 
1 
4 


= 0 , 



1 
3 
4 
4 
1 
4 
2 
−5 
−4 


= 0 



Se dice que dos matrices A, B ∈ Mm×n(R) son equivalentes si 


rang(A) = rang(B) 


Se denota escribiendo A ∼ B. 


Operaciones elementales 
En una matriz se consideran las siguientes operaciones 


Permutar dos l´ıneas. 


Sustituir una l´ınea por ella misma multiplicada por un n´umero distinto de 
cero. 


Sustituir una l´ınea por ella misma m´as una combinaci´on lineal de otras l´ıneas. 


Estas operaciones se llaman operaciones elementales. 


Dada una matriz A ∈ Mm×n(R), si realizamos operaciones elementales sobre 
sus l´ıneas, entonces obtenemos otra matriz B ∈ Mm×n(R) tal que A ∼ B. 
Una matriz se dice que es una matriz escalonada si el primer elemento no nulo 
en cada ﬁla est´a m´as a la derecha que el de la ﬁla anterior. Se puede comprobar que 
el rango de una matriz escalonada coincide con el n´umero de ﬁlas no nulas. 


Ejemplo 1.4.4 
La matriz A del Ejemplo 1.4.2 tiene rango 2 porque la matriz B as´ı obtenida 


A = 


 


 
 
 


3 
−2 
0 
1 
3 
4 
4 
1 
4 
2 
−5 
−4 


 


 
 
 ∼ 


 


 
 
 


1 
3 
4 
3 
−2 
0 
4 
1 
4 
2 
−5 
−4 


 


 
 
 ∼ 


 


 
 
 


1 
3 
4 
0 
11 
12 
0 
11 
12 
0 
11 
12 


 


 
 
 ∼ 


 


 
 
 


1 
3 
4 
0 
11 
12 
0 
0 
0 
0 
0 
0 


 


 
 
 = B 


es una matriz escalonada con dos ﬁlas no nulas, equivalente a la matriz A. Es decir 


rang(A) = rang(B) = 2 



El proceso seguido en este ´ultimo ejemplo para obtener la matriz escalonada B 
a partir de la matriz A, realizando operaciones elementales sobre las ﬁlas de A, se 
llama escalonar la matriz A. 
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Ejercicio 1.4.2 
Calcula el rango de las siguientes matrices: 


A = 


 


 
1 
1 
1 
2 
−1 
3 
1 
−2 
2 


 


 
B = 
sin α 
cos α 
− cos α 
sin α 


C = 


 


 
1 
−1 
0 
2 
0 
1 
0 
−1 
3 


 


 


D = 


 


 
 
 


1 
1 
1 
1 
0 
−1 
2 
0 
0 
0 
−1 
3 
1 
2 
−1 
1 


 


 
 
 
E = 


 


 
1 
2 
0 
3 
0 
−1 
−1 
1 
1 
0 
−2 
5 


 


 
F = 


 


 
 
 
 
 


1 
2 
0 
−2 
0 
3 
−1 
2 
3 
−3 
−2 
3 
2 
4 
0 


 


 
 
 
 
 



Dada una matriz A ∈ Mn(R), llamamos menor principal de orden p de la 
matriz A a todo menor de orden p de la matriz A, cuya diagonal principal est´e con- 
tenida en la diagonal principal de A. 


Ejemplo 1.4.5 
Dada la matriz 


A = 


 


 
 
 


1 
0 
−2 
3 
4 
7 
9 
−5 
6 
−3 
8 
0 
−5 
2 
−6 
5 


 


 
 
 


se tienen los siguientes menores principales: 


Menores principales de orden 1: 1, 7, 8, 5 


Menores principales de orden 2: 
1 
0 
4 
7 


, 
1 
−2 
6 
8 


, 
1 
3 
−5 
5 


, 
7 
9 
−3 
8 


, 
7 
−5 
2 
5 


, 
8 
0 
−6 
5 



Menores principales de orden 3: 


1 
0 
−2 
4 
7 
9 
6 
−3 
8 


, 



1 
0 
−3 
4 
7 
−5 
−5 
2 
5 


, 



1 
−2 
3 
6 
8 
0 
−5 
−6 
5 


, 



7 
9 
−5 
−3 
8 
0 
2 
−6 
5 



Menores principales de orden 4: |A| 



Llamamos menor principal conducente de orden p de la matriz A, al menor 
principal de A que se obtiene eligiendo las p primeras ﬁlas y las p primeras colum- 
nas de A. 
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Ejemplo 1.4.6 
Los menores principales conducentes de la matriz del Ejemplo 1.4.5 son 


A1 = 1 , 
A2 = 
1 
0 
4 
7 


, 
A3 = 



1 
0 
−2 
4 
7 
9 
6 
−3 
8 


, 
A4 = |A| 



1.5. 
Sistemas de ecuaciones lineales 


En esta secci´on vamos a estudiar los sistemas de ecuaciones lineales desde un 
punto de vista matricial. De esta forma se tiene una herramienta mucho m´as eﬁ- 
ciente a la hora de obtener la soluci´on (soluciones) de un sistema de ecuaciones 
lineales que con los m´etodos estudiados en el Bachillerato: igualaci´on, sustituci´on 
y reducci´on. 
Los sistemas de ecuaciones lineales son una herramienta imprescindible a la 
hora de abordar el estudio de los tres siguientes temas: geometr´ıa del espacio tridi- 
mensional, espacios vectoriales y diagonalizaci´on de matrices. 


1.5.1. 
Deﬁniciones 


Un sistema de m ecuaciones lineales con n inc´ognitas es un conjunto de ex- 
presiones de la forma: 
 
 
 
 
 


 
 
 
 


a11x1 
+ 
a12x2 
+ 
· · · 
+ 
a1nxn 
= 
b1 
a21x1 
+ 
a22x2 
+ 
· · · 
+ 
a2nxn 
= 
b2 
... 
... 
... 
... 
... 
... 
... 
... 
... 
am1x1 
+ 
am2x2 
+ 
· · · 
+ 
amnxn 
= 
bm 


(1.2) 


donde 
aij ∈ R , 
i = 1, 2, . . . , m , j = 1, 2, . . . , n 


bi ∈ R , 
i = 1, 2, . . . , m 


son valores conocidos, y 
xj , 
j = 1, 2, . . . , n 


son las inc´ognitas del sistema. 


Ejemplo 1.5.1 
Sistema de 2 ecuaciones lineales y 3 inc´ognitas: 
2x 
− 
3y 
+ 
z 
= 
5 
−5x 
+ 
7y 
+ 
6z 
= 
23 
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Una soluci´on del sistema lineal (1.2), es un vector1 (x∗ 
1, x∗ 
2, . . . , x∗ 
n) ∈ Rn tal 
que 
a11x∗ 
1 + a12x∗ 
2 + · · · + a1nx∗ 
n = b1 
a21x∗ 
1 + a22x∗ 
2 + · · · + a2nx∗ 
n = b2 
... 
... 
... 
... 
... 
... 
... 
... 
... 
am1x∗ 
1 + am2x∗ 
2 + · · · + amnx∗ 
n = bm 


Ejemplo 1.5.2 
El vector (−4, −3, 4) es una soluci´on del sistema del Ejemplo 1.5.1, ya que 


2 · (−4) − 3 · (−3) + 1 · 4 = 5 
−5 · (−4) + 7 · (−3) + 6 · 4 = 23 



Ejercicio 1.5.1 
Demuestra que, para cada z ∈ R, el vector (−104 + 25z, −71 + 17z, z) es una 
soluci´on del sistema del Ejemplo 1.5.1. 


Dos sistemas de ecuaciones lineales son equivalentes si tienen las mismas 
soluciones. 


El sistema (1.2) se puede escribir en forma matricial como 


A · X = B 
(1.3) 


donde 


A = 


 


 
 
 
 


a11 
a12 
· · · 
a1n 
a21 
a22 
· · · 
a2n 
... 
... 
... 
... 
am1 
am2 
· · · 
amn 


 


 
 
 
 


es la matriz de coeﬁcientes del sistema, y 


X = 


 


 
 
 
 


x1 
x2 
... 
xn 


 


 
 
 
 , 
B = 


 


 
 
 
 


b1 
b2 
... 
bm 


 


 
 
 
 


son, respectivamente, la matriz inc´ognita y la matriz de los t´erminos indepen- 
dientes. 


V´ease aqu´ı la analog´ıa se˜nalada en la introducci´on de este tema entre las ecua- 
ciones (1.3) y (1.1). 


1El concepto de vector y el de espacio vectorial Rn lo estudiaremos en el siguiente tema. 
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La matriz 


AB = 


 


 
 
 
 


a11 
a12 
· · · 
a1n 
b1 
a21 
a22 
· · · 
a2n 
b2 
... 
... 
... 
... 
... 
am1 
am2 
· · · 
amn 
bm 


 


 
 
 
 


es la matriz ampliada del sistema. No confundir la matriz ampliada AB con el 
producto de la matrices A y B, A · B, el cual s´olo tendr´ıa sentido si m = n. 


Ejemplo 1.5.3 
Dado el sistema de ecuaciones lineales del Ejemplo 1.5.1, se tiene que la matriz de 
coeﬁcientes y la matriz ampliada del sistema son, respectivamente 


A = 
2 
−3 
1 
−5 
7 
6 


, 
AB = 
2 
−3 
1 
5 
−5 
7 
6 
23 



La matriz de los t´erminos independientes y la matriz inc´ognita son, respectivamente 


B = 
5 
23 


, 
X = 


 


 
x 
y 
z 


 


 



Si realizamos operaciones elementales sobre las ﬁlas de la matriz ampliada del 
sistema (1.3), AB, entonces obtenemos otra matriz 


A′B′ = 


 


 
 
 
 


a′ 
11 
a′ 
12 
· · · 
a′ 
1n 
b′ 
1 
a′ 
21 
a′ 
22 
· · · 
a′ 
2n 
b′ 
2 
... 
... 
... 
... 
... 
a′ 
m1 
a′ 
m2 
· · · 
a′ 
mn 
b′ 
m 


 


 
 
 
 


tal que, el sistema 
A′ · X = B′ 
(1.4) 


siendo 


A′ = 


 


 
 
 
 


a′ 
11 
a′ 
12 
· · · 
a′ 
1n 
a′ 
21 
a′ 
22 
· · · 
a′ 
2n 
... 
... 
... 
... 
a′ 
m1 
a′ 
m2 
· · · 
a′ 
mn 


 


 
 
 
 , 
B′ = 


 


 
 
 
 


b′ 
1 
b′ 
2... 
b′ 
m 


 


 
 
 
 


las matrices de coeﬁcientes y ampliada, respectivamente, del sistema (1.4), es equi- 
valente al sistema (1.3). 


Ejercicio 1.5.2 
Comprueba que el sistema de ecuaciones lineales dado en el Ejemplo 1.5.1 y el 
sistema 
2x 
− 
3y 
+ 
z 
= 
5 
y 
− 
17z 
= 
−71 


son equivalentes. 
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1.5.2. 
Clasiﬁcaci´on de los sistemas de ecuaciones lineales 


Un sistema es incompatible si no tiene soluci´on; un sistema es compatible 
determinado si s´olo tiene una soluci´on; un sistema es compatible indeterminado 
si tiene m´as de una soluci´on (inﬁnitas). 


Un sistema es homog´eneo si B = 0, y es inhomog´eneo si B = 0. 


Un sistema homog´eneo siempre tiene como soluci´on 


x1 = 0, x2 = 0, . . . , xn = 0 


que se llama soluci´on trivial. 


Ejemplo 1.5.4 
El sistema de ecuaciones lineales dado en el Ejemplo 1.5.1 es compatible indeter- 
minado, ya que tiene inﬁnitas soluciones: 


(−104 + 25z, −71 + 17z, z) , 
z ∈ R 


y es inhomog´eneo. 


Teorema 1.5.1 (de Rouch´e-Fr¨obenius) 
Dado el sistema de m ecuaciones lineales y n inc´ognitas 


A · X = B 


se tiene: 


1. B = 0 


a) Si rang(A) < rang(AB), entonces el sistema es incompatible. 


b) Si rang(A) = rang(AB) < n, entonces el sistema es compatible inde- 
terminado. 


c) Si rang(A) = rang(AB) = n, entonces el sistema es compatible deter- 
minado. 


2. B = 0 


a) Si rang(A) < n, entonces el sistema es compatible indeterminado. 


b) Si rang(A) = n, entonces el sistema s´olo tiene la soluci´on trivial. 



Ejemplo 1.5.5 
El sistema de ecuaciones lineales 
 
 


 


x 
− 
y 
+ 
z 
= 
4 
−x 
+ 
2y 
− 
z 
= 
6 
2x 
+ 
z 
= 
0 
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es compatible determinado porque 


AB = 


 


 
1 
−1 
1 
4 
−1 
2 
−1 
6 
2 
0 
1 
0 


 


 ∼ 


 


 
1 
−1 
1 
4 
0 
1 
0 
10 
0 
−2 
1 
8 


 


 ∼ 


 


 
1 
−1 
1 
4 
0 
1 
0 
10 
0 
0 
1 
28 


 


 


y as´ı 
rang(A) = 3 = rang(AB) = n´umero de inc´ognitas 


Notar que en el proceso de escalonar la matriz AB tambi´en se obtiene que 


A = 


 


 
1 
−1 
1 
−1 
2 
−1 
2 
0 
1 


 


 ∼ 


 


 
1 
−1 
1 
0 
1 
0 
0 
−2 
1 


 


 ∼ 


 


 
1 
−1 
1 
0 
1 
0 
0 
0 
1 


 


 



Ejemplo 1.5.6 
El sistema de ecuaciones lineales 
 
 


 


3x 
− 
y 
+ 
4z 
= 
5 
2x 
− 
2y 
+ 
z 
= 
1 
−4x 
+ 
8y 
+ 
3z 
= 
5 


es compatible indeterminado porque 


AB = 


 


 
3 
−1 
4 
5 
2 
−2 
1 
1 
−4 
8 
3 
5 


 


 ∼ 


 


 
3 
−1 
4 
5 
0 
4 
5 
7 
0 
20 
25 
35 


 


 ∼ 


 


 
3 
−1 
4 
5 
0 
4 
5 
7 
0 
0 
0 
0 


 


 


y as´ı 
rang(A) = 2 = rang(AB) < 3 = n´umero de inc´ognitas 



Ejemplo 1.5.7 
El sistema de ecuaciones lineales 
 
 


 


x 
+ 
y 
= 
−1 
2x 
− 
y 
= 
7 
x 
+ 
2y 
= 
4 


es incompatible porque 


AB = 


 


 
1 
1 
−1 
2 
−1 
7 
1 
2 
4 


 


 ∼ 


 


 
1 
1 
−1 
0 
3 
−9 
0 
−1 
−5 


 


 ∼ 


 


 
1 
1 
−1 
0 
3 
−9 
0 
0 
−24 


 


 


y as´ı 
rang(A) = 2 < 3 = rang(ABB) 
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Deﬁnici´on 1.5.1 (Sistema de Cramer) 
Un sistema de m ecuaciones lineales y n inc´ognitas 


A · X = B 


es un sistema de Cramer, si 


1. Tiene el mismo n´umero de ecuaciones que de inc´ognitas (m = n). 


2. La matriz de coeﬁcientes del sistema, A, es regular, es decir |A| = 0. 



Por tanto, todo sistema de Cramer es compatible determinado. 


Ejemplo 1.5.8 
El sistema de ecuaciones lineales dado en el Ejemplo 1.5.5 es un sistema de Cramer, 
ya que tiene 3 ecuaciones, 3 inc´ognitas y rang(A) = 3. 


Teorema 1.5.2 (F´ormula de Cramer) 
La soluci´on del sistema de Cramer A · X = B, donde 


A = 


 


 
 
 
 


a11 
a12 
· · · 
a1n 
a21 
a22 
· · · 
a2n 
... 
... 
... 
... 
an1 
an2 
· · · 
ann 


 


 
 
 
 , 
B = 


 


 
 
 
 


b1 
b2 
... 
bn 


 


 
 
 
 


es 


xi = |Mi| 


|A| , 
i = 1, 2, . . . , n 


siendo 


Mi = 


 


 
 
 
 


a11 
a12 
· · · 
a1 i−1 
b1 
a1 i+1 
· · · 
a1n 
a21 
a22 
· · · 
a2 i−1 
b2 
a2 i+1 
· · · 
a2n 
... 
... 
... 
... 
... 
... 
... 
... 
an1 
an2 
· · · 
an i−1 
bn 
an i+1 
· · · 
ann 


 


 
 
 
 



Ejemplo 1.5.9 
La soluci´on del sistema de Cramer del Ejemplo 1.5.5 es: 


x = 



4 
−1 
1 
6 
2 
−1 
0 
0 
1 




1 
−1 
1 
−1 
2 
−1 
2 
0 
1 



= 14 


−1 = −14 
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y = 



1 
4 
1 
−1 
6 
−1 
2 
0 
1 




1 
−1 
1 
−1 
2 
−1 
2 
0 
1 



= −10 


−1 = 10 


z = 



1 
−1 
4 
−1 
2 
6 
2 
0 
0 




1 
−1 
1 
−1 
2 
−1 
2 
0 
1 



= −28 


−1 = 28 



Cuando un sistema es compatible indeterminado, se consideran las ecuaciones 
que forman parte del menor de mayor orden no nulo de la matriz de coeﬁcientes y 
se toman como inc´ognitas las correspondientes a las columnas de dicho menor no 
nulo, mientras que las dem´as inc´ognitas se consideran como par´ametros. Este nuevo 
sistema en el que aparecen par´ametros es, entonces, de Cramer, tal como se muestra 
en el siguiente 


Ejemplo 1.5.10 
El sistema de ecuaciones lineales del Ejemplo 1.5.6 
 
 


 


3x 
− 
y 
+ 
4z 
= 
5 
2x 
− 
2y 
+ 
z 
= 
1 
−4x 
+ 
8y 
+ 
3z 
= 
5 


hemos visto que es compatible indeterminado, ya que 


rang(A) = 2 = rang(AB) < 3 = n´umero de inc´ognitas 


Como 
3 
−1 
2 
−2 


= −4 = 0 


es un menor de orden 2 distinto de cero, ese sistema se puede transformar en el 
sistema de Cramer 
3x 
− 
y 
= 
5 − 4z 
2x 
− 
2y 
= 
1 − z 


de 2 ecuaciones, 2 inc´ognitas (x e y) y 1 par´ametro (z). 
La soluci´on del sistema, que depende del par´ametro z, es: 


x = 


5 − 4z 
−1 
1 − z 
−2 



3 
−1 
2 
−2 



= −9 + 7z 


−4 
= 9 


4 − 7 


4z 
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y = 


3 
5 − 4z 
2 
1 − z 



3 
−1 
2 
−2 



= −7 + 5z 


−4 
= 7 


4 − 5 


4z 


Por tanto, las inﬁnitas soluciones del sistema original son de la forma 


(9/4 − 7/4z, 7/4 − 5/4z, z) , 
z ∈ R 


o equivalentemente 


(4 − 7λ, 3 − 5λ, −1 + 4λ) , 
λ ∈ R 



Ejemplo 1.5.11 
Discute, en funci´on del par´ametro real a, la compatibilidad del sistema de ecuacio- 
nes lineales 
 
 


 


a y 
− 
z 
= 
a 
a x 
+ 
2z 
= 
0 
x 
+ 
3y 
+ 
z 
= 
5 


Resu´elvelo cuando sea posible. 
Las matrices de coeﬁcientes y ampliada del sistema son: 


A = 


 


 
0 
a 
−1 
a 
0 
2 
1 
3 
1 


 


 , 
AB = 


 


 
0 
a 
−1 
a 
a 
0 
2 
0 
1 
3 
1 
5 


 


 


Calculamos el determinante de la matriz de coeﬁcientes: 


|A| = 



0 
a 
−1 
a 
0 
2 
1 
3 
1 


= 2a − 3a − a2 = −a − a2 = −a(a + 1) 


Por tanto 
|A| = 0 ⇐⇒ a(a + 1) = 0 ⇐⇒ a = 0 
´o 
a = −1 


Esto da lugar a los tres siguientes casos: 


1 a = −1, 0 


En este caso, |A| = 0. De aqu´ı 


rang(A) = 3 = rang(AB) 


y como el n´umero de inc´ognitas tambi´en es igual a 3, obtenemos que el siste- 
ma es compatible determinado. 
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Resolviendo el sistema por el m´etodo de Cramer, se tiene 


x = 



a 
a 
−1 
0 
0 
2 
5 
3 
1 




0 
a 
−1 
a 
0 
2 
1 
3 
1 



= 10a − 6a 


−a(a + 1) = 
4a 


−a(a + 1) = − 
4 


a + 1 


y = 



0 
a 
−1 
a 
0 
2 
1 
5 
1 




0 
a 
−1 
a 
0 
2 
1 
3 
1 



= 2a − 5a − a2 


−a(a + 1) 
= −3a − a2 


−a(a + 1) = −a(a + 3) 


−a(a + 1) = a + 3 


a + 1 


z = 



0 
a 
a 
a 
0 
0 
1 
3 
5 




0 
a 
−1 
a 
0 
2 
1 
3 
1 



= 3a2 − 5a2 


−a(a + 1) = 
−2a2 


−a(a + 1) = 
2a 


a + 1 


2 a = −1. 


Escalonamos la matriz de coeﬁcientes y la matriz ampliada del sistema: 


AB = 


 


 
0 
−1 
−1 
−1 
−1 
0 
2 
0 
1 
3 
1 
5 


 


 ∼ 


 


 
1 
3 
1 
5 
0 
−1 
−1 
−1 
−1 
0 
2 
0 


 


 


∼ 


 


 
1 
3 
1 
5 
0 
−1 
−1 
−1 
0 
3 
3 
5 


 


 ∼ 


 


 
1 
3 
1 
5 
0 
−1 
−1 
−1 
0 
0 
0 
2 


 


 


De aqu´ı 
rang(A) = 2 < 3 = rang(AB) 


Por tanto, el sistema es incompatible. 


3 a = 0. 


Escalonamos la matriz de coeﬁcientes y la matriz ampliada del sistema: 


AB = 


 


 
0 
0 
−1 
0 
0 
0 
2 
0 
1 
3 
1 
5 


 


 ∼ 


 


 
1 
3 
1 
5 
0 
0 
1 
0 
0 
0 
2 
0 


 


 ∼ 


 


 
1 
3 
1 
5 
0 
0 
1 
0 
0 
0 
0 
0 


 


 


De aqu´ı 
rang(A) = 2 = rang(AB) 
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y como el n´umero de inc´ognitas es igual a 3, obtenemos que el sistema es 
compatible indeterminado. 


Por ´ultimo, resolvemos el sistema compatible indeterminado: 
x 
+ 
3y 
+ 
z 
= 
5 
z 
= 
0 


Sustituyendo la segunda ecuaci´on, z = 0, en la primera ecuaci´on de este 
sistema, obtenemos: 
x = 5 − 3y 


es decir, las inﬁnitas soluciones del sistema son los vectores (5 − 3y, y, 0), 
donde y ∈ R. Dicho de otra forma, las soluciones del sistema son los elemen- 
tos del conjunto 
{(5 − 3y, y, 0) : y ∈ R} 


Tal y como veremos en el siguiente tema, los elementos de este conjunto 
forman una recta en el espacio tridimensional. 



Ejercicio 1.5.3 
Discute, en funci´on de los par´ametros reales, la compatibilidad de los siguientes 
sistemas de ecuaciones lineales. Resu´elvelos cuando sea posible. 


(a) 


 
 
 
 


 
 
 


x 
+ 
2y 
+ 
z 
= 
2 
x 
+ 
y 
+ 
2z 
= 
3 
x 
+ 
3y 
+ 
a z 
= 
1 
x 
+ 
2y 
+ 
z 
= 
b 


(b) 


 
 


 


x 
+ 
y 
= 
2 
a x 
+ 
y 
= 
1 
x 
− 
y 
= 
a 


(c) 


 
 
 
 


 
 
 


a x 
+ 
y 
+ 
z 
= 
a 
x 
+ 
a y 
− 
z 
= 
1 
3x 
+ 
y 
+ 
b z 
= 
2 
x 
− 
y 
− 
z 
= 
1 


(d) 


 
 


 


(a + 1) x 
+ 
y 
+ 
z 
= 
0 
x 
+ 
(a + 1) y 
+ 
z 
= 
0 
x 
+ 
y 
+ 
(a + 1) z 
= 
0 


(e) 


 
 


 


x 
+ 
2y 
+ 
a z 
= 
6 
x 
+ 
(a − 1) y 
− 
z 
= 
−2 
2x 
− 
a z 
= 
−1 


(f ) 


 
 


 


x 
+ 
a y 
− 
z 
= 
−4 
a x 
+ 
z 
= 
8 
−x 
+ 
2y 
− 
a z 
= 
−a − 10 
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(g) 


 
 
 
 


 
 
 


x 
+ 
y 
+ 
z 
+ 
t 
= 
10 
x 
− 
z 
+ 
t 
= 
2 
y 
+ 
z 
= 
a 
2x 
+ 
y 
+ 
2t 
= 
0 


(h) 


 
 
 
 


 
 
 


2a x 
− 
y 
+ 
z 
= 
0 
x 
− 
2a z 
= 
−1 
− 
a y 
+ 
5z 
= 
2 
3a x 
− 
y 
− 
z 
= 
−a 


(i) 


 
 


 


2x 
+ 
2y 
+ 
a z 
= 
1 
−x 
+ 
z 
= 
a 
x 
+ 
2a y 
+ 
2z 
= 
2 
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Tema 2 


Geometr´ıa del plano y del espacio 


2.1. 
Introducci´on 


En este tema estudiamos la geometr´ıa del plano y del espacio, centr´andonos en 
el espacio tridimensional. 
Empezamos introduciendo los espacios vectoriales R2 y R3, dot´andolos de es- 
tructura eucl´ıdea, y a continuaci´on presentamos el espacio af´ın tridimensional con 
los problemas de incidencia y paralelismo entre rectas y planos. 
A lo largo de este tema aparecer´an conceptos que introduciremos con m´as de- 
talle en el siguiente tema. Creemos conveniente hacerlo de esta forma con el ﬁn de 
que el estudiante vaya asimilando paulatinamente el contenido de esta primera parte 
de la asignatura; por esta raz´on la exposici´on ir´a de lo particular a lo general. 


2.2. 
Los espacios vectoriales R2 y R3 


Consideramos el cuerpo de los n´umeros reales, R, a cuyos elementos llamare- 
mos escalares. El concepto de cuerpo lo introduciremos en el tema 3. 


2.2.1. 
El espacio vectorial R2 


En el conjunto R2, formado por todos los pares ordenados de n´umeros reales, 


R2 = {(x, y) : x, y ∈ R} 


se pueden deﬁnir dos operaciones: suma y producto por un escalar. 
Suma 
Dados (u1, u2), (v1, v2) ∈ R2: 


(u1, u2) + (v1, v2) = (u1 + v1, u2 + v2) 


Con esta operaci´on el conjunto R2 tiene estructura de grupo conmutativo. Esta es- 
tructura algebraica la introduciremos en el tema 3. 
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Ejemplo 2.2.1 
Sean (2, −1), (−3, 4) ∈ R2. Entonces 


(2, −1) + (−3, 4) = (−1, 3) 



Producto por un escalar 
Dados λ ∈ R y (v1, v2) ∈ R2: 


λ · (v1, v2) = (λ v1, λ v2) 


Ejemplo 2.2.2 
Sean 2 ∈ R y (−1, 4) ∈ R2. Entonces 


2 · (−1, 4) = (−2, 8) 



El producto por un escalar cumple las siguientes propiedades: 


λ · 
(u1, u2) + (v1, v2) 
= λ · (u1, u2) + λ · (v1, v2), para todo λ ∈ R y para 
todo (u1, u2), (v1, v2) ∈ R2. 


(λ + µ) · (v1, v2) = λ · (v1, v2) + µ · (v1, v2), para todo λ, µ ∈ R y para todo 
(v1, v2) ∈ R2. 


λ· 
µ·(v1, v2) 
= (λ µ)·(v1, v2), para todo λ, µ ∈ R y para todo (v1, v2) ∈ R2. 


1 · (v1, v2) = (v1, v2), para todo (v1, v2) ∈ R2. 


Con estas dos operaciones el conjunto R2 tiene estructura de espacio vectorial so- 
bre el cuerpo R. Esta estructura algebraica tambi´en ser´a introducida en el siguiente 
tema. 
Los elementos de R2 se llaman vectores y los denotamos por 


u = (u1, u2) , 
v = (v1, v2) , 
. . . 


El vector nulo lo denotamos por 0 = (0, 0). 


El espacio vectorial R2 lo representamos en el plano considerando dos ejes per- 
pendiculares. 
Los vectores u = (3, 2),v = (−2, 3) ∈ R2 los representamos as´ı: 


1 
2 
3 
−1 
−2 
−3 


1 


2 


3 


−1 


u 
v 


X 


Y 
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La suma de vectores de R2 cumple la ley del paralelogramo. As´ı, dados dos 
vectores u y v de R2, la suma de los vectores u y v, u + v, es la diagonal del 
paralelogramo que determinan u y v: 


u 


v 


u + v 


X 


Y 


Un sistema de vectores de R2 es una colecci´on de vectores de R2, los cuales se 
pueden repetir. 


Ejemplo 2.2.3 
El conjunto S = {(−1, 2), (0, −3), (5, 1)} es un sistema de vectores del espacio 
vectorial R2. 


Una base del espacio vectorial R2 es un sistema de vectores de R2, formado por 
dos vectores {u,v} tal que 
u1 
u2 
v1 
v2 


= 0 


siendo u = (u1, u2) y v = (v1, v2). 


Ejemplo 2.2.4 
El sistema de vectores B = {(1, −1), (2, 3)} de R2 es una base del espacio vectorial 
R2, ya que 
1 
−1 
2 
3 


= 5 = 0 



El sistema de vectores {ı, }, dondeı = (1, 0) y  = (0, 1), es una base de R2, 
llamada base can´onica de R2. 


Se puede demostrar que dada una base B = {e1,e2} de R2 y un vector v ∈ R2, 
existen dos ´unicos escalares λ1, λ2 ∈ R tal que 


v = λ1 · e1 + λ2 · e2 


El par (λ1, λ2) se dice que son las coordenadas del vector v en la base B (respecto 
a la base B). 
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Ejemplo 2.2.5 
Las coordenadas del vector v = (−2, −13) en la base B = {(1, −1), (2, 3)} del 
espacio vectorial R2 son (4, −3), ya que 


v = (−2, −13) = 4 · (1, −1) + (−3) · (2, 3) 



Las coordenadas de un vector en la base can´onica de R2 coinciden son sus com- 
ponentes: 
v = (v1, v2) = v1 · (1, 0) + v2 · (0, 1) = v1 ·ı + v2 ·  


Ejercicio 2.2.1 
Calcula las coordenadas del vector (−3, 1) respecto a la base B = {(1, 2), (1, 3)} 
del espacio vectorial R2. 


2.2.2. 
El espacio vectorial R3 


En el conjunto R3, formado por todas las ternas ordenadas de n´umeros reales, 


R3 = {(x, y, z) : x, y, z ∈ R} 


se pueden deﬁnir dos operaciones: suma y producto por un escalar. 
Suma 
Dados (u1, u2, u3), (v1, v2, v3) ∈ R3: 


(u1, u2, u3) + (v1, v2, v3) = (u1 + v1, u2 + v2, u3 + v3) 


Con esta operaci´on el conjunto R3 tiene estructura de grupo conmutativo. 


Ejemplo 2.2.6 
Sean (2, −1, 3), (1, 4, −2) ∈ R3. Entonces 


(2, −1, 3) + (1, 4, −2) = (3, 3, 1) 



Producto por un escalar 
Dados λ ∈ R y (v1, v2, v3) ∈ R3: 


λ · (v1, v2, v3) = (λ v1, λ v2, λ v3) 


Ejemplo 2.2.7 
Sean 3 ∈ R y (2, −1, 3) ∈ R3. Entonces 


3 · (2, −1, 3) = (6, −3, 9) 
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El producto por un escalar cumple las mismas propiedades que en el caso del 
espacio vectorial R2: 


λ · 
(u1, u2, u3) + (v1, v2, v3) 
= λ · (u1, u2, u3) + λ · (v1, v2, v3), para todo 
λ ∈ R y para todo (u1, u2, u3), (v1, v2, v3) ∈ R3. 


(λ + µ) · (v1, v2, v3) = λ · (v1, v2, v3) + µ · (v1, v2, v3), para todo λ, µ ∈ R y 
para todo (v1, v2, v3) ∈ R3. 


λ· 
µ · (v1, v2, v3) 
= (λ µ) · (v1, v2, v3), para todo λ, µ ∈ R y para todo 
(v1, v2, v3) ∈ R3. 


1 · (v1, v2, v3) = (v1, v2, v3), para todo (v1, v2, v3) ∈ R3. 


Con estas dos operaciones el conjunto R3 tiene estructura de espacio vectorial 
sobre el cuerpo R. 


Los elementos de R3 se llaman vectores y los denotamos por 


u = (u1, u2, u3) , 
v = (v1, v2, v3) , 
. . . 


El vector nulo lo denotamos por 0 = (0, 0, 0). 
El espacio vectorial R3 lo representamos en el espacio considerando tres ejes 
perpendiculares. 
El vector v = (2, 3, 4) ∈ R3 lo representamos as´ı: 


X 
Y 


Z 


1.0 


2.0 


3.0 


1.0 
2.0 
3.0 


1.0 


2.0 


3.0 


4.0 


v 


La suma de vectores de R3 tambi´en cumple la regla del paralelogramo. 
Un sistema de vectores de R3 es una colecci´on de vectores de R3, los cuales se 
pueden repetir. 


Ejemplo 2.2.8 
El conjunto S = {(2, 1, 3), (4, −5, 1)} es un sistema de vectores del espacio vecto- 
rial R3. 
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Dos vectores u = (u1, u2, u3) y v = (v1, v2, v3) son linealmente independien- 
tes si 


rang 
u1 
u2 
u3 
v1 
v2 
v3 


= 2 


En otro caso se dice que los dos vectores son linealmente dependientes o que los 
dos vectores tienen la misma direcci´on. 
Si los vectores u = (u1, u2, u3) y v = (v1, v2, v3) son linealmente dependientes 
y no nulos, entonces existe un escalar λ = 0 tal que u = λ · v, por tanto 


(u1, u2, u3) = λ · (v1, v2, v3) 


es decir, si v1 = 0, v2 = 0 y v3 = 0, entonces 
u1 
v1 
= u2 


v2 
= u3 


v3 
= λ 


Si por el contrario, existe alg´un i ∈ {1, 2, 3} tal que vi = 0, entonces ui = 0. 


Ejemplo 2.2.9 
Los vectores (1, 1, 2) y (2, 3, −1) son linealmente independientes, ya que 


rang 
1 
1 
2 
2 
3 
−1 


= 2 


mientras que los vectores (2, 0, 4) y (6, 0, 12) son linealmente dependientes, ya que 


rang 
2 
0 
4 
6 
0 
12 


= 1 = 2 


En este caso 
(6, 0, 12) = 3 · (2, 0, 4) 


o lo que es lo mismo 
6 
2 = 12 


4 = 3 



Tres vectores u = (u1, u2, u3), v = (v1, v2, v3) y w = (w1, w2, w3) son lineal- 
mente independientes si 


u1 
u2 
u3 
v1 
v2 
v3 
w1 
w2 
w3 


= 0 


Se dice entonces que los vectores u, v y w son no coplanarios. 
En otro caso se dice que los tres vectores son linealmente dependientes o co- 
planarios. 


Ejemplo 2.2.10 
Los vectores (1, 2, 1), (2, 1, 0) y (3, 1, 1) son linealmente independientes, ya que 


1 
2 
1 
2 
1 
0 
3 
1 
1 


= −4 = 0 
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mientras que los vectores (1, 2, 1), (2, 1, 0) y (3, 3, 1) son linealmente dependientes, 
ya que 


1 
2 
1 
2 
1 
0 
3 
3 
1 


= 0 



Una base del espacio vectorial R3 es un sistema de vectores de R3 formado por 
tres vectores linealmente independientes. 


Ejemplo 2.2.11 
El sistema de vectores B = {(1, 2, 1), (2, 1, 0), (3, 1, 1)} del espacio vectorial R3 es 
una base de R3 ya que, como hemos visto en el Ejemplo 2.2.10, sus tres vectores 
son linealmente independientes. 


El sistema de vectores {ı,,k}, dondeı = (1, 0, 0),  = (0, 1, 0) y k = (0, 0, 1), 
es una base de R3, llamada base can´onica de R3. Hay que hacer notar aqu´ı que los 
s´ımbolos ı y  tambi´en se han utilizado anteriormente para los vectores de la base 
can´onica de R2. 


Dada una base B = {e1,e2,e3} de R3 y un vector v ∈ R3, existen tres ´unicos 
escalares λ1, λ2, λ3 ∈ R tal que 


v = λ1 · e1 + λ2 · e2 + λ3 · e3 


La terna (λ1, λ2, λ3) se dice que son las coordenadas del vector v en la base B. 


Ejemplo 2.2.12 
Las coordenadas del vector v = (5, 9, 3) respecto a la base 


B = {(1, 2, 1), (2, 1, 0), (3, 1, 1)} 


del espacio vectorial R3 son (4, 2, −1), ya que 


v = (5, 9, 3) = 4 · (1, 2, 1) + 2 · (2, 1, 0) + (−1) · (3, 1, 1) 



Ejercicio 2.2.2 
Calcula las coordenadas del vector (9, −3, 5) respecto a la base 


B = {(1, −1, 2), (0, 2, −1), (2, 1, 0)} 


del espacio vectorial R3. 


Las coordenadas de un vector en la base can´onica de R3 coinciden son sus com- 
ponentes: 


v = (v1, v2, v3) = v1 · (1, 0, 0) + v2 · (0, 1, 0) + v3 · (0, 0, 1) 


= v1 ·ı + v2 ·  + v3 · k 
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2.3. 
Los espacios vectoriales eucl´ıdeos R2 y R3 


En esta secci´on introducimos en los espacios vectoriales R2 y R3 el producto 
escalar usual, el cual nos permite deﬁnir los conceptos de m´odulo (longitud) de un 
vector y ´angulo entre dos vectores. 


2.3.1. 
Producto escalar 


Llamamos producto escalar de los vectores u = (u1, u2, u3) y v = (v1, v2, v3) 
de R3, al n´umero real 
u · v = u1v1 + u2v2 + u3v3 


Ejemplo 2.3.1 
El producto escalar de los vectores u = (2, 1, 3) y v = (−1, 1, 5) de R3 es 


u · v = 2 · (−1) + 1 · 1 + 3 · 5 = −2 + 1 + 15 = 14 



Hay que hacer notar que el s´ımbolo “ · ” tanto lo empleamos para designar el 
producto de n´umeros reales, como para el producto de un escalar y un vector o para 
el producto escalar de vectores. 


El producto escalar cumple las propiedades: 


u · v = v · u, para todo u,v ∈ R3. 


(λu + µv) · w = λu · w + µv · w, para todo u,v, w ∈ R3 y λ ∈ R. 


v · v > 0, para todo v ∈ R3, v = 0. 


De esta forma el espacio vectorial R3 es un espacio vectorial eucl´ıdeo. 


De forma similar podemos deﬁnir el espacio vectorial eucl´ıdeo R2. 


2.3.2. 
M´odulo de un vector 


Llamamos m´odulo del vector v = (v1, v2, v3) ∈ R3, al n´umero real 


||v|| = 
√ 


v · v = 


(v1)2 + (v2)2 + (v3)2 


Ejemplo 2.3.2 
El m´odulo del vector v = (2, 3, −1) es 


||v|| = 


22 + 32 + (−1)2 = 
√ 


4 + 9 + 1 = 
√ 


14 



El m´odulo de un vector cumple las siguientes propiedades: 
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||v|| > 0, para todo v ∈ R3, v = 0. Adem´as, ||0|| = 0. 


||λv|| = |λ| ||v||, para todo v ∈ R3 y λ ∈ R. 


|u · v| ||u|| ||v||, para todo u,v ∈ R3. 


||u + v|| ||u|| + ||v||, para todo u,v ∈ R3. 


De la misma forma podemos deﬁnir el m´odulo de un vector de R2. 
Un vector v ∈ R3 es un vector unitario si ||v|| = 1. 


Ejemplo 2.3.3 
El vector v = (1/ 
√ 


14, 2/ 
√ 


14, 3/ 
√ 


14) es un vector unitario, ya que 


||v|| = 


(1/ 
√ 


14)2 + (2/ 
√ 


14)2 + (3/ 
√ 


14)2 


= 


1/14 + 4/14 + 9/14 = 


14/14 = 1 



Dado un vector no nulo v ∈ R3, entonces los vectores 


v 


||v|| , 
− v 


||v|| 


son unitarios y tienen la misma direcci´on que v. 


Los vectores de la base can´onica de R3 


ı = (1, 0, 0) , 
 = (0, 1, 0) , 
k = (0, 0, 1) 


son vectores unitarios. 


Igualmente podemos deﬁnir el concepto de vector unitario en R2 y observar que 
los vectores de la base can´onica de R2 son vectores unitarios. 


2.3.3. 
´Angulo entre dos vectores 


Sean dos vectores no nulos u,v ∈ R3. El ´angulo que forman los vectores u y v 
viene dado por su coseno: 


cos α = 
u · v 


||u|| ||v|| 


donde α ∈ [0, π]. 


α 


u 


v 
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Ejemplo 2.3.4 
El ´angulo entre los vectores u = (0, −3, 5) y v = (1, 1, 10) es α ≈ 0.65 radianes, 
ya que 


cos α = 
u · v 


||u|| ||v|| = 
(0, −3, 5) · (1, 1, 10) 


||(0, −3, 5)|| ||(1, 1, 10)|| = 
47 
√ 


34 
√ 


102 = 
47 


34 
√ 


3 ≈ 0.8 


y as´ı 
α = arc cos 
47/(34 
√ 


3) 
≈ 0.65 radianes 



Dos vectores no nulos u,v ∈ R3 son ortogonales, si u · v = 0. 


Ejemplo 2.3.5 
Los vectores u = (2, −1, 3) y v = (1, 5, 1) son ortogonales, ya que 


u · v = (2, −1, 3) · (1, 5, 1) = 2 − 5 + 3 = 0 



Los vectores de la base can´onica de R3 son ortogonales. Entonces se dice que la 
base can´onica de R3 es una base ortonormal, ya que sus vectores son ortogonales 
y unitarios. 


De forma an´aloga podemos deﬁnir el ´angulo entre dos vectores de R2. 


2.3.4. 
Producto vectorial 


Dados dos vectores u = (u1, u2, u3),v = (v1, v2, v3) ∈ R3, llamamos producto 
vectorial de u y v al vector que viene deﬁnido por el desarrollo formal del siguiente 
determinante por la primera ﬁla: 


u × v = 



ı 
 
k 
u1 
u2 
u3 
v1 
v2 
v3 


= (u2v3 − u3v2)ı − (u3v1 − u1v3) + (u1v2 − u2v1)k 


= (u2v3 − u3v2, −u3v1 + u1v3, u1v2 − u2v1) 


Ejemplo 2.3.6 
El producto vectorial de los vectores u = (0, −3, 5) y v = (3, 1, 10) es 


u × v = 



ı 
 
k 
0 
−3 
5 
3 
1 
10 


= −35ı + 5 + 3k = (−35, 5, 3) 



El producto vectorial cumple las siguientes propiedades: 


u × v = 0 si, y s´olo si, u y v son linealmente dependientes. 
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u × v = −v × u, para todo u,v ∈ R3. 


u × v es ortogonal a los vectores u y v, para todo u,v ∈ R3. 


u × (v + w) = u × v + u × w, para todo u,v, w ∈ R3. 


(u + v) × w = u × w + v × w, para todo u,v, w ∈ R3. 


(λu) × v = u × (λv) = λ(u × v), para todo λ ∈ R y u,v ∈ R3. 


Ejercicio 2.3.1 
Sean u,v y w tres vectores en R3 no coplanarios. Se˜nala en cada apartado si los vec- 
tores que se dan tienen la misma direcci´on, son ortogonales o si no puede asegurarse 
nada: 


(a) u × v, u 


(b) u × v, u × w 


(c) u · v, u 


(d) u − v, u × v 


(e) u − v, u + v 


(f ) (u + 2v) × u, u × v 



Adem´as: 
||u × v|| = ||u|| ||v|| sen α 


es el ´area del paralelogramo determinado por los vectores u y v, donde α ∈ [0, π] 
es el ´angulo que forman los vectores u,v ∈ R3. 


α 


u 


v 


Ejemplo 2.3.7 
El ´area del paralelogramo determinado por el vector u = (2, 3, 1) y por el vector 
v = (−1, −2, 0) es 
||u × v|| = 
√ 


6 


ya que 


u × v = 



ı 
 
k 
2 
3 
1 
−1 
−2 
0 


= (2, −1, −1) 


y as´ı 
||u × v|| = 
√ 


4 + 1 + 1 = 
√ 


6 
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2.3.5. 
Producto mixto 


Dados tres vectores u = (u1, u2, u3),v = (v1, v2, v3), w = (w1, w2, w3) ∈ R3, 
llamamos producto mixto de u, v y w al n´umero real 


[u,v, w] = u · (v × w) = 



u1 
u2 
u3 
v1 
v2 
v3 
w1 
w2 
w3 



Las propiedades del producto mixto se deducen de las propiedades de los productos 
escalar y vectorial, y de las propiedades de los determinantes. 
Adem´as: 
|[u,v, w]| 


es el volumen del paralelep´ıpedo cuyas aristas vienen determinadas por los vectores 
u, v y w. 


u 


v 


w 


Ejemplo 2.3.8 
El volumen del paralelep´ıpedo cuyas aristas vienen determinadas por los tres vec- 
tores u = (2, −1, 1), v = (−1, −2, 0) y w = (1, −1, 2) es 


|[u,v, w]| = 7 


ya que 


[u,v, w] = 



2 
−1 
1 
−1 
−2 
0 
1 
−1 
2 


= −7 



2.4. 
Geometr´ıa af´ın y eucl´ıdea 


En este apartado nos centramos en el espacio af´ın tridimensional, considerando 
las variedades lineales aﬁnes (rectas y planos) y estudiando las posiciones relativas 
entre ellas. Mediante la estructura eucl´ıdea vista en la secci´on anterior, introducimos 
el ´angulo y la distancia entre dos variedades lineales aﬁnes. 
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2.4.1. 
El espacio af´ın tridimensional 


El espacio af´ın tridimensional est´a formado por un conjunto E3 (a cuyos ele- 
mentos se les llama puntos), el espacio vectorial R3 y una aplicaci´on 


+ : E3 × R3 
E3 


(P,v) 
P + v 


que cumple las siguientes propiedades: 


P + (u + v) = (P + u) + v, para todo P ∈ E3 y para todo (u,v) ∈ R3 × R3. 


Dados (P, Q) ∈ E3 × E3 existe un ´unico vector v ∈ R3 tal que P + v = Q. 


Se llama vector ﬁjo a un par (P, Q) de puntos de E3 y se representa por −→ 
PQ. 
Se dice que P es el origen y Q el extremo del vector ﬁjo −→ 
PQ. 
Por la segunda propiedad anterior se tiene que dado el vector ﬁjo −→ 
PQ, es decir 
el par de puntos (P, Q), existe un ´unico vector v ∈ R3 tal que Q = P + v. 
El vector v se llama vector libre y se conviene en decir que el vector ﬁjo −→ 
PQ es un 
representante del vector v. 
Cuando se escribe 
v = −→ 
PQ 


se est´a identiﬁcando el vector ﬁjo −→ 
PQ con el vector v ∈ R3 que representa. 


P 


Q 


v = −→ 
P Q 


Con esta identiﬁcaci´on la primera propiedad se escribe 


−→ 
PQ + −→ 
QR = −→ 
PR 


para cada P, Q, R ∈ E3. 


P 


Q 


R 


u = −→ 
P Q 


u + v = −→ 
P Q + −→ 
QR 


v = −→ 
QR 
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Dados un punto O (origen) de E3 y una base {e1,e2,e3} de R3, se dice que 
{O;e1,e2,e3} es una referencia cartesiana de E3. 
Cuando la base es ortonormal, a la referencia se le llama referencia cartesiana 
rectangular. 
Se llaman coordenadas cartesianas (rectangulares) de un punto P ∈ E3 res- 
pecto a una referencia cartesiana (rectangular) {O;e1,e2,e3}, a las coordenadas 
(λ1, λ2, λ3) del vector −→ 
OP en la base {e1,e2,e3}. 


A partir de aqu´ı consideraremos la referencia cartesiana rectangular {O;ı,,k}, 
donde O ∈ E3 es un origen cualquiera. 
Dado P ∈ E3, si 
−→ 
OP = (p1, p2, p3) = p1ı + p2 + p3k 


escribiremos P = (p1, p2, p3), indicando que (p1, p2, p3) son las coordenadas del 
punto P en la referencia cartesiana rectangular {O;ı,,k}. 
La primera coordenada, p1, se llama abscisa. 
La segunda coordenada, p2, se llama ordenada. 
La tercera coordenada, p3, se llama cota. 
Si P = (p1, p2, p3) y v = (v1, v2, v3), entonces 


Q = P + v = (p1 + v1, p2 + v2, p3 + v3) 


Aqu´ı representamos gr´aﬁcamente el punto P = (2, 3, 4): 


X 
Y 


Z 


1.0 


2.0 


3.0 


1.0 
2.0 
3.0 


1.0 


2.0 


3.0 


4.0 


P 


Dados P, Q ∈ E3, se llama distancia entre los puntos P y Q al n´umero real 


d(P, Q) = ||−→ 
PQ|| = 


(q1 − p1)2 + (q2 − p2)2 + (q3 − p3)2 


siendo P = (p1, p2, p3) y Q = (q1, q2, q3); evidentemente, la distancia entre P y Q 
la podemos deﬁnir porque R3 es un espacio vectorial eucl´ıdeo. 


Ejemplo 2.4.1 
La distancia entre los puntos P = (3, −1, 2) y Q = (5, −4, 1) es 


d(P, Q) = 


(5 − 3)2 + (−4 − (−1))2 + (1 − 2)2 = 
√ 


22 + 32 + 1 = 
√ 


14 
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Entonces, se dice que el espacio af´ın tridimensional tiene estructura de espacio 
af´ın eucl´ıdeo. 


La distancia cumple las siguientes propiedades: 


d(P, Q) = 0 si, y s´olo si, P = Q. 


d(P, Q) = d(Q, P), para todo P, Q ∈ E3. 


d(P, R) d(P, Q) + d(Q, R), para todo P, Q, R ∈ E3. 


2.4.2. 
Variedades lineales aﬁnes 


Entenderemos como variedades lineales aﬁnes del espacio af´ın tridimensional 
las rectas y los planos. 
Recta 
Dado un punto P = (x0, y0, z0) y un vector no nulo v = (v1, v2, v3), se llama recta 
que pasa por P y tiene la direcci´on de v al conjunto formado por los puntos Q, tales 
que −→ 
PQ = λv, con λ ∈ R, es decir, tales que 


Q = P + λv , 
con λ ∈ R 


X 


Y 


Z 


P 


Q 


v 


El vector v se llama vector director de la recta. 


Si Q = (x, y, z), entonces 


(x, y, z) = (x0, y0, z0) + λ (v1, v2, v3) 


La ecuaci´on anterior se llama ecuaci´on vectorial de la recta. 


Las ecuaciones 
 
 


 


x = x0 + λ v1 
y = y0 + λ v2 
z = z0 + λ v3 
se llaman ecuaciones param´etricas de la recta. 


La expresi´on 
x − x0 


v1 
= y − y0 


v2 
= z − z0 


v3 
se llama ecuaci´on continua de la recta. 
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Ejemplo 2.4.2 
La recta que pasa por el punto P = (1, −2, 3) y tiene como vector director al vector 
v = (2, 4, 0), tiene como ecuaci´on vectorial 


(x, y, z) = (1, −2, 3) + λ(2, 4, 0) 


como ecuaciones param´etricas 
 
 


 


x = 
1 + 2λ 
y = −2 + 4λ 
z = 
3 


y como ecuaci´on continua 


x − 1 


2 
= y + 2 


4 
= z − 3 


0 


La expresi´on formal 
z − 3 


0 
en la ecuaci´on continua de la recta s´olo signiﬁca que los puntos de dicha recta tienen 
su tercera coordenada constante e igual a 3 (z = 3). 


Plano 
Dado un punto P = (x0, y0, z0) y dos vectores u = (u1, u2, u3) y v = (v1, v2, v3) 
linealmente independientes, se llama plano que pasa por el punto P y tiene a los 
vectores u y v como vectores directores, al conjunto formado por los puntos Q, tales 
que −→ 
PQ = λ u + µv, con λ, µ ∈ R, es decir, tales que 


Q = P + λ u + µv , 
con λ, µ ∈ R 


X 


Y 


Z 


P 
Q 


u 


v 


Si Q = (x, y, z), entonces 


(x, y, z) = (x0, y0, z0) + λ (u1, u2, u3) + µ (v1, v2, v3) 


Fundamentos Matem´aticos de la Ingenier´ıa – 503 – 1r. semestre 
2009/2010 
44 
c 
UJI 


52 
M. Barreda Rochera / J. A. López Ortí - ISBN: 978-84-692-9833-6 
Fundamentos Matemáticos de la Ingeniería. Parte I: Álgebra Lineal - UJI 


La ecuaci´on anterior se llama ecuaci´on vectorial del plano. 


Las ecuaciones 
 
 


 


x = x0 + λ u1 + µ v1 
y = y0 + λ u2 + µ v2 
z = z0 + λ u3 + µ v3 
se llaman ecuaciones param´etricas del plano. 


Ejemplo 2.4.3 
El plano que pasa por el punto P = (−1, 3, 1) y tiene como vectores directores a 
los vectores u = (2, 5, −1) y v = (−3, 2, −1) tiene como ecuaci´on vectorial 


(x, y, z) = (−1, 3, 1) + λ (2, 5, −1) + µ (−3, 2, −1) 


y como ecuaciones param´etricas 
 
 


 


x = −1 + 2λ − 3µ 
y = 
3 + 5λ + 2µ 
z = 
1 − λ − µ 



Como los vectores u = (u1, u2, u3) y v = (v1, v2, v3) son linealmente indepen- 
dientes, el sistema 
 
 


 


λ u1 + µ v1 = x − x0 
λ u2 + µ v2 = y − y0 
λ u3 + µ v3 = z − z0 
es compatible determinado (las inc´ognitas son λ y µ), es decir 


rang 


 


 
u1 
v1 
u2 
v2 
u3 
v3 


 


 = rang 


 


 
u1 
v1 
x − x0 
u2 
v2 
y − y0 
u3 
v3 
z − z0 


 


 = 2 


y por tanto 


u1 
v1 
x − x0 
u2 
v2 
y − y0 
u3 
v3 
z − z0 


= 0 


Desarrollando el determinante obtenemos la llamada ecuaci´on general del plano: 


A x + B y + C z + D = 0 


Ejemplo 2.4.4 
La ecuaci´on general del plano del Ejemplo 2.4.3 es 


2 
−3 
x + 1 
5 
2 
y − 3 
−1 
−1 
z − 1 


= 0 


es decir 
3x − 5y − 19z + 37 = 0 
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Dado el punto P = (x0, y0, z0) y los vectores u = (u1, u2, u3) y v = (v1, v2, v3) 
linealmente independientes, consideramos el vector 


n = u × v 


el cual es ortogonal a los vectores u y v, y se llama vector normal al plano. 


Ejemplo 2.4.5 
El vector normal del plano del Ejemplo 2.4.3 es 


n = u × v = (−3, 5, 19) 



Si Q = (x, y, z) es un punto del plano determinado por el punto P y los vectores 
u y v 


P 
Q 


n 


v 


u 


entonces los vectores −→ 
PQ y n son ortogonales. 
Si n = (A, B, C), entonces 


0 = −→ 
PQ·n = (x−x0, y−y0, z−z0)·(A, B, C) = A(x−x0)+B(y−y0)+C(z−z0) 


de donde 
A x + B y + C z + D = 0 


siendo D = −A x0 − B y0 − C z0. 


2.4.3. 
Posiciones relativas entre variedades lineales aﬁnes 


Posiciones relativas entre dos rectas 
Sean las rectas 
r ≡ P + λ u , 
s ≡ Q + µv 


donde P = (a, b, c), u = (u1, u2, u3), Q = (a′, b′, c′) y v = (v1, v2, v3). 


1 Si los vectores u y v son linealmente dependientes 


rang 
u1 
u2 
u3 
v1 
v2 
v3 


= 1 


entonces las rectas r y s son paralelas, en cuyo caso lo designamos por r s, 
o son coincidentes, en cuyo caso escribimos r = s. 
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a Si el punto P pertenece a la recta s o, equivalentemente, el punto Q 
pertenece a la recta r, entonces r = s. 


P 
Q 
u 
v 


Ejemplo 2.4.6 
Sean las rectas 


r ≡ x + 1 


2 
= 5 − y = z + 3 
y 
s ≡ 


 
 


 


x = 7 − 2λ 
y = 1 + 
λ 
z = 1 − 
λ 


La recta r pasa por el punto P = (−1, 5, −3) y tiene como vector director al vector 
u = (2, −1, 1), mientras que la recta s pasa por el punto Q = (7, 1, 1) y tiene como 
vector director al vector v = (−2, 1, −1). Por una parte 


v = (−2, 1, −1) = −u 


es decir, los vectores u y v son linealmente dependientes. Por otra parte 


7 + 1 


2 
= 5 − 1 = 1 + 3 


esto es, Q ∈ s. Por tanto, las rectas r y s son coincidentes. 


b Si el punto P no pertenece a la recta s o, equivalentemente, el punto Q 
no pertenece a la recta r, entonces r s. 


P 


Q 


u 


v 


r 


s 


Ejemplo 2.4.7 
Sean las rectas 


r ≡ 


 
 


 


x = −1 + 2λ 
y = 
5 − 
λ 
z = −3 + 
λ 
y 
s ≡ 5 − x 


2 
= y − 1 = 1 − z 


La recta r pasa por el punto P = (−1, 5, −3) y tiene como vector director al vector 
u = (2, −1, 1); la recta s pasa por el punto Q = (5, 1, 1) y tiene como vector 
director al vector v = (−2, 1, −1). Por una parte 


u = (2, −1, 1) = −v 
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es decir, los vectores u y v son linealmente dependientes. Por otra parte 


5 − (−1) 


2 
= 3 = 4 = 5 − 1 = 1 − (−3) 


esto es, P ∈ s. Por tanto, las rectas r y s son paralelas. 


2 Si los vectores u y v son linealmente independientes y los vectores u, v y −→ 
PQ 
son linealmente dependientes 


rang 
u1 
u2 
u3 
v1 
v2 
v3 


= 2 , 
u,v, −→ 
PQ 
= 



u1 
u2 
u3 
v1 
v2 
v3 
a′ − a 
b′ − b 
c′ − c 


= 0 


entonces las rectas r y s se cortan. 


Q 


P 
u 


v 


r 


s 


Ejemplo 2.4.8 
Sean las rectas 


r ≡ x − 1 


2 
= −y − 2 = z − 1 


0 
y 
s ≡ 


 
 


 


x = 
1 + 
λ 
y = 
1 − 2λ 
z = −5 + 3λ 


La recta r pasa por el punto P = (1, −2, 1) y tiene como vector director al vector 
u = (2, −1, 0); la recta s pasa por el punto Q = (1, 1, −5) y tiene como vector 
director al vector v = (1, −2, 3). Por una parte, los vectores u = (2, −1, 0) y 
v = (1, −2, 3) son linealmente independientes, ya que 


rang 
2 
−1 
0 
1 
−2 
3 


= 2 


al tener esta matriz un menor de orden 2 distinto de cero: 
2 
−1 
1 
−2 


= −3 = 0 


Por otra parte, 


u,v, −→ 
PQ 
= 



2 
−1 
0 
1 
−2 
3 
1 − 1 
1 − (−2) 
−5 − 1 


= 



2 
−1 
0 
1 
−2 
3 
0 
3 
−6 


= 0 
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Por tanto, las rectas r y s se cortan (en un punto). 
El punto de corte se obtiene, por ejemplo, resolviendo el sistema 


1 + λ − 1 


2 
= −(1 − 2λ) − 2 = −5 + 3λ − 1 


0 
es decir 
λ 
2 = 2λ − 3 = −6 + 3λ 


0 
Por tanto, λ = 2. Entonces, el punto de corte es (3, −3, 1). 
Aqu´ı hay que hacer notar que la ecuaci´on 


2λ − 3 = −6 + 3λ 


0 
signiﬁca que −6 + 3λ = 0. 


3 Si los vectores u, v y −→ 
PQ son linealmente independientes 


u,v, −→ 
PQ 
= 



u1 
u2 
u3 
v1 
v2 
v3 
a′ − a 
b′ − b 
c′ − c 


= 0 


entonces las rectas r y s se cruzan. 


P 
u 
r 


Q 
v 


s 


Ejemplo 2.4.9 
Sean las rectas 


r ≡ x − 1 


2 
= −y − 2 = z − 1 


0 
y 
s ≡ 


 
 


 


x = −3 + 
λ 
y = 
2 − 2λ 
z = −5 + 3λ 


La recta r pasa por el punto P = (1, −2, 1) y tiene como vector director al vector 
u = (2, −1, 0); la recta s pasa por el punto Q = (−3, 2, −5) y tiene como vector 
director al vector v = (1, −2, 3). Entonces, el vector −→ 
PQ = (−4, 4, −6). 
Adem´as 


u,v, −→ 
PQ 
= 



2 
−1 
0 
1 
−2 
3 
−4 
4 
−6 


= 



2 
−1 
0 
1 
−2 
3 
−4 
4 
−6 


= 6 = 0 


Por tanto, las rectas r y s se cruzan. 
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Posiciones relativas entre dos planos 
Sean los planos 


π ≡ A x + B y + C z + D = 0 , 
π′ ≡ A′x + B′y + C′z + D′ = 0 


y sean P ∈ π, P ′ ∈ π′, y n = (A, B, C), n′ = (A′, B′, C′) los vectores normales a 
los planos π y π′, respectivamente. 


1 Si los vectores n y n′ son linealmente dependientes 


rang 
A 
B 
C 
A′ 
B′ 
C′ 


= 1 


entonces los planos π y π′ son paralelos, en cuyo caso lo designamos por 
π π′, o son coincidentes, en cuyo caso escribimos π = π′. 


a Si el punto P pertenece al plano π′ o, equivalentemente, el punto P ′ 


pertenece al plano π, entonces π = π′. 


P 
P ′ 


n 


n′ 


π = π′ 


Ejemplo 2.4.10 
Veamos que los planos 


π ≡ 2x − y + z − 1 = 0 , 
π′ ≡ 


 
 


 


x = 1 + 2λ − 3µ 
y = 4 − 
λ + 2µ 
z = 3 − 5λ + 8µ 


son coincidentes. 
Por una parte, el vector n = (2, −1, 1) es el vector normal al plano π y el punto 
P = (0, 0, 1) ∈ π. 
Por otra parte, como 


ı 
 
k 
2 
−1 
−5 
−3 
2 
8 


= 2ı −  + k 


se tiene que n′ = (2, −1, 1) es el vector normal al plano π′ y, adem´as, el punto 
Q = (1, 4, 3) ∈ π′. 
Entonces, los vectores n y n′ son linealmente independientes (son iguales) y Q ∈ π. 
Por tanto, los planos π y π′ son coincidentes. 
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b Si el punto P no pertenece al plano π′ o, equivalentemente, el punto P ′ 


no pertenece al plano π, entonces π π′. 


P 


n 


π 


P ′ 
π′ 


n′ 


Ejemplo 2.4.11 
Veamos que los planos 


π ≡ 2x − 4y + 5z − 2 = 0 , 
π′ ≡ 4x − 8y + 10z + 8 = 0 


son paralelos. 


Los vectores n = (2, −4, 5) y n′ = (4, −8, 10) son los vectores normales a los 
planos π y π′, respectivamente; estos vectores son linealmente dependientes, ya que 
n′ = 2n, y el punto P = (1, 0, 0) ∈ π, pero no pertenece al plano π′ (P /∈ π′). 
Por tanto, los planos π y π′ son paralelos. 


2 Si los vectores n y n′ son linealmente independientes 


rang 
A 
B 
C 
A′ 
B′ 
C′ 


= 2 


entonces los planos π y π′ se cortan (en una recta r). 


P 
P ′ 


n 
n′ 


π 
π′ 


r 


La recta r queda determinada por las ecuaciones: 
A x + B y + C z + D = 0 
A′x + B′y + C′z + D′ = 0 


Estas ecuaciones se llaman ecuaciones impl´ıcitas de la recta r. 
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Ejemplo 2.4.12 
Veamos que los planos 


π ≡ 2x − 3y + 4z + 1 = 0 , 
π′ ≡ x + y − z + 2 = 0 


se cortan (en una recta r). 
Los vectores n = (2, −3, 4) y n′ = (1, 1, −1) son los vectores normales a los planos 
π y π′, respectivamente. Estos vectores son linealmente independientes, ya que 


rang 
2 
−3 
4 
1 
1 
−1 


= rang 
2 
−3 
4 
0 
−1 
6 


= 2 


Por tanto, los planos π y π′ se cortan. 
La recta r tiene como ecuaciones impl´ıcitas 


r ≡ 
2x − 3y + 4z + 1 = 0 
x + 
y − 
z + 2 = 0 



Posiciones relativas entre una recta y un plano 
Sea la recta 
r ≡ P + λv 


donde P = (a, b, c) y v = (v1, v2, v3), y sea el plano 


π ≡ A x + B y + C z + D = 0 


siendo n = (A, B, C) el vector normal al plano π. 


1 Si los vectores v y n son ortogonales, entonces la recta r y el plano π son 
paralelos, en cuyo caso lo designamos por r π, o la recta r est´a contenida 
en el plano π, en cuyo caso escribimos r ⊂ π. 


a Si P ∈ π, entonces r ⊂ π. 


P 


n 


π 


r 
v 
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Ejemplo 2.4.13 
Veamos que la recta 


r ≡ x − 3 = y + 2 


2 
= −z + 4 


est´a contenida en el plano π ≡ x − y − z − 1 = 0. 
La recta r pasa por el punto P = (3, −2, 4) y tiene como vector director al vector 
v = (1, 2, −1). El plano π tiene vector normal n = (1, −1, −1). 
Entonces 
v · n = (1, 2, −1) · (1, −1, −1) = 1 − 2 + 1 = 0 


Adem´as, P ∈ π, ya que 3 − (−2) − 4 − 1 = 0. 
Por tanto, la recta r est´a contenida en el plano π. 


b Si P /∈ π, entonces r π. 


P 


n 


π 


r 


v 


Q 


Ejemplo 2.4.14 
Veamos que la recta 


r ≡ x − 3 = y + 2 


2 
= −z + 4 


es paralela al plano π ≡ x − y − z + 2 = 0. 
La recta r pasa por el punto P = (3, −2, 4) y tiene como vector director al vector 
v = (1, 2, −1). El plano π tiene vector normal n = (1, −1, −1). 
Entonces 
v · n = (1, 2, −1) · (1, −1, −1) = 1 − 2 + 1 = 0 


Adem´as, P /∈ π, ya que 3 − (−2) − 4 + 2 = 3 = 0. 
Por tanto, la recta r es paralela al plano π. 


2 Si los vectores v y n no son ortogonales (v ⊥ n), entonces la recta r y el 
plano π se cortan (en un punto). 
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P 


n 


π 


r 


v 


Q 


Ejemplo 2.4.15 
Veamos que la recta 


r ≡ 


 
 


 


x = 
1 + 2λ 
y = −1 − 
λ 
z = 
2 + 3λ 


y el plano π ≡ 2x + y + 4z + 10 = 0 se cortan (en un punto). 
La recta r pasa por el punto P = (1, −1, 2) y tiene como vector director al vector 
v = (2, −1, 3). El plano π tiene vector normal n = (2, 1, 4). 
Entonces 
v · n = (2, −1, 3) · (2, 1, 4) = 4 − 1 + 12 = 15 = 0 


Por tanto, la recta r y el plano π se cortan (en un punto). 
Para obtener el punto de corte resolvemos la ecuaci´on: 


2(1 + 2λ) + (−1 − λ) + 4(2 + 3λ) + 10 = 0 


es decir 
19 + 15λ = 0 


Por tanto 
λ = −19 


15 
y as´ı 


x = 1 − 2 19 


15 = −23 


15 , 
y = −1 + 19 


15 = − 4 


15 , 
z = 2 − 3 19 


15 = −27 


15 = −9 


5 


es decir, el punto de corte es Q = (−23/15, −4/15, −9/5). 
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Ejercicio 2.4.1 
Estudia la posici´on relativa de los siguientes pares de variedades lineales aﬁnes: 


(a) El plano 3x + 2z − 1 = 0 y la recta 
 
 


 


x = 
1 
+ 4λ 
y = −λ 
z = 
2 


(b) Las rectas 
x − 1 


2 
= y − 1 


2 
= z − 3 


3 
, 
x − z = 0 
y = 3 


(c) Los planos 2x − 3y + z − 1 = 0 y x + y − 2z + 1 = 0. 



Ejercicio 2.4.2 
Estudia, en funci´on del par´ametro real a, la posici´on relativa de los siguientes pares 
de variedades lineales aﬁnes: 


(a) 


 
 


 


x 
= 
1 
+ 
α 
− 
β 
y 
= 
2 
− 
α 
+ 
β 
z 
= 
1 
+ 
α 
, 


 
 


 


x 
= 
1 
− 
a α 
+ 
β 
y 
= 
α 
− 
β 
z 
= 
3 
+ 
a α 
+ 
2β 


(b) 
2x 
+ 
3y 
= 
7 
a x 
− 
3z 
= 
−5 
, 
6x − 4y − 3z + 3 = 0 


(c) −2x + y + 4 = 0 , 


 
 


 


x 
= 
3 
− 
α 
y 
= 
2 
− 
2α 
z 
= 
1 
+ 
a β 


(d) 3x − y − z = 0 , 
−6x + a y + 2z + 1 = 0. 


(e) 2x − y + 1 = 0 , 
x + 4y + a z − 3 = 0. 



2.4.4. 
´Angulo entre variedades lineales aﬁnes 


´Angulo entre dos rectas que se cortan 
Sean dos rectas que se cortan (en un punto P): 


r ≡ P + λ u , 
s ≡ P + µv 


P 


u 


v 


r 


s 


α 
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Entonces 


cos α = 
|u · v| 


||u|| ||v|| 


Ejemplo 2.4.16 
Sean las rectas r y s del Ejemplo 2.4.8, las cuales, como ya vimos, se cortan. 
Los vectores directores de las rectas r y s son, respectivamente, u = (2, −1, 0) y 
v = (1, −2, 3). Entonces, el coseno del ´angulo que forman las rectas r y s es 


cos α = 
|u · v| 


||u|| ||v|| = 
|(2, −1, 0) · (1, −2, 3)| 


||(2, −1, 0)|| ||(1, −2, 3)|| = 
4 
√ 


5 
√ 


14 


Por tanto 
α = arc cos(4/ 
√ 


70) ≈ 1.07 radianes 



´Angulo entre dos planos que se cortan 
Sean dos planos que se cortan (en una recta r): 


π ≡ A x + B y + C z + D = 0 , 
π′ ≡ A′x + B′y + C′z + D′ = 0 


siendo n = (A, B, C) y n′ = (A′, B′, C′) los vectores normales a los planos π y π′, 
respectivamente. 


n 
n′ 


π 
π′ 


r 


α 


Entonces 


cos α = 
|n · n′| 


||n|| ||n′|| 


Ejemplo 2.4.17 
Sean los planos π y π′ del Ejemplo 2.4.12, los cuales, como ya vimos, se cortan. 
Los vectores normales de los planos π y π′ son, respectivamente, n = (2, −3, 4) y 
n′ = (1, 1, −1). Entonces, el coseno del ´angulo que forman los planos π y π′ es 


cos α = 
|n · n| 


||n|| ||n|| = 
|(2, −3, 4) · (1, 1, −1)| 


||(2, −3, 4)|| ||(1, 1, −1)|| = 
5 
√ 


29 
√ 


3 


Por tanto 
α = arc cos(5/ 
√ 


87) ≈ 1.005 radianes 
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´Angulo entre un plano y una recta que se cortan 
Sean un plano y una recta que se cortan (en un punto P): 


π ≡ A x + B y + C z + D = 0 , 
r ≡ P + λv 


donde n = (A, B, C) es el vector normal al plano π. 


n 


v 


r 


π 
P 


β 


α 


Entonces 


sen α = 
|v · n| 


||v|| ||n|| , 
β = π 


2 − α 


Ejemplo 2.4.18 
Sean la recta r y el plano π del Ejemplo 2.4.15, los cuales, como ya vimos, se cortan. 
El vector director de la recta r es v = (2, −1, 3) y el vector normal al plano π es 
n = (2, 1, 4). Entonces, el seno del ´angulo que forman la recta r y el plano π es 


sen α = 
|v · n| 


||v|| ||n|| = 
|(2, −1, 3) · (2, 1, 4)| 


||(2, −1, 3)|| ||(2, 1, 4)|| = 
15 
√ 


14 
√ 


21 


Por tanto 
α = arc cos(14/ 
√ 


294) ≈ 1.065 radianes 



2.4.5. 
Distancia entre variedades lineales aﬁnes 


Distancia de un punto a una recta 
Sea una recta r ≡ P + λv y sea un punto Q /∈ r. 


P 


Q 


r 
α 


−→ 
P Q 


d 


v 
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Entonces, la distancia entre el punto Q y la recta r, que denotamos por d(Q, r), es: 


d(Q, r) = ||v × −→ 
PQ|| 


||v|| 


Ejemplo 2.4.19 
Sean el punto Q = (1, −2, 3) y la recta determinada por el punto P = (−1, 3, 4) y 
el vector (director) v = (2, −1, 7) 


r ≡ x + 1 


2 
= −y + 3 = z − 4 


7 


Entonces, como −→ 
PQ = (2, −5, −1), obtenemos 


v × −→ 
PQ = 



ı 
 
k 
2 
−1 
7 
2 
−5 
−1 


= 36ı + 16 − 8k 


Por tanto 


d(Q, r) = ||v × −→ 
PQ|| 


||v|| 
= 4 
√ 


101 


3 
√ 


6 
≈ 5.47 



Distancia de un punto a un plano 
Sea un plano π ≡ A x + B y + C z + D = 0 y sea un punto Q = (x0, y0, z0) /∈ π. 


n 
−→ 
P Q 


Q 


π 
P 


d 
α 


Entonces, la distancia entre el punto Q y el plano π, que denotamos por d(Q, π), es: 


d(Q, π) = |A x0 + B y0 + C z0 + D| 
√ 


A2 + B2 + C2 


Ejemplo 2.4.20 
Sean el punto Q = (1, −1, 2) y el plano π ≡ 2x − 5y + 3z + 5 = 0. Entonces, la 
distancia entre el punto Q y el plano π es 


d(Q, π) = 18 
√ 


38 ≈ 2.92 
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Distancia entre dos rectas paralelas 
Sean dos rectas paralelas 


r ≡ P + λv , 
r ≡ Q + µv , 


Q 


P 


s 


r 


d 
d 


v 


v 


Entonces, la distancia entre las rectas r y s, que denotamos por d(r, s), es: 


d(r, s) = d(P, s) = d(Q, r) = ||v × −→ 
PQ|| 


||v|| 


Ejemplo 2.4.21 
Sean las rectas r y s del Ejemplo 2.4.7, las cuales, como ya vimos, son paralelas. 
La recta r est´a determinada por el punto P = (−1, 5, −3) y el vector (director) 
u = (2, −1, 1); la recta s est´a determinada por el punto Q = (5, 1, 1) y el vector 
(director) u = (2, −1, 1). Entonces, el vector −→ 
PQ = (6, −4, 4) y as´ı 


v × −→ 
PQ = 



ı 
 
k 
2 
−1 
1 
6 
−4 
4 


= −2 − 2k 


Por tanto 


d(r, s) = d(P, s) = ||v × −→ 
PQ|| 


||v|| 
= 2 
√ 


2 
√ 


6 
= 2 
√ 


3 
≈ 1.1547 



Distancia entre dos rectas que se cruzan 
Sean dos rectas que se cruzan 


r ≡ P + λ u , 
s ≡ Q + µv 


es decir, los vectores u, v y −→ 
PQ son linealmente independientes. 


P 
u 
r 


Q 
v 


s 
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Entonces 


d(r, s) = 


u,v, −→ 
PQ 


||u × v|| 


Ejemplo 2.4.22 
Sean las rectas r y s del Ejemplo 2.4.9, las cuales, como ya vimos, se cruzan. 
La recta r pasa por el punto P = (1, −2, 1) y tiene como vector director al vector 
u = (2, −1, 0); la recta s est´a determinada por el punto Q = (−3, 2, −5) y el vector 
(director) u = (1, −2, 3). Adem´as, tal y como calculamos en el Ejemplo 2.4.9: 
u,v, −→ 
PQ 
= 6 


Por otra parte 


u × v = 



ı 
 
k 
2 
−1 
0 
1 
−2 
3 


= −3ı − 6 − 3k 


Por tanto 


d(r, s) = 


u,v, −→ 
PQ 


||u × v|| 
= 
6 
√ 


54 = 2 
√ 


6 ≈ 0.816 



Distancia entre dos planos paralelos 
Sean dos planos paralelos 


π ≡ A x + B y + C z + D = 0 , 
π′ ≡ A x + B y + C z + D′ = 0 


π 


π′ 
d 
n 


n 


Entonces 


d(π, π′) = 
|D′ − D| 
√ 


A2 + B2 + C2 
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Ejemplo 2.4.23 
Sean los planos del Ejemplo 2.4.11, los cuales, como ya indicamos, son paralelos. 
Entonces 


d(π, π′) = 
|4 − (−2)| 


22 + (−4)2 + 32 = 2 
√ 


5 ≈ 0.89 



Distancia entre una recta y un plano paralelos 
Sean una recta y un plano paralelos 


r ≡ P + λv , 
π ≡ A x + B y + C z + D = 0 


es decir, los vectores n = (A, B, C) y v son ortogonales. 


n 


P 
r 
v 


π 


d 


Si P = (x0, y0, z0), entonces 


d(r, π) = d(P, π) = |A x0 + B y0 + C z0 + D| 
√ 


A2 + B2 + C2 


Ejemplo 2.4.24 
Sean la recta r y el plano π del Ejemplo 2.4.14. 
La recta r pasa por el punto P = (3, −2, 4) y el plano π tiene como ecuaci´on 
general: x − y − z + 2 = 0. Entonces 


d(r, π) = d(P, π) = 
|3 + 2 − 4 + 2| 


12 + (−1)2 + (−1)2 = 3 
√ 


3 = 
√ 


3 ≈ 1.73 



Ejercicio 2.4.3 
Dado el plano 3x + 4z + a = 0 y la recta 


x − 1 


4 
= y 


−1 = z 


b 


halla el valor de a y b de forma que la distancia de la recta al plano sea de 2 unidades. 
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Ejercicio 2.4.4 
Determina la posici´on relativa entre las siguientes variedades lineales aﬁnes: 


2x − 3y + 
z = 11 
x + 
y + 3z = 
8 , 


 
 


 


x 
= 
1 
+ 
α 
+ 
β 
y 
= 
1 
− 
α 
+ 
8β 
z 
= 
5 
− 
α 
+ 
5β 


Si son paralelas, calcula la distancia entre ambas y si no lo son, calcula el ´angulo 
que forman. 


Ejercicio 2.4.5 
Obt´en la ecuaci´on de un plano que dista 4 unidades del origen de coordenadas, es 
perpendicular a la recta x = y = z y corta al eje OZ en un punto de cota negativa. 


Ejercicio 2.4.6 
Considera las rectas 


r ≡ 
x + y + z = 2 
x + 2y − 3z = 8 , 
s ≡ 
ax − y − z = 1 
x − y + z = −2 


¿Para qu´e valor del par´ametro a se cortan las rectas r y s? Para dicho valor, calcula 
el punto de corte de ambas rectas y el ´angulo que forman. 
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Tema 3 


Espacios vectoriales 


En este tema estudiamos el concepto de espacio vectorial, centr´andonos en los 
espacios vectoriales reales de dimensi´on ﬁnita, y m´as concretamente en el espacio 
vectorial Rn. No obstante, creemos necesario introducir otras estructuras algebrai- 
cas, con el ﬁn de situar los espacios vectoriales en el contexto adecuado. 


3.1. 
Introducci´on 


En este primer apartado vamos a formalizar algunos conceptos que ya hemos 
utilizado en los dos temas anteriores, como son las leyes de composici´on (operacio- 
nes) sobre un conjunto, las cuales nos permitir´an introducir la estructura de cuerpo, 
con la que tambi´en hemos trabajado en el caso particular del cuerpo de los n´umeros 
reales. 


3.1.1. 
Leyes de composici´on 


Ley de composici´on interna 
Sea A un conjunto no vac´ıo. Una ley de composici´on interna sobre A es una apli- 
caci´on de A × A en A; por tanto, a cada (a, b) ∈ A × A le hacemos corresponder 
un ´unico elemento c ∈ A: 
A × A −→ A 
(a, b) −→ c 


Ejemplo 3.1.1 
La suma de n´umeros naturales es una ley de composici´on interna sobre N, ya que 
a cada par de n´umeros naturales (m, n) le hacemos corresponder el n´umero natural 
m + n: 
+ : N × N −→ N 
(m, n) −→ m + n 



Ejemplo 3.1.2 
La resta de n´umeros naturales no es una ley de composici´on interna sobre N, ya 
que, por ejemplo: 
2 − 3 = −1 /∈ N 
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En cambio, la resta de n´umeros enteros s´ı que es una ley de composici´on interna 
sobre Z, ya que a cada (m, n) ∈ Z le hacemos corresponder el ´unico n´umero entero 
m − n: 
− : Z × Z −→ Z 
(m, n) −→ m − n 



Ley de composici´on externa 
Sean un conjunto no vac´ıo A y otro conjunto no vac´ıo K, al que llamaremos do- 
minio de operadores. Una ley de composici´on externa sobre A con dominio de 
operadores K es una aplicaci´on de K × A en A; por tanto, a cada (k, a) ∈ K × A 
le hacemos corresponder un ´unico elemento b ∈ A: 


K × A −→ A 
(k, a) −→ b 


Ejemplo 3.1.3 
El producto de un n´umero natural y un n´umero entero es una ley de composici´on 
externa sobre Z con dominio de operadores N, ya que a cada (m, n) ∈ N × Z le 
hacemos corresponder el n´umero entero m · n: 


· : N × Z −→ Z 
(m, n) −→ m · n 


En cambio, el producto de un n´umero entero y un n´umero natural no es una ley de 
composici´on externa sobre N con dominio de operadores Z, ya que, por ejemplo: 


−2 · 3 = −6 /∈ N 



Las leyes de composici´on, en ocasiones, tambi´en se llaman operaciones. 


3.1.2. 
Estructuras algebraicas 


Una estructura algebraica est´a formada por un conjunto en el que hay deﬁnidas 
ciertas operaciones cumpliendo unas determinadas propiedades. 


Un grupo es una estructura algebraica formada por un conjunto G y una ley de 
composici´on interna 
∗ : G × G −→ G 
(a, b) −→ a ∗ b 


cumpliendo las siguientes propiedades: 


Dados cualesquiera a, b, c ∈ G, se tiene que a ∗ (b ∗ c) = (a ∗ b) ∗ c. 


(Propiedad asociativa) 


Existe un e ∈ G tal que, dado cualquier a ∈ G se tiene que a ∗ e = a = e ∗ a. 


(Existencia de elemento neutro) 


Fundamentos Matem´aticos de la Ingenier´ıa – 503 – 1r. semestre 
2009/2010 
64 
c 
UJI 


72 
M. Barreda Rochera / J. A. López Ortí - ISBN: 978-84-692-9833-6 
Fundamentos Matemáticos de la Ingeniería. Parte I: Álgebra Lineal - UJI 


Dado cualquier a ∈ G, existe un a′ ∈ G tal que a ∗ a′ = e = a′ ∗ a. 


(Existencia de elemento sim´etrico) 


Un grupo es, por tanto, un par (G; ∗), donde G es un conjunto y ∗ es una ley de 
composici´on interna sobre G cumpliendo las tres propiedades anteriores. 


Ejercicio 3.1.1 
Sea (G; ∗) un grupo. Entonces: 


El elemento neutro e ∈ G es ´unico. 


Dado cualquier a ∈ G, se tiene que el elemento sim´etrico a′ ∈ G es ´unico. 



Un grupo (G; ∗) se dice que es un grupo abeliano o grupo conmutativo si para 
cualesquiera a, b ∈ G se tiene que 


a ∗ b = b ∗ a 
(Propiedad conmutativa) 


Ejemplo 3.1.4 
El conjunto de los n´umeros enteros, Z, con la suma tiene estructura de grupo abe- 
liano, ya que la suma es asociativa y conmutativa, el elemento neutro es el n´umero 
0 y el elemento sim´etrico de m ∈ Z es −m ∈ Z. Se denota por (Z; +). 
Lo mismo sucede con el conjunto de los n´umeros racionales, Q, y con el conjunto 
de los n´umeros reales, R. Estos grupos se representan por (Q; +) y (R; +), respec- 
tivamente. 
En cambio, el conjunto de los n´umeros enteros, Z, con el producto no tiene estruc- 
tura de grupo, ya que los ´unicos n´umeros enteros que tienen sim´etrico son el 1 y el 
−1. Evidentemente, el elemento neutro es el n´umero 1. 


Un anillo es una estructura algebraica formada por un conjunto A y dos leyes de 
composici´on interna sobre A: suma (+) y producto (·), cumpliendo las siguientes 
propiedades: 


(A; +) es un grupo abeliano. 


Dados cualesquiera a, b, c ∈ A, se tiene que 


a · (b · c) = (a · b) · c 
(Propiedad asociativa) 


Dados cualesquiera a, b, c ∈ A, se tiene que 


a · (b + c) = a · b + a · c 
(b + c) · a = b · a + c · a 


(Propiedades distributivas a izquierda y derecha, respectivamente, del produc- 
to respecto de la suma) 
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Un anillo es, por tanto, una terna (A; +, ·), donde (A; +) es un grupo abeliano, 
y el producto es una operaci´on asociativa y distributiva respecto de la suma. 
El elemento neutro de la suma se representa por el s´ımbolo 0 (no confundir con 
el n´umero 0) y el elemento sim´etrico de a ∈ A (respecto de la suma) se representa 
por −a y se llama elemento opuesto de a. 


Un anillo se dice que es un anillo unitario si el producto tiene elemento neutro. 
En este caso, el elemento neutro del producto tambi´en es ´unico; se denota por el 
s´ımbolo 1 (no confundir con el n´umero 1). 
Adem´as, se dice que un anillo es un anillo conmutativo si el producto es con- 
mutativo. 


Ejercicio 3.1.2 
Sea (A; +, ·) un anillo. Entonces: 


Para cualquier a ∈ A, se tiene que 


a · 0 = 0 = 0 · a 


Para cualesquiera a, b ∈ A, se tiene que 


(−a) · b = a · (−b) = −(a · b) 


Para cualesquiera a, b ∈ A, se tiene que 


(−a) · (−b) = a · b 



Ejemplo 3.1.5 
Los conjuntos num´ericos Z, Q y R son, con las operaciones suma y producto usua- 
les, anillos unitarios y conmutativos, ya que como hemos visto en el Ejemplo 3.1.4 
estos conjuntos con la suma tienen estructura de grupo y, adem´as, el producto tiene 
elemento neutro (el n´umero 1), es asociativo, es conmutativo y se cumple la propie- 
dad distributiva del producto respecto de la suma. 
Se representan por (Z; +, ·), (Q; +, ·) y (R; +, ·), respectivamente. 


Un anillo (K; +, ·) se dice que es un cuerpo si 


(K; +, ·) es un anillo unitario y conmutativo. 


Para cualquier a ∈ K, a = 0, existe a−1 ∈ K tal que 


a · a−1 = 1 = a−1 · a 


El elemento a−1 es el elemento sim´etrico de a respecto del producto y se 
llama inverso del elemento a. Al igual que suced´ıa con la suma, el elemento 
inverso a−1 de un elemento a ∈ K es ´unico. 
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Los elementos de un cuerpo se suelen llamar escalares. 


Ejemplo 3.1.6 
Los conjuntos num´ericos Q y R son, con las operaciones suma y producto, cuerpos 
conmutativos, ya que como hemos visto en el Ejemplo 3.1.5 estos conjuntos tienen 
estructura de anillo unitario y conmutativo y, adem´as, todo n´umero racional (real) 
no nulo tiene inverso. 
En cambio el anillo (Z; +, ·) no es un cuerpo, ya que como hemos visto en el Ejem- 
plo 3.1.4 los ´unicos n´umero enteros que tienen inverso son el 1 y el −1. 


Ejemplo 3.1.7 
En el conjunto 
R2 = R × R = {(x, y) : x ∈ R , y ∈ R} 


se consideran las operaciones 


Suma : 
(a, b) + (c, d) = (a + c, b + d) 
Producto : 
(a, b) · (c, d) = (ac − bd, ad + bc) 


Veamos que con estas dos operaciones el conjunto R2 tiene estructura de cuerpo; se 
denota por C y se llama el cuerpo de los n´umeros complejos. Los elementos de C 
se llaman n´umeros complejos. 


Es f´acil comprobar que la suma y el producto de n´umeros complejos es asociativa 
y conmutativa. El elemento neutro de la suma es (0, 0) y el elemento neutro del 
producto es (1, 0). El elemento opuesto de (a, b) es (−a, −b) y elemento inverso de 
(a, b) = (0, 0) es 


(a, b)−1 = 
a 


a2 + b2, 
−b 


a2 + b2 



Si identiﬁcamos el conjunto R × {0} con R, podemos considerar el conjunto de 
los n´umeros reales como un subconjunto del conjunto de los n´umeros complejos. 
De esta forma identiﬁcamos el elemento (a, 0) con el n´umero real a. Adem´as, si 
denotamos el n´umero complejo (0, b) como b i, podemos representar todo n´umero 
complejo en forma bin´omica: 


(a, b) = a + b i 


El n´umero real a se llama parte real del n´umero complejo a + b i, y el n´umero real 
b se llama parte imaginaria del n´umero complejo a + b i. 
De esta forma i = (0, 1) y as´ı 


i2 = i · i = (0, 1) · (0, 1) = (−1, 0) = −1 


es decir, i es una ra´ız de la ecuaci´on z2 + 1 = 0; la otra ra´ız es −i. 
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El concepto de matriz que hemos visto en el Tema 1 (matriz real) se puede 
generalizar a cualquier cuerpo K (cuando no hay confusi´on se suelen omitir las 
operaciones a la hora de designar un cuerpo). 
As´ı, dado un cuerpo (K; +, ·), llamamos matriz de tama˜no m×n sobre el cuerpo 
K a una colecci´on de m · n escalares, que representamos en una tabla rectangular 
de m ﬁlas y n columnas 


A = 


 


 
 
 
 


a11 
a12 
· · · 
a1n 
a21 
a22 
· · · 
a2n 
... 
... 
... 
... 
am1 
am2 
· · · 
amn 


 


 
 
 
 


donde aij ∈ K, con i = 1, 2, . . . , m, y j = 1, 2, . . . , n. 
El conjunto de todas las matrices de tama˜no m × n sobre el cuerpo K lo repre- 
sentamos por Mm×n(K). 
Todos los conceptos introducidos en el Tema 1, los cuales se deﬁnieron para 
el cuerpo de los n´umeros reales, se pueden generalizar cuando K es un cuerpo 
arbitrario. 
Evidentemente, algunos de esos conceptos requieren de unas deﬁniciones m´as 
precisas como, por ejemplo, la noci´on de determinante de una matriz. 


Ejemplo 3.1.8 
El conjunto Mm×n(K) con la operaci´on suma de matrices, tiene estructura de grupo 
conmutativo, ya que la suma de matrices es una ley de composici´on interna, es aso- 
ciativa, tiene elemento neutro (la matriz nula), toda matriz tiene elemento sim´etrico 
(su matriz opuesta) y la suma de matrices es conmutativa. 
Se representa por (Mm×n(K); +). 


El conjunto de todas las matrices de tama˜no n × n sobre el cuerpo K lo repre- 
sentamos por Mn(K). 


Ejemplo 3.1.9 
El conjunto Mn(K) con las operaciones suma y producto de matrices, tiene estruc- 
tura de anillo unitario, ya que como hemos visto en el Ejemplo 3.1.8 el conjunto 
Mn(K) con la suma de matrices tiene estructura de grupo conmutativo; adem´as, 
el producto de matrices es asociativo, tiene elemento neutro (la matriz identidad) 
y se cumple la propiedad distributiva del producto respecto de la suma. Por lo tan- 
to, el conjunto Mn(K), con las operaciones suma y producto de matrices, tiene 
estructura de anillo unitario. Se representa por (Mn(K); +, ·). 
En cambio, no tiene estructura de cuerpo porque s´olo las matrices que tienen deter- 
minante distinto de cero (elemento neutro de la suma en K) tienen inversa. 


Consideremos un cuerpo K y n ∈ N ∪ {0}, y sean a0, a1, . . . , an ∈ K, con 
an = 0. Decimos que la aplicaci´on 


p : K −→ K 
x −→ p(x) = a0 + a1 x + · · · + an xn 
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es un polinomio de grado n en la indeterminada x con coeﬁcientes en el cuerpo K; 
se dice que los escalares a1, a2, . . . , an son los coeﬁcientes del polinomio p. 
Utilizamos la notaci´on x2 = x x, x3 = xxx, . . . 
El conjunto de todos los polinomios con coeﬁcientes en el cuerpo K se denota 
por K[x]. 


Ejemplo 3.1.10 
La aplicaci´on p(x) = x3 − 2x + 3 es un polinomio de grado 3 con coeﬁcientes en 
el cuerpo de los n´umeros reales (racionales). 
Evidentemente 
p(x) = 1x3 + 0x2 + (−2)x + 3 


Un polinomio p ∈ K[x] se dice que es un monomio si es de la forma 


p(x) = axn 


donde el escalar a ∈ K, a = 0, se llama coeﬁciente del monomio y el n´umero n se 
llama grado del monomio. 


Ejemplo 3.1.11 
El polinomio p(x) = 4x3 con coeﬁcientes en el cuerpo de los n´umeros reales (ra- 
cionales) es un monomio de grado 3. 


En el conjunto K[x] se puede deﬁnir una ley de composici´on interna 


+ : K[x] × K[x] −→ K[x] 
(p, q) 
−→ p + q 


donde 
(p + q)(x) = p(x) + q(x) 


Ejercicio 3.1.3 
Sea K un cuerpo. Demuestra que (K[x]; +) es un grupo conmutativo. 


En el conjunto K[x] se puede deﬁnir otra ley de composici´on interna 


· : K[x] × K[x] −→ K[x] 
(p, q) 
−→ p · q 


donde 
(p · q)(x) = p(x) q(x) 


Ejercicio 3.1.4 
Sea K un cuerpo. Demuestra que (K[x]; +, ·) es un anillo conmutativo y unitario. 


Es evidente que K[x] no tiene estructura de cuerpo; el monomio p(x) = x no 
tiene inverso, ya que f(x) = 1/x, que es una aplicaci´on deﬁnida en K − {0}, no es 
un polinomio y es la ´unica aplicaci´on que cumple: 


p(x) f(x) = 1 , 
x = 0 
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3.2. 
Espacios vectoriales. El espacio vectorial Rn 


En el Tema 1 hemos hablado de vectores al referirnos a la soluci´on o solucio- 
nes de un sistema de ecuaciones lineales. En esta secci´on vamos a introducir for- 
malmente el concepto de espacio vectorial para despu´es particularizarlo al espacio 
vectorial Rn. 


3.2.1. 
Deﬁnici´on de espacio vectorial 


Deﬁnici´on 3.2.1 (Espacio vectorial) 
Sea E un conjunto no vac´ıo y sea K un cuerpo. Se dice que el conjunto E tiene 
estructura de espacio vectorial sobre el cuerpo K si, y s´olo, si 


1. Existe una ley de composici´on interna, (+), sobre E tal que (E; +) es un 
grupo conmutativo. 


El elemento neutro lo representamos por 0 y lo llamamos vector nulo; dado 
cualquier u ∈ E, al elemento opuesto de u lo representamos por −u. 


A esta operaci´on la llamamos suma. 


2. Existe una ley de composici´on externa, (·), sobre E con dominio de operado- 
res K 
· : K × E −→ E 
(λ, u) −→ λ · u 


cumpliendo las siguientes propiedades 


a) λ · (u + v) = λ · u + λ · v, ∀λ ∈ K, ∀u,v ∈ E. 


b) (λ + µ) · u = λ · u + µ · u, ∀λ, µ ∈ K, ∀u ∈ E. 


c) (λ µ) · u = λ · (µ · u), ∀λ, µ ∈ K, ∀u ∈ E. 


d) 1 · u = u, ∀u ∈ E. 


A esta operaci´on la llamamos producto por un escalar. 



Los elementos de E se llaman vectores y los elementos de K se llaman, como 
indicamos anteriormente, escalares. 
Un espacio vectorial es, por tanto, una terna [(E; +), (K; +, ·), ·] donde (E; +) 
es un grupo conmutativo, (K; +, ·) es un cuerpo y (·) es una ley de composici´on 
externa sobre E con dominio de operadores K, veriﬁcando las propiedades anterio- 
res. Para simpliﬁcar se suele decir que E es un espacio vectorial sobre un cuerpo K 
o que E es un K-espacio vectorial. 
Hay que hacer notar aqu´ı hemos designado con el mismo s´ımbolo, (+), la suma 
de escalares y la suma de vectores, y con el mismo s´ımbolo, (·), el producto de 
escalares y la ley de composici´on externa (producto de un escalar y un vector). 
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Ejemplo 3.2.1 
Todo cuerpo es un espacio vectorial sobre ´el mismo. Para ello basta tomar como ley 
de composici´on externa el producto de escalares. 
En particular, el cuerpo de los n´umeros reales, R, y el cuerpo de los n´umeros com- 
plejos, C, tienen estructura de espacio vectorial, el primero sobre el cuerpo de los 
n´umeros reales y el segundo sobre el cuerpo de los n´umeros complejos (y sobre el 
cuerpo de los n´umeros reales). 


Ejemplo 3.2.2 
El conjunto Mm×n(K) tiene estructura de espacio vectorial sobre el cuerpo K, ya 
que, como hemos visto en el Ejemplo 3.1.8, (Mm×n(K); +) es un grupo conmu- 
tativo y, adem´as, el producto de un escalar y una matriz es una ley de composici´on 
externa sobre Mm×n(K) con dominio de operadores K, que cumple las propieda- 
des: 


(a) λ · (A + B) = λ · A + λ · B, ∀λ ∈ K, ∀A, B ∈ Mm×n(K). 


(b) (λ + µ) · A = λ · A + µ · A, ∀λ, µ ∈ K, ∀A ∈ Mm×n(K). 


(c) (λ µ) · A = λ · (µ · A), ∀λ, µ ∈ K, ∀A ∈ Mm×n(K). 


(d) 1 · A = A, ∀A ∈ Mm×n(K). 



En el conjunto de los polinomios con coeﬁcientes en un cuerpo K, K[x], se puede 
deﬁnir una ley de composici´on externa con dominio de operadores K 


· : K × K[x] 
K[x] 
(λ, p) 
λ · p 


donde 
(λ · p)(x) = λ p(x) 


Ejercicio 3.2.1 
Dado un cuerpo K, demuestra que K[x] tiene estructura de K-espacio vectorial. 


Ejercicio 3.2.2 
Sea E un espacio vectorial sobre un cuerpo K. Entonces, para cualesquiera que 
sean los vectores u,v ∈ E y los escalares λ, µ ∈ K se veriﬁca: 


(a) 0 · u = 0. 


(b) λ ·0 = 0. 


(c) Si λ · u = 0, entonces λ = 0 ´o u = 0. 
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3.2.2. 
El espacio vectorial Rn 


Recordemos que en el Apartado 2.2.1 introdujimos el espacio vectorial R2. Aho- 
ra vamos a formalizar esta estructura algebraica y, a su vez, esta formalizaci´on nos 
servir´a para aﬁanzar el concepto de espacio vectorial. 


Consideremos el conjunto (no vac´ıo) 


R2 = {(x, y) : x, y ∈ R} 


y el cuerpo de los n´umeros reales, R. 
En el conjunto R2 se puede deﬁnir la ley de composici´on interna 


+ : R2 × R2 
R2 


((u1, u2), (v1, v2)) 
(u1 + v1, u2 + v2) 


la cual, evidentemente, dota a R2 de estructura de grupo conmutativo. 
Adem´as, podemos deﬁnir la ley de composici´on externa sobre R2 con dominio 
de operadores R 


· : R × R2 
R2 


(λ , (u, v)) 
(λ u, λ v) 


la cual, evidentemente, junto con la ley de composici´on interna anterior, dota a R2 


de estructura de espacio vectorial sobre el cuerpo R. 


Ejercicio 3.2.3 
Demuestra que R2, con las dos operaciones anteriores, tiene estructura de espacio 
vectorial sobre el cuerpo de los n´umeros reales. 


An´alogamente, dado el conjunto (no vac´ıo) 


R3 = {(x, y, z) : x, y, z ∈ R} 


y el cuerpo de los n´umeros reales, R, podemos deﬁnir una ley de composici´on inter- 
na sobre R3 y una ley de composici´on externa sobre R3 con dominio de operadores 
R, generalizando las operaciones anteriores deﬁnidas en R2, dotando a R3 de es- 
tructura de espacio vectorial sobre el cuerpo R. 


Ejercicio 3.2.4 
Demuestra que R3, con las dos operaciones anteriores, tiene estructura de espacio 
vectorial sobre el cuerpo de los n´umeros reales. 


En general, 


Rn = {(x1, x2, . . . , xn) : xi ∈ R , i = 1, 2, . . . , n} 


con la ley de composici´on interna 


(u1, u2, . . . , un) + (v1, v2, . . . , vn) = (u1 + v1, u2 + v2, . . . , un + vn) 
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y la ley de composici´on externa con dominio de operadores R 


λ · (u1, u2, . . . , un) = (λ u1, λ u2, . . . , λ un) 


tiene estructura de espacio vectorial sobre el cuerpo de los n´umeros reales. 


Ejercicio 3.2.5 
Demuestra que Rn, con las dos operaciones anteriores, tiene estructura de espacio 
vectorial sobre el cuerpo de los n´umeros reales. 


Notar que si K es un cuerpo, entonces al conjunto 


Kn = {(x1, x2, . . . , xn) : xi ∈ K , i = 1, 2, . . . , n} 


se le puede dotar de estructura de K-espacio vectorial de forma an´aloga a Rn. 


3.3. 
Subespacios vectoriales 


Consideremos un espacio vectorial E sobre un cuerpo K y un subconjunto no 
vac´ıo S de E (S ⊆ E). Queremos determinar qu´e condiciones debe cumplir el 
subconjunto S para que con las operaciones deﬁnidas en E, el subconjunto S ten- 
ga estructura de espacio vectorial sobre el cuerpo K. En tal caso, diremos que el 
subconjunto S es un subespacio vectorial de E. 


3.3.1. 
Deﬁnici´on de subespacio vectorial 


Deﬁnici´on 3.3.1 (Subespacio vectorial) 
Sea E un espacio vectorial sobre un cuerpo K y S ⊆ E, S = ∅. 
Decimos que S es un subespacio vectorial de E si S es, con las mismas operaciones 
de E, un espacio vectorial sobre K. 


Teorema 3.3.1 (Caracterizaci´on de subespacio vectorial) 
Sea E un espacio vectorial sobre un cuerpo K y S ⊆ E, S = ∅. 
S es un subespacio vectorial de E si, y s´olo si, dados cualesquiera u,v ∈ S y 
cualesquiera λ, µ ∈ R, entonces λ · u + µ · v ∈ S. 


Ejemplo 3.3.1 
Consideremos el espacio vectorial R3 y veamos que el plano 


S = {(x, y, z) ∈ R3 : x − 4y + z = 0} 


es un subespacio vectorial de R3. 


Para ello tomamos cualquier par de vectores u = (u1, u2, u3) y v = (v1, v2, v3) de 
S y cualquier par de escalares λ, µ ∈ R, y calculamos la combinaci´on lineal: 


λ · u + µ · v = (λu1 + µv1, λu2 + µv2, λu3 + µv3) 
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Pero los vectores u,v ∈ S, por tanto 


u1 − 4u2 + u3 = 0 , 
v1 − 4v2 + v3 = 0 


De aqu´ı 


λu1 + µv1 − 4(λu2 + µv2) + λu3 + µv3 = 
λ(u1 − 4u2 + u3) + µ(v1 − 4v2 + v3) = 0 


Entonces 
λ · u + µ · v ∈ S 


y as´ı demostramos que el plano S es un subespacio vectorial de R3. 


Ejercicio 3.3.1 
Sea K un cuerpo, n ∈ N y 


Kn[x] = {a0 + a1 x + · · · + an xn : ai ∈ K , i = 0, 1, . . . , n} 


el conjunto de polinomios de grado menor o igual que n con coeﬁcientes en el 
cuerpo K. Demuestra que Kn[x] es un subespacio vectorial de K[x]. 


Un espacio vectorial E tiene como subespacios vectoriales al propio espacio 
vectorial E y al conjunto formado por el vector nulo, {0}. Tales subespacios se lla- 
man subespacios impropios de E. Los subespacios propios de E son los subes- 
pacios vectoriales de E distintos de E y {0}. 


Ejemplo 3.3.2 
Consideremos el espacio vectorial R3 y veamos que el plano 


S = {(x, y, z) ∈ R3 : x − 4y + z = 3} 


no es un subespacio vectorial de R3. 


Elegimos los vectores u = (1, 0, 2) y v = (0, 0, 3) de S, y los escalares λ = 1 y 
µ = 1. Entonces 


λ · u + µ · v = 1 · (1, 0, 2) + 1 · (0, 0, 3) = (1, 0, 5) /∈ S 


por tanto, el plano S no es un subespacio vectorial de R3. 


Como veremos m´as adelante, los ´unicos subespacios propios de R3 son las rec- 
tas y los planos que pasan por el origen. Tambi´en comprobaremos que, en general, 
los ´unicos subespacios propios de Rn son los correspondientes a las soluciones de 
sistemas de ecuaciones lineales con n inc´ognitas homog´eneos compatibles indeter- 
minados. 
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3.3.2. 
Suma e intersecci´on de subespacios 


Teorema 3.3.2 (Suma de subespacios) 
Sean S, T dos subespacios vectoriales de un K-espacio vectorial E. 
Llamamos suma de S y T al conjunto 


S + T = {u + v : u ∈ S , v ∈ T} 


La suma de S y T es un subespacio vectorial de E. 


Ejercicio 3.3.2 
Demuestra el Teorema 3.3.2. 


Teorema 3.3.3 (Intersecci´on de subespacios) 
Sean S, T dos subespacios vectoriales de un K-espacio vectorial E. 
La intersecci´on de S y T 


S ∩ T = {u ∈ E : u ∈ S , u ∈ T} 


es un subespacio vectorial de E. 


Ejercicio 3.3.3 
Demuestra el Teorema 3.3.3. 


En general, la uni´on de subespacios vectoriales no es un subespacio vectorial. 


Ejemplo 3.3.3 
Sean los subespacios vectoriales de R3 


S = {(x, y, z) ∈ R3 : x − y = 0} , 
T = {(x, y, z) ∈ R3 : x + z = 0} 


Veamos que la uni´on de S y T, S ∪ T, no es un subespacio vectorial de R3. 


Para ello consideramos los vectores 


u = (1, 1, 0) ∈ S ⊂ S ∪ T , 
v = (1, 0, −1) ∈ T ⊂ S ∪ T 


y los escalares λ = 1 y µ = 1. Entonces 


λ · u + µ · v = 1 · (1, 1, 0) + 1 · (1, 0, −1) = (2, 1, −1) /∈ S ∪ T 


ya que (2, 1, −1) /∈ S y (2, 1, −1) /∈ T. 


Deﬁnici´on 3.3.2 (Suma directa) 
Sean S, T dos subespacios vectoriales de un K-espacio vectorial E. 
Decimos que la suma de S y T es directa, y lo denotamos por S ⊕ T, si cada vector 
de S + T se puede expresar de forma ´unica como la suma de un vector de S y un 
vector de T. 
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Teorema 3.3.4 (Caracterizaci´on de suma directa) 
Sean S, T dos subespacios vectoriales de un K-espacio vectorial E. 
La suma de S y T es directa si, y s´olo si, S ∩ T = {0}. 


Ejemplo 3.3.4 
Dados los subespacios vectoriales de R4 


S = {(x, y, z, t) ∈ R4 : x = y = z = t} 
T = {(x, y, z, t) ∈ R4 : y = 2x , t = 2z} 


veamos que la suma de S y T es directa. 


La intersecci´on de S y T es 


S ∩ T = {(x, y, z, t) ∈ R4 : x = y = z = t , y = 2x , t = 2z} 


es decir, sus vectores son soluci´on del sistema de ecuaciones lineales 
 
 
 
 
 
 


 
 
 
 
 


x − y = 0 
x − z = 0 
x − t = 0 
2x − y = 0 
2z − t = 0 


Escalonando la matriz de coeﬁcientes del sistema, obtenemos 
 


 
 
 
 
 


1 
−1 
0 
0 
1 
0 
−1 
0 
1 
0 
0 
−1 
2 
−1 
0 
0 
0 
0 
2 
−1 


 


 
 
 
 
 
∼ 


 


 
 
 
 
 


1 
−1 
0 
0 
0 
1 
−1 
0 
0 
1 
0 
−1 
0 
1 
0 
0 
0 
0 
2 
−1 


 


 
 
 
 
 
∼ 


 


 
 
 
 
 


1 
−1 
0 
0 
0 
1 
−1 
0 
0 
0 
1 
−1 
0 
0 
1 
0 
0 
0 
2 
−1 


 


 
 
 
 
 


∼ 


 


 
 
 
 
 


1 
−1 
0 
0 
0 
1 
−1 
0 
0 
0 
1 
−1 
0 
0 
0 
1 
0 
0 
0 
1 


 


 
 
 
 
 
∼ 


 


 
 
 
 
 


1 
−1 
0 
0 
0 
1 
−1 
0 
0 
0 
1 
−1 
0 
0 
0 
1 
0 
0 
0 
0 


 


 
 
 
 
 


es decir, el sistema anterior es equivalente a 
 
 
 
 


 
 
 


x − y = 0 
y − z = 0 
z − t = 0 
t = 0 


el cual, tiene como ´unica soluci´on t = 0, z = 0, y = 0, x = 0, esto es, la soluci´on 
trivial (0, 0, 0, 0). Por tanto 


S ∩ T = {(0, 0, 0, 0)} 


es decir, la suma de S y T es directa. 
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3.4. 
Dependencia e independencia lineal. Sistema ge- 
nerador de un subespacio 


En este apartado introducimos el concepto de sistema de vectores y los clasi- 
ﬁcamos en sistemas libres y sistemas ligados. Adem´as, mediante la combinaci´on 
lineal de vectores, podemos deﬁnir la noci´on de subespacio vectorial generado por 
un sistema de vectores, que ser´a una de las dos formas en la que expresaremos un 
subespacio vectorial. 


3.4.1. 
Dependencia e independencia lineal 


Deﬁnici´on 3.4.1 (Sistema de vectores) 
Un sistema de vectores de un K-espacio vectorial E es una colecci´on de vectores 
de E, los cuales se pueden repetir. 


Ejemplo 3.4.1 
En el espacio vectorial R3, el conjunto S = {(1, −1, 0), (−2, 3, 1)} es un sistema 
de vectores de R3. 
En el espacio vectorial M2×3(R), el conjunto de matrices 
1 
−1 
3 
2 
−4 
1 


, 
5 
−2 
0 
1 
7 
−1 


, 
−3 
5 
−6 
1 
0 
9 



es un sistema de vectores de M2×3(R). 


Deﬁnici´on 3.4.2 (Combinaci´on lineal) 
Sea E un K-espacio vectorial. Dado un vector u ∈ E, diremos que u es combina- 
ci´on lineal del sistema de vectores S = {v1,v2, . . . ,vn} (combinaci´on lineal de los 
vectores v1,v2, . . . ,vn), si existen n escalares λ1, λ2, . . . , λn tales que 


u = λ1 · v1 + λ2 · v2 + · · · + λn · vn 


Los escalares λ1, λ2, . . . , λn se llaman coeﬁcientes de la combinaci´on lineal. 


Ejemplo 3.4.2 
En el espacio vectorial R3, el vector (−1, 5, 0) es combinaci´on lineal de los vectores 
del sistema 
S = {(1, 1, −1), (2, 0, 1), (0, 2, 1)} 


ya que 
(−1, 5, 0) = 1 · (1, 1, −1) + (−1) · (2, 0, 1) + 2 · (0, 2, 1) 


y, por lo tanto, los coeﬁcientes de la combinaci´on lineal son 1, −1 y 2. 


Notar aqu´ı que la Deﬁnici´on 3.4.2 generaliza el concepto de combinaci´on lineal 
de l´ıneas de una matriz utilizado en el apartado 1.3.2 al considerar, por ejemplo, las 
ﬁlas de una matriz de tama˜no 4 × 3 como 4 vectores de R3. 
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Deﬁnici´on 3.4.3 (Sistema ligado) 
En un K-espacio vectorial E, un sistema de vectores S = {v1,v2, . . . ,vn} es ligado 
(sus vectores son linealmente dependientes), si existen n escalares λ1, λ2, . . . , λn, 
no todos nulos, tales que 


λ1 · v1 + λ2 · v2 + · · · + λn · vn = 0 


Si S es un sistema formado por un n´umero inﬁnito de vectores, se dice que el siste- 
ma S es ligado si S tiene un subconjunto ﬁnito de vectores ligado. 


Ejemplo 3.4.3 
En el espacio vectorial R3, el sistema de vectores 


S = {(1, 1, 1), (1, 0, 1), (0, 1, 0)} 


es ligado porque 


1 · (1, 1, 1) + (−1) · (1, 0, 1) + (−1) · (0, 1, 0) = (0, 0, 0) 



Este ejemplo nos muestra que cuando un sistema es ligado, un vector del sistema 
siempre se puede poner como combinaci´on lineal del resto de vectores del sistema; 
en este caso 


(1, 1, 1) = (1, 0, 1) + (0, 1, 0) = 1 · (1, 0, 1) + 1 · (0, 1, 0) 


Ejercicio 3.4.1 
Sean E un K-espacio vectorial S un sistema de vectores de E. Si 0 ∈ S, entonces 
el sistema S es ligado. 


Deﬁnici´on 3.4.4 (Sistema libre) 
En un K-espacio vectorial E, un sistema de vectores S = {v1,v2, . . . ,vn} es libre 
(sus vectores son linealmente independientes), si cualquier combinaci´on lineal 
igualada al vector nulo 


λ1 · v1 + λ2 · v2 + · · · + λn · vn = 0 


implica necesariamente que λ1 = λ2 = · · · = λn = 0. 
Si S es un sistema formado por un n´umero inﬁnito de vectores, se dice que el siste- 
ma S es libre si todo subconjunto ﬁnito de S es libre. 


Ejemplo 3.4.4 
En el espacio vectorial R3, el sistema de vectores 


S = {(1, 1, −1), (2, 0, 1), (0, 2, 1)} 
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es libre, ya que si consideramos cualquier combinaci´on lineal de los vectores del 
sistema S igualada al vector nulo 


λ · (1, 1, −1) + µ · (2, 0, 1) + ν · (0, 2, 1) = (0, 0, 0) 


obtenemos que los coeﬁcientes de la combinaci´on lineal son soluci´on del sistema 
 
 


 


λ + 2µ 
= 0 
λ + 
+ 2ν = 0 
−λ + 
µ + 
ν = 0 


La matriz de coeﬁcientes de este sistema es 


A = 


 


 
1 
2 
0 
1 
0 
2 
−1 
1 
1 


 


 


Observemos que la columnas de esta matriz corresponden a las componentes de los 
vectores del sistema S. 
Es f´acil comprobar que la matriz A tiene rango igual a 3, por lo tanto, aplicando el 
teorema de Rouch´e-Fr¨obenius, la ´unica soluci´on del sistema de ecuaciones lineales 
anterior es la soluci´on trivial λ = 0, µ = 0, ν = 0. De aqu´ı, el sistema S es libre. 


Ejercicio 3.4.2 
En el espacio vectorial R4, determina si los siguientes sistemas de vectores son 
libres o ligados: 


(a) S = {(1, −1, 1, −1), (2, −1, 3, 4), (3, −2, 4, 3)}. 


(b) T = {(1, 3, 2, −1), (0, 1, −2, 1), (−1, 7, −2, −4), (1, 0, 0, 1)}. 



Al resolver este ejercicio tomando como base el Ejemplo 3.4.4, se puede inferir 
que un sistema formado por m vectores en el espacio vectorial Rn es libre si, y s´olo 
si, al poner esos vectores como ﬁlas o columnas de una matriz (de tama˜no m×n o de 
tama˜no n×m), esa matriz tiene rango m. Por tanto, en el espacio vectorial Rn todo 
sistema formado por m´as de n vectores es ligado; esta consecuencia ser´a enunciada 
en general cuando introduzcamos los espacios vectoriales de tipo ﬁnito. 


Ejercicio 3.4.3 
Sea K un cuerpo. Demuestra que el sistema de vectores 


S = {1, x, x2, x3, . . .} 


del espacio vectorial K[x] es un sistema libre. 
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Ejemplo 3.4.5 
El sistema de vectores de Rn 


C = {(1, 0, . . . , 0), (0, 1, . . ., 0), . . . , (0, 0, . . . , 1)} 


es un sistema libre, ya que 


1 
0 
· · · 
0 
0 
1 
· · · 
0 
... 
... 
... 
... 
0 
0 
· · · 
1 



= 1 = 0 


Es evidente que este sistema de vectores est´a formado por n vectores. 


Deﬁnici´on 3.4.5 (Rango de un sistema de vectores) 
Se llama rango de un sistema de vectores al n´umero m´aximo de vectores lineal- 
mente independientes del sistema. 


Como consecuencia de lo expuesto anteriormente, el rango de un sistema de 
vectores coincide con el rango de la matriz que se obtiene al poner esos vectores 
como ﬁlas o columnas. 


Ejemplo 3.4.6 
El rango del sistema de vectores 


S = {(1, −1, 2, 3), (2, 1, 0, −1), (3, 0, 2, 2)} 


es igual a 2, ya que 


rang 


 


 
1 
−1 
2 
3 
2 
1 
0 
−1 
3 
0 
2 
2 


 


 = rang 


 


 
1 
−1 
2 
3 
0 
−3 
4 
7 
0 
−3 
4 
7 


 


 


= rang 


 


 
1 
−1 
2 
3 
0 
−3 
4 
7 
0 
0 
0 
0 


 


 = 2 



Ejercicio 3.4.4 
Calcula el rango de los siguientes sistemas de vectores 


(a) {(1, 1, −1), (2, 0, 1), (0, 2, 1)} 


(b) {(1, 1, 1), (1, 0, 1), (0, 1, 0)} 


(c) {3, 1, 4), (2, −1, 1), (−5, 3, 2)} 


(d) {(1, 3, 2), (5, 4, −1), (−3, 0, 3)} 


(e) {(1, −1, 1, −1), (2, −1, 3, 4), (3, −2, 4, 3)} 


(f ) {(1, 3, 2, −1), (0, 1, −2, 1), (−1, 7, −2, −4), (1, 0, 0, 1)} 
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3.4.2. 
Sistema generador de un subespacio 


Teorema 3.4.1 (Subespacio generado por un sistema de vectores) 
Sea E un K-espacio vectorial. Dado un sistema de vectores S = {v1,v2, . . . ,vn}, 
el conjunto de combinaciones lineales que se pueden hacer con los vectores de S 
es un subespacio vectorial de E, que denotamos por S = v1,v2, . . . ,vn y que 
llamamos subespacio generado por S, es decir 


S = {λ1 v1 + · · · + λn vn : λ1, . . . , λn ∈ K} 


Si S es un sistema formado por un n´umero inﬁnito de vectores, el subespacio gene- 
rado por S es 


S = {λ1 v1 + · · · + λn vn : λ1, . . . , λn ∈ K , v1, . . . ,vn ∈ S, n ∈ N} 


Adem´as, S es el menor subespacio vectorial que contiene a S. 


Ejercicio 3.4.5 
Demuestra el Teorema 3.4.1. 


Deﬁnici´on 3.4.6 (Sistema generador de un subespacio) 
Decimos que un sistema de vectores {v1,v2, . . . ,vn} de un K-espacio vectorial E 
es un sistema generador de un subespacio vectorial S de E, si 


S = v1,v2, . . . ,vn 


Si T es un sistema formado por un n´umero inﬁnito de vectores de E, decimos que 
T es un sistema generador de un subespacio vectorial S de E si S = T. 


Ejemplo 3.4.7 
Veamos que el sistema de vectores S = {(1, 2, −1), (−1, 0, 2), (1, 2, 1)} es un sis- 
tema generador de R3. 


Por la Deﬁnici´on 3.4.6 hemos de demostrar que R3 = S. 
Esta igualdad entre estos dos conjuntos la vamos a comprobar viendo que un con- 
junto es un subconjunto del otro conjunto y viceversa: 


⊆ Sea (x, y, z) ∈ R3. Queremos probar que (x, y, z) ∈ S. Pero (x, y, z) ∈ S 
si, y s´olo si, existen tres escalares λ, µ, ν ∈ R tales que 


(x, y, z) = λ · (1, 2, −1) + µ · (−1, 0, 2) + ν · (1, 2, 1) 


Entonces, (x, y, z) ∈ S si, y s´olo si, el sistema 
 
 


 


λ − 
µ + 
ν = x 
2λ 
+ 2ν = y 
−λ + 2µ + 
ν = z 
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es compatible. Pero el determinante de la matriz de coeﬁcientes de este siste- 
ma es 


|A| = 



1 
−1 
1 
2 
0 
2 
−1 
2 
1 


= 4 = 0 


Aplicando el teorema de Rouch´e-Fr¨obenius, el sistema anterior es compatible 
determinado. Su ´unica soluci´on es 


λ = −x + 3 


4y − 1 


2z , 
µ = −x + 1 


2y , 
ν = x − 1 


4y + 1 


2z 


Por tanto, R3 ⊆ S. 


⊇ Sea (x, y, z) ∈ S. Entonces, existen λ, µ, ν ∈ R tal que 


(x, y, z) = λ · (1, 2, −1) + µ · (−1, 0, 2) + ν · (1, 2, 1) 
= (λ − µ + ν, 2λ + 2ν, −λ + 2µ + ν) ∈ R3 


es decir, S ⊆ R3. 



Ejercicio 3.4.6 
Demuestra que el sistema de vectores 


C = {(1, 0, . . . , 0), (0, 1, . . ., 0), . . . , (0, 0, . . . , 1)} 


es un sistema generador de Rn, vi´endolo primero para n = 2 y n = 3, y despu´es en 
general para cualquier n ∈ N. 


Ejercicio 3.4.7 
Demuestra que 
M2×3(R) = A1, A2, A3, A4, A5, A6 


donde 


A11 = 
1 
0 
0 
0 
0 
0 


, 
A12 = 
0 
1 
0 
0 
0 
0 


, 
A13 = 
0 
0 
1 
0 
0 
0 


, 


A21 = 
0 
0 
0 
1 
0 
0 


, 
A22 = 
0 
0 
0 
0 
1 
0 


, 
A23 = 
0 
0 
0 
0 
0 
1 




En general, se puede demostrar que 


S = {A11, A12, . . . , A1n, A21, A22, . . . , A2n, . . . , Am1, Am2, . . . , Amn} 


donde Aij es la matriz de tama˜no m × n, tal que tiene todos sus elementos nulos 
salvo el que ocupa el lugar ij que es igual a 1, es un sistema generador del espacio 
vectorial Mm×n(R). 
En el ejercicio anterior, el cuerpo de los n´umeros reales se puede sustituir por 
un cuerpo arbitrario K. 
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Ejercicio 3.4.8 
Sea K un cuerpo. Demuestra que el sistema de vectores 


S = {1, x, x2, x3, . . .} 


es un sistema generador del espacio vectorial K[x]. 


Ejercicio 3.4.9 
Sean E un K-espacio vectorial y 


S = u1, u2, . . . , um , 
T = v1,v2, . . . ,vn 


donde ui ∈ E, i = 1, 2, . . . , m, y vj ∈ E, j = 1, 2, . . . , n. 
Entonces 
S + T = u1, u2, . . . , um,v1,v2, . . . ,vn 



Teorema 3.4.2 
Sean E un K-espacio vectorial y S = v1,v2, . . . ,vn el subespacio generado por 
los vectores v1, v2, . . . , vn de E. Si el vector vi, para alg´un i ∈ {1, 2, . . . , n}, es 
combinaci´on lineal de los vectores v1,v2, . . . ,vi−1,vi+1, . . . ,vn, entonces 


S = v1,v2, . . . ,vi−1,vi+1, . . . ,vn 



Ejemplo 3.4.8 
El subespacio generado por el sistema de vectores dado en el Ejemplo 3.4.6 es 


S = (1, −1, 2, 3), (2, 1, 0, −1) 


ya que al calcular el rango en el Ejemplo 3.4.6 se puede comprobar tambi´en que el 
tercer vector del sistema S, (3, 0, 2, 2), es combinaci´on lineal de los otros dos. 
Adem´as 
S = (1, −1, 2, 3), (0, −3, 4, 7) 



3.5. 
Bases y dimensi´on. Cambios de base 


El concepto de base es uno de los m´as importantes dentro de la estructura de 
espacio vectorial. Aqu´ı s´olo vamos a considerar espacios vectoriales de dimensi´on 
ﬁnita, los cuales son los m´as comunes en Ingenier´ıa. 
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3.5.1. 
Bases y dimensi´on 


Deﬁnici´on 3.5.1 (Espacio vectorial de tipo ﬁnito) 
Un espacio vectorial es de tipo ﬁnito, si admite un sistema generador con un n´umero 
ﬁnito de vectores. 


Ejemplo 3.5.1 
El espacio vectorial R3 es, como hemos visto en el Ejemplo 3.4.4, de tipo ﬁnito. 


Ejercicio 3.5.1 
Demuestra que el espacio vectorial Rn es de tipo ﬁnito. 


Ejemplo 3.5.2 
El espacio vectorial M2×3(R) es un espacio vectorial de tipo ﬁnito, como conse- 
cuencia del Ejercicio 3.4.6. 
En general, y como consecuencia de la observaci´on posterior al Ejercicio 3.4.6, el 
espacio vectorial Mm×n(R) es de tipo ﬁnito. 


Ejercicio 3.5.2 
Sea K un cuerpo. Demuestra que el espacio vectorial K[x] no es de tipo ﬁnito. 


A partir de aqu´ı, salvo que indiquemos lo contrario, s´olo consideraremos espa- 
cios vectoriales de tipo ﬁnito. 


Deﬁnici´on 3.5.2 (Base de un espacio vectorial) 
Un sistema de vectores de un K-espacio vectorial E se dice que es una base de E, 
si es un sistema libre y generador de E. 


Teorema 3.5.1 (Existencia de base) 
Todo K-espacio vectorial E de tipo ﬁnito, E = {0}, admite una base formada por 
un n´umero ﬁnito de vectores. 
En general, todo K-espacio vectorial E , E = {0}, admite una base. 


Ejemplo 3.5.3 
Como consecuencia del Ejemplo 3.4.5 y del Ejercicio 3.4.6, el sistema de vectores 


C = {(1, 0, . . . , 0), (0, 1, . . ., 0), . . . , (0, 0, . . . , 1)} 


es una base de Rn. 
Esta base de Rn se llama base can´onica de Rn. 
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Ejercicio 3.5.3 
Comprueba que el sistema de vectores 


B = {(2, −1, 3), (1, 4, 0), (3, 2, −1)} 


es una base de R3. 


Teorema 3.5.2 (Caracterizaci´on de base) 
Sea E un K-espacio vectorial y B un sistema de vectores de E. 
B es una base de E si, y s´olo si, todo vector de E se puede expresar de forma ´unica 
como combinaci´on lineal de vectores de B. 


Ejercicio 3.5.4 
Demuestra el Teorema 3.5.2. 


Deﬁnici´on 3.5.3 
Sean E un K-espacio vectorial, B = {v1,v2, . . . ,vn} una base de E y v ∈ E. 
Entonces, en virtud del Teorema 3.5.2, existen n ´unicos escalares λ1, λ2, . . . , λn 
tales que 
v = λ1 · v1 + λ2 · v2 + · · · + λn · vn 


Estos n ´unicos escalares se llaman coordenadas del vector v en la base B. 
En ocasiones las coordenadas se suelen escribir como un vector de Kn. Entonces 
decimos que (λ1, λ2, . . . , λn) ∈ Kn son las coordenadas del vector v en la base B. 


N´otese que la deﬁnici´on anterior requiere que cuando se consideren las coorde- 
nadas de un vector v de un espacio vectorial E con respecto a una base B de E, ´esta 
debe estar dada en un determinado orden. 


Ejercicio 3.5.5 
Demuestra que las coordenadas del vector (4, −12, 5) en la base B del Ejerci- 
cio 3.5.3 son (2, −3, 1). 


Teorema 3.5.3 
Si un espacio vectorial tiene una base de n vectores, entonces: 


1. Todo sistema libre de n vectores es base. 


2. Todo sistema generador de n vectores es base. 



Teorema 3.5.4 
En un espacio vectorial todas las bases tienen el mismo n´umero de vectores. 


Como consecuencia de este teorema, obtenemos la siguiente 
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Deﬁnici´on 3.5.4 (Dimensi´on de un espacio vectorial) 
Sea E un K-espacio vectorial, E = {0}. Llamamos dimensi´on de E, y lo denota- 
mos por dim (E), al n´umero de vectores de una cualquiera de sus bases. 


Convenimos en decir que 
dim ({0}) = 0 


Ejemplo 3.5.4 
El espacio vectorial Rn es un espacio vectorial de dimensi´on n, ya que como hemos 
visto en el Ejemplo 3.5.3, el espacio vectorial Rn tiene una base formada por n 
vectores. 


Teorema 3.5.5 
Sea S un subespacio vectorial de un K-espacio vectorial E. Entonces 


dim (S) dim (E) 


Adem´as, 
dim (S) = dim (E) ⇐⇒ S = E 



Ejercicio 3.5.6 
Dados los siguientes sistemas de vectores 


(a) {(1, 0, 0), (0, 1, 1), (1, 1, 1)}. 


(b) {(1, 0, 0, 4), (3, 2, 1, 0)}. 


(c) {(1, 1, 1, 1), (2, 0, 1, 0)}. 


(d) {(1, 0, 1, 0), (0, 2, 1, 0), (2, 1, 2, 1)}. 


(e) {(1, 2, 3), (0, 1, 2), (1, 3, 5)}. 


(f ) {(1, 2, 3), (0, 1, 2), (1, 3, 5), (1, 2, 2)}. 


(g) {(1, 1, 0), (1, 1, 2)}. 


halla una base y la dimensi´on de los subespacios generados por tales sistemas de 
vectores. 


Como consecuencia del ´ultimo teorema, en un espacio vectorial de dimensi´on 
n, todo sistema de m´as de n vectores es ligado. 


Ejercicio 3.5.7 
Sea E un K-espacio vectorial. Sean B = {u1, u2, . . . ,un} y V = {v1,v2, . . . ,vm} 
dos bases de E. Razona si las siguientes aﬁrmaciones son correctas o si son falsas: 


(a) u1 = vi, para alg´un i ∈ {1, 2, . . . , m}. 
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(b) El vector 0 pertenece a las dos bases. 


(c) n = m. 


(d) u3 no depende linealmente de v1,v2, . . . ,vm. 



Ejercicio 3.5.8 
Sea E un K-espacio vectorial de dimensi´on n. ¿Cu´ales de las siguientes aﬁrmacio- 
nes son correctas y cu´ales falsas? ¿Por qu´e? 


(a) El espacio vectorial E s´olo tiene n vectores y son linealmente independientes. 


(b) En el espacio vectorial E s´olo hay una base formada por n vectores. 


(c) Existen n vectores linealmente independientes v1,v2, . . . ,vn ∈ E tales que 
E = v1,v2, . . . ,vn. 


(d) Cualquier sistema de vectores del espacio vectorial E que contenga menos de 
n vectores nunca es libre. 


(e) Cualquier sistema de vectores del espacio vectorial E que contenga m´as de n 
vectores es siempre ligado. 



Tal y como hemos indicado al ﬁnal del apartado 3.3.1, los ´unicos subespacios 
vectoriales propios de Rn son los correspondientes a las soluciones de sistemas de 
ecuaciones lineales con n inc´ognitas homog´eneos y compatibles indeterminados. 
Estas ecuaciones lineales son las llamadas ecuaciones del subespacio vectorial. 
Veamos esta aﬁrmaci´on en general y despu´es consideraremos un ejemplo. 
Supongamos que S es un subespacio vectorial propio de Rn de dimensi´on m; 
por lo tanto, m < n. Entonces, existen m vectores v1,v2, . . . ,vm ∈ S linealmente 
independientes tales que 
S = v1,v2, . . . ,vm 


Supongamos que 


v1 = (v11, v12, . . . , v1n) , 
v2 = (v21, v22, . . . , v2n) , . . . , 
vm = (vm1, vm2, . . . , vmn) 


Sea x = (x1, x2, . . . , xn) ∈ Rn y consideremos las matrices 


A = 


 


 
 
 
 
 
 


x1 
x2 
· · · 
xn 
v11 
v12 
. . . 
v1n 
v21 
v22 
. . . 
v2n 
... 
... 
... 
... 
vm1 
vm2 
. . . 
vmn 


 


 
 
 
 
 
 
, 
B = 


 


 
 
 
 


v11 
v12 
. . . 
v1n 
v21 
v22 
. . . 
v2n 
... 
... 
... 
... 
vm1 
vm2 
. . . 
vmn 


 


 
 
 
 
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Evidentemente, la matriz B es una submatriz de A. Adem´as, rang(B) = m, ya que 
los vectores v1,v2, . . . ,vm son linealmente independientes. 
Por tanto 


x ∈ S ⇐⇒ ∃λ1, λ2, . . . , λm ∈ R : x = λ1 · v1 + λ2 · v2 + · · · λm · vm 
⇐⇒ rang(A) = m 


Eligiendo un menor de orden m distinto de cero de la matriz B, que tambi´en es 
un menor de orden m no nulo de la matriz A, obtenemos que x ∈ S si, y s´olo si, 
los n − m menores orlados de la matriz A, a partir del menor de orden m no nulo, 
deben ser todos iguales a cero. ´Estas son, pues, las n−m ecuaciones del subespacio 
vectorial S, las cuales son independientes. 
Es decir, si S es un subespacio vectorial propio de Rn, entonces 


dim (S) = n − n´umero de ecuaciones independientes de S 


Es evidente que el espacio vectorial Rn, consider´andolo como subespacio vec- 
torial de ´el mismo, no tiene ecuaciones. 


Ejemplo 3.5.5 
Vamos a calcular una base, la dimensi´on y las ecuaciones independientes del subes- 
pacio vectorial de R4: 


S = (2, 1, 0, −1), (−1, 2, 1, −2), (−4, 4, 1, −3) 


En primer lugar, debemos ver si el sistema generador de S est´a formado por vectores 
linealmente independientes. Para ello, consideramos la matriz cuyas ﬁlas son las 
componentes de esos vectores y la escalonamos 
 


 
2 
6 
1 
−1 
−1 
2 
1 
−2 
−4 
4 
1 
−3 


 


 ∼ 


 


 
2 
6 
1 
−1 
0 
4 
3 
−5 
0 
4 
3 
−5 


 


 ∼ 


 


 
2 
6 
1 
−1 
0 
4 
3 
−5 
0 
0 
0 
0 


 


 


Entonces 
S = (2, 0, 1, −1), (0, 4, 3, −5) 


Adem´as 
BS = {(2, 0, 1, −1), (0, 4, 3, −5)} 


es una base de S y, por tanto, dim(S) = 2. 
En segundo lugar 


(x, y, z, t) ∈ S ⇐⇒ rang 


 


 
x 
y 
z 
t 
2 
0 
1 
−1 
0 
4 
3 
−5 


 


 = 2 


Como 
2 
0 
0 
4 


= 8 = 0 
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se tiene que 


(x, y, z, t) ∈ S ⇐⇒ 



x 
y 
z 
2 
0 
1 
0 
4 
3 


= 0 , 



x 
y 
t 
2 
0 
−1 
0 
4 
−5 


= 0 


Pero 


x 
y 
z 
2 
0 
1 
0 
4 
3 


= 8z − (4x + 6y) = 8z − 4x − 6y = −2(2x + 3y − 4z) 



x 
y 
t 
2 
0 
−1 
0 
4 
−5 


= 8t − (−4x − 10y) = 8t + 4x + 10y = 2(2x + 5y + 4t) 


Por tanto 


S = {(x, y, z, t) ∈ R4 : 2x + 3y − 4z = 0 , 2x + 5y + 4t = 0} 



Ejercicio 3.5.9 
Calcula las ecuaciones de los subespacios generados por los sistemas de vectores 
del Ejercicio 3.5.6. 


Teorema 3.5.6 (Dimensi´on de la suma y la intersecci´on) 
Sean S y T dos subespacios vectoriales de un K-espacio vectorial E. Entonces 


dim (S + T) = dim (S) + dim (T) − dim (S ∩ T) 



Ejemplo 3.5.6 
Dados los subespacios vectoriales de R4: 


S = {(x, y, z, t) ∈ R4 : x − 2y + 3z − t = 0 , 
2x + y − z + 2t = 0 , 3x − y + 2z + t = 0} 


T = (0, 1, 1, −1), (0, 1, 0, 1), (0, 3, 1, 1) 


vamos a calcular una base y las ecuaciones de S + T y de S ∩ T. 


En primer lugar, vamos a obtener las ecuaciones independientes de S y una base. 
Para ello escalonamos la matriz de coeﬁcientes del sistema que determina el subes- 
pacio vectorial S: 
 


 
1 
−2 
3 
−1 
2 
1 
−1 
2 
3 
−1 
2 
1 


 


 ∼ 


 


 
1 
−2 
3 
−1 
0 
−5 
7 
−4 
0 
−5 
7 
−4 


 


 ∼ 


 


 
1 
−2 
3 
−1 
0 
−5 
7 
−4 
0 
0 
0 
0 


 


 
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Por tanto 


S = {(x, y, z, t) ∈ R4 : x − 2y + 3z − t = 0 , −5y + 7z − 4t = 0} 


Para obtener una base de S tenemos que resolver el sistema 
x − 2y + 3z − 
t = 0 
− 5y + 7z − 4t = 0 


La matriz de coeﬁcientes del sistema es 
1 
−2 
3 
−1 
0 
−5 
7 
−4 



y tiene un menor de orden 2 no nulo 
1 
−2 
0 
−5 


= −5 = 0 


Por tanto, el sistema anterior se puede transformar en el sistema de Cramer (de 2 
ecuaciones, 2 inc´ognitas, x e y, y 2 par´ametros, z y t): 
x − 2y + = t − 3z 
− 5y + = 4t − 7z 


cuya soluci´on es 


x = 


t − 3z 
−2 
4t − 7z 
−5 



1 
−2 
0 
−5 



= −5(t − 3z) + 2(4t − 7z) 


−5 
= −z 


5 − 3 


5 t 


y = 


1 
t − 3z 
0 
4t − 7z 



1 
−2 
0 
−5 



= 4t − 7z 


−5 
= 7 


5 z − 4 


5 t 


Entonces 


S = {(x, y, z, t) ∈ R4 : x = −1/5 z − 3/5 t , y = 7/5 z − 4/5 t} 
= {(−1/5 z − 3/5 t, 7/5 z − 4/5 t, z, t) : z, t ∈ R} 
= {z(−1/5, 7/5, 1, 0) + t(−3/5, −4/5, 0, 1) : z, t ∈ R} 
= (−1/5, 7/5, 1, 0), (−3/5, −4/5, 0, 1) 
= (−1, 7, 5, 0), (−3, −4, 0, 5) 


es decir 
S = (−1, 7, 5, 0), (−3, −4, 0, 5) 


Como el sistema de vectores 


BS = {(−1, 7, 5, 0), (−3, −4, 0, 5)} 
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es libre, obtenemos que BS es una base de S. De aqu´ı, dim(S) = 2. 
En segundo lugar, vamos a obtener una base de T y sus ecuaciones independientes. 
Para ello escalonamos la matriz cuyas ﬁlas son los vectores que generan el subes- 
pacio vectorial T: 
 


 
0 
1 
1 
−1 
0 
1 
0 
1 
0 
3 
1 
1 


 


 ∼ 


 


 
0 
1 
1 
−1 
0 
0 
−1 
2 
0 
0 
2 
−4 


 


 ∼ 


 


 
0 
1 
1 
−1 
0 
0 
−1 
2 
0 
0 
0 
0 


 


 


Entonces 
T = (0, 1, 1, −1), (0, 0, −1, 2) 


y el sistema de vectores 


BT = {(0, 1, 1, −1), (0, 0, −1, 2)} 


es un sistema libre, por tanto BT es una base de T y, as´ı, dim(T) = 2. 
Adem´as 


(x, y, z, t) ∈ T ⇐⇒ rang 


 


 
x 
y 
z 
t 
0 
1 
1 
−1 
0 
0 
−1 
2 


 


 = 2 


Como 
1 
1 
0 
−1 


= −1 = 0 


es un menor de orden 2 no nulo, se tiene que 


(x, y, z, t) ∈ S ⇐⇒ 



x 
y 
z 
0 
1 
1 
0 
0 
−1 


= 0 , 



x 
y 
t 
0 
1 
−1 
0 
−1 
2 


= 0 


Pero 


x 
y 
z 
0 
1 
1 
0 
0 
−1 


= −x 



x 
y 
t 
0 
1 
−1 
0 
−1 
2 


= 2y − t − (y + 2z) = y − 2z − t 


Por tanto 
T = {(x, y, z, t) ∈ R4 : x = 0 , y − 2z − t = 0} 


En tercer lugar, debido al Ejercicio 3.4.9 obtenemos que 


S + T = (−1, 7, 5, 0), (−3, −4, 0, 5), (0, 1, 1, −1), (0, 0, −1, 2) 
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Para obtener una base de S + T escalonamos la matriz cuyas ﬁlas son los vectores 
del sistema generador de S + T: 
 


 
 
 


−1 
7 
5 
0 
−3 
−4 
0 
5 
0 
1 
1 
−1 
0 
0 
−1 
2 


 


 
 
 ∼ 


 


 
 
 


−1 
7 
5 
0 
0 
25 
15 
−5 
0 
1 
1 
−1 
0 
0 
−1 
2 


 


 
 
 ∼ 


 


 
 
 


−1 
7 
5 
0 
0 
25 
15 
−5 
0 
0 
10 
−20 
0 
0 
−1 
2 


 


 
 
 


∼ 


 


 
 
 


−1 
7 
5 
0 
0 
5 
3 
−1 
0 
0 
1 
−2 
0 
0 
−1 
2 


 


 
 
 ∼ 


 


 
 
 


−1 
7 
5 
0 
0 
5 
3 
−1 
0 
0 
1 
−2 
0 
0 
0 
0 


 


 
 
 


Entonces 
S + T = (−1, 7, 5, 0), (0, 5, 3, −1), (0, 0, 1, −2) 


y, adem´as, el sistema de vectores 


BS+T = {(−1, 7, 5, 0), (0, 5, 3, −1), (0, 0, 1, −2)} 


es un sistema libre, por tanto BS+T es un base de S + T y dim(S + T) = 3. 
Para obtener las ecuaciones de S + T procedemos igual que hemos hecho para 
calcular las ecuaciones de T: 


(x, y, z, t) ∈ S + T ⇐⇒ rang 


 


 
 
 


x 
y 
z 
t 
−1 
7 
5 
0 
0 
5 
3 
−1 
0 
0 
1 
−2 


 


 
 
 = 3 


⇐⇒ 



x 
y 
z 
t 
−1 
7 
5 
0 
0 
5 
3 
−1 
0 
0 
1 
−2 



= 0 


ya que BS+T es una base. 
Pero 


x 
y 
z 
t 
−1 
7 
5 
0 
0 
5 
3 
−1 
0 
0 
1 
−2 



= 1 · (−1)4+3 



x 
y 
t 
−1 
7 
0 
0 
5 
−1 


+ (−2) · (−1)4+4 



x 
y 
z 
−1 
7 
5 
0 
5 
3 



= −(−7x − 5t − y) − 2(21x − 5z − (25x − 3y)) 
= 7x + 5t + y − 42 + 10z + 50x − 6y 
= 15x − 5y + 10z + 5t = 5(3x − y + 2z + t) 


Entonces 
S + T = {(x, y, z, t) ∈ R4 : 3x − y + 2z + t = 0} 
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Por ´ultimo 


S ∩ T = {(x, y, z, t) ∈ R4 : x − 2y + 3z − t = 0 , 
− 5y + 7z − 4t = 0 , x = 0 , y − 2z − t = 0} 


Para calcular las ecuaciones independientes de S ∩ T procedemos igual que hemos 
hecho con las ecuaciones de T; escalonamos la matriz de coeﬁcientes del sistema 
que determina S ∩ T: 


 


 
 
 


1 
−2 
3 
−1 
0 
−5 
7 
−4 
1 
0 
0 
0 
0 
1 
−2 
−1 


 


 
 
 ∼ 


 


 
 
 


1 
0 
0 
0 
0 
1 
−2 
−1 
0 
−5 
7 
−4 
1 
−2 
3 
−1 


 


 
 
 ∼ 


 


 
 
 


1 
0 
0 
0 
0 
1 
−2 
−1 
0 
−5 
7 
−4 
0 
−2 
3 
−1 


 


 
 
 


∼ 


 


 
 
 


1 
0 
0 
0 
0 
1 
−2 
−1 
0 
0 
−3 
−9 
0 
0 
−1 
−3 


 


 
 
 ∼ 


 


 
 
 


1 
0 
0 
0 
0 
1 
−2 
−1 
0 
0 
−3 
−9 
0 
0 
0 
0 


 


 
 
 ∼ 


 


 
 
 


1 
0 
0 
0 
0 
1 
−2 
−1 
0 
0 
1 
3 
0 
0 
0 
0 


 


 
 
 


Entonces 


S ∩ T = {(x, y, z, t) ∈ R4 : x = 0 , y − 2z − t = 0 , z + 3t = 0} 


Para obtener una base de S ∩ T resolvemos el sistema: 
 
 


 


x 
= 0 
y − 2z − 
t = 0 
z + 3t = 0 


Para resolverlo despejamos y en la segunda ecuaci´on del sistema: 


y = 2z + t 


y despejamos z en la tercera ecuaci´on del sistema: 


z = −3t 


Sustituyendo z = −3t en la ecuaci´on y = 2z + t, obtenemos: 


y = −2t 


Utilizando ahora la primera ecuaci´on del sistema, x = 0, se tiene que las inﬁnitas 
soluciones del mismo son: (0, −2t, −3t), con t ∈ R. 
Por tanto 


S ∩ T = {(x, y, z, t) ∈ R4 : x = 0 , y = −5t , z = −3t} 
= {(0, −5t, −3t, t) : t ∈ R} = {t(0, −5, −3, 1) : t ∈ R} 
= (0, −5, −3, 1) 
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es decir 
S ∩ T = (0, −5, −3, 1) 


De aqu´ı 
BS∩T = {(0, −5, −3, 1)} 


es una base de S ∩ T y, as´ı, dim(S ∩ T) = 1. 
Evidentemente se cumple la tesis del Teorema 3.5.6. 


Ejercicio 3.5.10 
Se consideran los siguientes subespacios vectoriales de R4: 


S = (0, 1, 1, −1), (0, 1, 0, 1), (0, 0, 1, −2) 
T = {(x, y, z, t) ∈ R4 : x − 2y + 3z − t = 0 , 2x + y − z + 2t = 0} 


Calcula: 


(a) Una base y las ecuaciones de S + T. 


(b) Una base y las ecuaciones de S ∩ T. 



Para terminar este apartado, vamos a caracterizar los subespacios vectoriales de 
R2 y R3. 
Subespacios vectoriales de R2 


Subespacios impropios 


• {(0, 0)}, de dimensi´on 0. 


• R2, de dimensi´on 2. 


Subespacios propios 


• Los ´unicos subespacios propios de R2 son los de dimensi´on 1, que son 
las rectas que pasan por el origen. 


Para demostrar esta aﬁrmaci´on consideramos un subespacio vectorial de 
R2 de dimensi´on 1. ´Este deber´a estar generado por un vector no nulo: 


S = (u, v) 


con (u, v) = (0, 0). 
Entonces 


(x, y) ∈ S ⇐⇒ rang 
x 
y 
u 
v 


= 1 ⇐⇒ v x − u y = 0 


Es decir 
S = {(x, y) ∈ R2 : v x − u y = 0} 


es una recta que pasa por el origen. 
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Subespacios vectoriales de R3 


Subespacios impropios 


• {(0, 0, 0)}, de dimensi´on 0. 


• R3, de dimensi´on 3. 


Subespacios propios 


• R3 tiene subespacios propios de dimensi´on 1 (rectas que pasan por el 
origen) y de dimensi´on 2 (planos que pasan por el origen). 


Para demostrar esta aﬁrmaci´on, en primer lugar, consideramos un subes- 
pacio vectorial de R3 de dimensi´on 1. ´Este deber´a estar generado por un 
vector no nulo: 
S = (u, v, w) 


con (u, v, w) = (0, 0, 0). Supongamos que, por ejemplo, u = 0. 
Entonces 


(x, y, z) ∈ S ⇐⇒ rang 
x 
y 
z 
u 
v 
w 


= 1 


⇐⇒ 
x 
y 
u 
v 


= 0 , 
x 
z 
u 
w 


= 0 


⇐⇒ v x − u y = 0 , w x − u z = 0 


es decir 


S = {(x, y, z) ∈ R3 : v x − u y = 0 , w x − u z = 0} 


es una recta (intersecci´on de dos planos) que pasa por el origen. 
En segundo lugar, consideramos un subespacio vectorial de R3 de di- 
mensi´on 2. ´Este deber´a estar generado por dos vectores linealmente in- 
dependientes: 
S = (u1, u2, u3), (v1, v2, v3) 


tal que 


rang 
u1 
u2 
u3 
v1 
v2 
v3 


= 2 


Entonces 


(x, y, z) ∈ S ⇐⇒ rang 


 


 
x 
y 
z 
u1 
u2 
u3 
v1 
v2 
v3 


 


 = 2 ⇐⇒ 



x 
y 
z 
u1 
u2 
u3 
v1 
v2 
v3 


= 0 


Pero 


x 
y 
z 
u1 
u2 
u3 
v1 
v2 
v3 


= A x + B y + C z 
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siendo 


A = 
u2 
u3 
v2 
v3 


, 
B = − 
u1 
u3 
v1 
v3 


, 
C = 
u1 
u2 
v1 
v2 



Entonces 


S = {(x, y, z) ∈ R3 : A x + B y + C z = 0} 


es un plano que pasa por el origen. 


3.5.2. 
Cambios de base 


Sean E un K-espacio vectorial, B = {v1,v2, . . . ,vn} una base de E y v ∈ E. 
Entonces 
v = λ1 · v1 + λ2 · v2 + · · · + λn · vn 
donde (λ1, λ2, . . . , λn) son las coordenadas del vector v en la base B. 
Sea ahora B′ = {v ′ 
1,v ′ 
2, . . . ,v ′ 
n} otra base de E. Entonces, los vectores de la 
base B tendr´an unas coordenadas en la base B′: 


vi = a1i · v ′ 
1 + a2i · v ′ 
2 + · · · + ani · v ′ 
n , 
i = 1, 2, . . . , n 


es decir, (a1i, a2i, . . . , ani) son las coordenadas de vi en la base B′. 
Si (λ′ 
1, λ′ 
2, . . . , λ′ 
n) son las coordenadas del vector v en la base B′, entonces 


v = λ′ 
1 · v ′ 
1 + λ′ 
2 · v ′ 
2 + · · · + λ′ 
n · v ′ 
n 


Por tanto, como 


λ′ 
1 · v ′ 
1 + λ′ 
2 · v ′ 
2 + · · · + λ′ 
n · v ′ 
n = v = λ1 · v1 + λ2 · v2 + · · · + λn · vn 


obtenemos que 


λ′ 
1 · v ′ 
1 + λ′ 
2 · v ′ 
2 + · · · + λ′ 
n · v ′ 
n = 
λ1 · (a11 · v ′ 
1 + a21 · v ′ 
2 + · · · + an1 · v ′ 
n)+ 
λ2 · (a12 · v ′ 
1 + a22 · v ′ 
2 + · · · + an2 · v ′ 
n) + · · · + 
λn · (a1n · v ′ 
1 + a2n · v ′ 
2 + · · · + ann · v ′ 
n) 


y as´ı 
λ′ 
1 − (λ1 · a11 + λ2 · a12 + · · · + λn · a1n) 
· v ′ 
1+ 
λ′ 
2 − (λ1 · a21 + λ2 · a22 + · · · + λn · a2n) 
· v ′ 
2 + · · · + 
λ′ 
n − (λ1 · an1 + λ2 · an2 + · · · + λn · ann) 
· v ′ 
n = 0 


De aqu´ı, como B′ es un sistema libre (sus vectores son linealmente independientes) 
por ser B′ una base de E, se tiene 


λ′ 
1 = λ1 · a11 + λ2 · a12 + · · · + λn · a1n 
λ′ 
2 = λ1 · a21 + λ2 · a22 + · · · + λn · a2n 
· · · · · · · · · · · · · · · · · · ·· · · · · · · · · · · · · · · 
λ′ 
n = λ1 · an1 + λ2 · an2 + · · · + λn · ann 
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es decir 
 


 
 
 
 


λ′ 
1 
λ′ 
2... 
λ′ 
n 


 


 
 
 
 = 


 


 
 
 
 


a11 
a12 
· · · 
a1n 
a21 
a22 
· · · 
a2n 
... 
... 
· · · 
... 
an1 
an2 
· · · 
ann 


 


 
 
 
 · 


 


 
 
 
 


λ1 
λ2 
... 
λn 


 


 
 
 
 


La matriz 


MB′B = 


 


 
 
 
 


a11 
a12 
· · · 
a1n 
a21 
a22 
· · · 
a2n 
... 
... 
· · · 
... 
an1 
an2 
· · · 
ann 


 


 
 
 
 


se llama matriz cambio de base de la base B a la base B′. 
Si llamamos 


XB = 


 


 
 
 
 


λ1 
λ2 
... 
λn 


 


 
 
 
 , 
XB′ = 


 


 
 
 
 


λ′ 
1 
λ′ 
2... 
λ′ 
n 


 


 
 
 
 


tenemos que 
XB′ = MB′B · XB 


Las ecuaciones anteriores referidas se llaman ecuaciones de cambio de base de la 
base B a la base B′. 
An´alogamente se pueden obtener las ecuaciones de cambio de base de la base 
B′ a la base B: 
XB = MBB′ · X′ 
B 
donde MBB′ es la matriz cambio de base de la base B′ a la base B. Entonces 


XB′ = MB′B · XB = MB′B · MBB′ · X′ 
B 


Como el vector v es arbitrario, es f´acil comprobar que 


MB′B · MBB′ = I 


Igualmente, tambi´en se puede comprobar que 


MBB′ · MB′B = I 


Por tanto, la matriz MB′B es una matriz regular y 


(MB′B)−1 = MBB′ 


Ejercicio 3.5.11 
Demuestra que: 


(a) MB′B · MBB′ = I. 


(b) MBB′ · MB′B = I. 
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Ejemplo 3.5.7 
Sabiendo que las coordenadas del vector v ∈ R3 en la base 


B = {(1, 0, −1), (2, 1, 3), (3, 1, −2)} 


son (2, 1, −4), vamos a calcular las coordenadas del vector v en la base 


B′ = {(7, 2, −1), (−4, −1, 7), (3, 1, 2)} 


utilizando la matriz cambio de base. 


Para calcular la matriz cambio de base de la base B a la base B′, tenemos que 
obtener las coordenadas de los vectores de la base B respecto a la base B′: 


(1, 0, −1) = a11 · (7, 2, −1) + a21 · (−4, −1, 7) + a31 · (3, 1, 2) 
(2, 1, 3) = a12 · (7, 2, −1) + a22 · (−4, −1, 7) + a32 · (3, 1, 2) 
(3, 1, −2) = a13 · (7, 2, −1) + a23 · (−4, −1, 7) + a33 · (3, 1, 2) 


Este sistema de ecuaciones lineales cuyos coeﬁcientes y los t´erminos independien- 
tes son vectores, da lugar a tres sistemas de ecuaciones lineales, de forma que el 
determinante de la matriz de coeﬁcientes (de los tres sistemas) es: 


7 
−4 
3 
2 
−1 
1 
−1 
7 
2 


= −14 + 4 + 42 − (3 + 49 − 16) = −4 = 0 


Por tanto, estos sistemas son de Cramer. 
De la primera ecuaci´on vectorial obtenemos el sistema 
 
 


 


7a11 − 4a21 + 3a31 = 
1 
2a11 − 
a21 + 
a31 = 
0 
−a11 + 7a21 + 2a31 = −1 


cuya soluci´on es 


a11 = 



1 
−4 
3 
0 
−1 
1 
−1 
7 
2 




7 
−4 
3 
2 
−1 
1 
−1 
7 
2 



= −2 + 4 − (3 + 7) 


−4 
= −8 


−4 = 2 


a21 = 



7 
1 
3 
1 
0 
1 
−1 
−1 
2 




7 
−4 
3 
2 
−1 
1 
−1 
7 
2 



= −1 − 6 − (−7 + 4) 


−4 
= −4 


−4 = 2 
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a31 = 



7 
−4 
1 
1 
−1 
0 
−1 
7 
−1 




7 
−4 
3 
2 
−1 
1 
−1 
7 
2 



= 7 + 14 − (1 + 8) 


−4 
= 12 


−4 = −3 


De la segunda ecuaci´on vectorial obtenemos el sistema 
 
 


 


7a12 − 4a22 + 3a32 = 2 
2a12 − 
a22 + 
a32 = 1 
−a12 + 7a22 + 2a32 = 3 


cuya soluci´on es 


a12 = 



2 
−4 
3 
1 
−1 
1 
3 
7 
2 




7 
−4 
3 
2 
−1 
1 
−1 
7 
2 



= −4 − 12 + 21 − (−9 + 14 − 8) 


−4 
= 8 


−4 = −2 


a22 = 



7 
2 
3 
1 
1 
1 
−1 
3 
2 




7 
−4 
3 
2 
−1 
1 
−1 
7 
2 



= 14 − 2 + 18 − (−3 + 21 + 8) 


−4 
= 4 


−4 = −1 


a32 = 



7 
−4 
2 
1 
−1 
1 
−1 
7 
3 




7 
−4 
3 
2 
−1 
1 
−1 
7 
2 



= −21 + 4 + 28 − (2 + 49 − 24) 


−4 
= −16 


−4 = 4 


De la tercera ecuaci´on vectorial obtenemos el sistema 
 
 


 


7a13 − 4a23 + 3a33 = 
3 
2a13 − 
a23 + 
a33 = 
1 
−a13 + 7a23 + 2a33 = −2 


cuya soluci´on es 


a13 = 



3 
−4 
3 
1 
−1 
1 
−2 
7 
2 




7 
−4 
3 
2 
−1 
1 
−1 
7 
2 



= −6 + 8 + 21 − (6 + 21 − 8) 


−4 
= 4 


−4 = −1 
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a23 = 



7 
3 
3 
1 
1 
1 
−1 
−2 
2 




7 
−4 
3 
2 
−1 
1 
−1 
7 
2 



= 14 − 3 − 12 − (−3 − 14 + 12) 


−4 
= 4 


−4 = −1 


a33 = 



7 
−4 
3 
1 
−1 
1 
−1 
7 
−2 




7 
−4 
3 
2 
−1 
1 
−1 
7 
2 



= 14 + 4 + 42 − (3 + 49 + 16) 


−4 
= −8 


−4 = 2 


Entonces, la matriz cambio de base de la base B a la base B′ es 


MB′B = 


 


 
2 
−2 
−1 
1 
−1 
−1 
−3 
4 
2 


 


 


De aqu´ı, las coordenadas del vector v en la base B′ son 
 


 
2 
−2 
−1 
1 
−1 
−1 
−3 
4 
2 


 


 · 


 


 
2 
1 
−4 


 


 = 


 


 
6 
5 
−10 


 


 


es decir, (6, 5, −10). 
El vector v es 


v = 2 · (1, 0, −1) + 1 · (2, 1, 3) + (−4) · (3, 1, −2) 
= 6 · (7, 2, −1) + 5 · (−4, −1, 7) + (−10) · (3, 1, 2) 
= (−8, −3, 9) 



Ejercicio 3.5.12 
Se consideran las siguientes bases B y B′ de R3: 


B = {(0, 1, 1), (1, 0, 1), (0, 1, 0)} , 
B′ = {(1, 1, 1), (1, 1, 0), (1, 0, 0)} 


Sea v un vector de R3 cuyas coordenadas en la base B son (1, 2, 3). Calcula las 
coordenadas del vector v en la base B′ aplicando la matriz cambio de base. 
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Tema 4 


Diagonalizaci´on de matrices 


4.1. 
Introducci´on 


El objetivo de este tema es, dada una matriz cuadrada A ∈ Mn(R), ver si existe 
una matriz regular P ∈ Mn(R) tal que 


D = P −1 · A · P 


sea una matriz diagonal. 
Todas las deﬁniciones y propiedades enunciadas en este tema se pueden gene- 
ralizar sustituyendo el cuerpo de los n´umeros reales por un cuerpo arbitrario K, a 
excepci´on del apartado dedicado a la diagonalizaci´on de matrices reales sim´etricas. 


4.2. 
Valores y vectores propios 


En este primer apartado vamos a introducir los conceptos de valor propio y 
vector propio de una matriz cuadrada. Como veremos en un ejemplo, no todas las 
matrices tienen valores propios; con la deﬁnici´on de la ecuaci´on caracter´ıstica en la 
siguiente secci´on, veremos que las ra´ıces de esta ecuaci´on son los valores propios 
de la matriz, por lo que, por ejemplo, demostrar que una matriz no tiene valores 
propios ser´a mucho m´as sencillo. 


Deﬁnici´on 4.2.1 (Valor propio) 
Dada A ∈ Mn(R), decimos que λ ∈ R es un valor propio o autovalor de A si 
existe, al menos, un vector no nulo v = (v1, v2, . . . , vn) ∈ Rn tal que 


A · 


 


 
 
 
 


v1 
v2 
... 
vn 


 


 
 
 
 = λ 


 


 
 
 
 


v1 
v2 
... 
vn 


 


 
 
 
 



Fundamentos Matem´aticos de la Ingenier´ıa – 503 – 1r. semestre 
2009/2010 
101 
c 
UJI 


09 
M. Barreda Rochera / J. A. López Ortí - ISBN: 978-84-692-9833-6 
Fundamentos Matemáticos de la Ingeniería. Parte I: Álgebra Lineal - UJI 


Ejemplo 4.2.1 
Veamos que la matriz 


A = 
0 
1 
−1 
0 



no tiene valores propios reales. 


Supongamos que λ ∈ R sea un valor propio de la matriz A. Entonces, existe un 
vector no nulo (x, y) ∈ R2 tal que 
0 
1 
−1 
0 


· 
x 
y 


= λ 
x 
y 



De aqu´ı, el vector (x, y) ∈ R2 tiene que ser soluci´on del sistema 
y = λ x 
−x = λ y 


Como (x, y) = (0, 0), se tiene que x = 0 ´o y = 0. Supongamos que x = 0. 
Sustituyendo la primera ecuaci´on en la segunda, obtenemos 


−x = λ2 x 


Por tanto 
x(1 + λ2) = 0 ⇐⇒ 1 + λ2 = 0 


La ecuaci´on 1 + λ2 = 0 no tiene ra´ıces reales, por lo tanto la matriz A no tiene 
valores propios reales. 


Deﬁnici´on 4.2.2 (Vector propio) 
Sea A ∈ Mn(R) y λ un valor propio de A. 
Decimos que un vector no nulo v = (v1, v2, . . . , vn) ∈ Rn es un vector propio o 
autovector de A correspondiente al valor propio λ si 


A · 


 


 
 
 
 


v1 
v2 
... 
vn 


 


 
 
 
 = λ 


 


 
 
 
 


v1 
v2 
... 
vn 


 


 
 
 
 



Deﬁnici´on 4.2.3 (Espectro) 
Dada A ∈ Mn(R), llamamos espectro de A, y lo representamos por σ(A), al 
conjunto de todos los valores propios de A; es decir: 


σ(A) = {λ ∈ R : λ es un valor propio de A} 
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Ejemplo 4.2.2 
El espectro de la matriz A del Ejemplo 4.2.1 es el conjunto vac´ıo, es decir 


σ(A) = ∅ 



Teorema 4.2.1 (Subespacio propio) 
Sean A ∈ Mn(R) y λ ∈ σ(A). El conjunto de los vectores propios de A corres- 
pondientes al valor propio λ constituyen, junto con el vector nulo, un subespacio 
vectorial de Rn, llamado subespacio propio de A correspondiente al valor pro- 
pio λ. 
El subespacio propio de A correspondiente al valor propio λ se denota por Eλ, es 
decir 


Eλ = {(v1, v2, . . . , vn) ∈ Rn : A · 


 


 
 
 
 


v1 
v2 
... 
vn 


 


 
 
 
 = λ 


 


 
 
 
 


v1 
v2 
... 
vn 


 


 
 
 
 


 
 
 
 
 


 
 
 
 



La demostraci´on de este teorema es inmediata ya que 


Eλ = {(v1, v2, . . . , vn) ∈ Rn : (A − λI) · 


 


 
 
 
 


v1 
v2 
... 
vn 


 


 
 
 
 = 


 


 
 
 
 


0 
0 
... 
0 


 


 
 
 
 


 
 
 
 
 


 
 
 
 


y, por tanto, Eλ es un subespacio vectorial de Rn cuya dimensi´on es: 


dim (Eλ) = n − rang(A − λI) 


Teorema 4.2.2 
Sean A ∈ Mn(R). Entonces: 


1. Dos valores propios distintos de A no tienen ning´un vector propio com´un, es 
decir 


λ, µ ∈ σ(A) , λ = µ 
=⇒ 
Eλ ∩ Eµ = {(0, 0, . . . , 0)} 


2. Si λ1, λ2, . . . , λr son un conjunto de r valores propios distintos de la matriz 
A y vi ∈ Eλi − {(0, 0, . . . , 0)}, entonces {v1,v2, . . . ,vr} es un sistema libre. 



Ejercicio 4.2.1 
Demuestra el primer apartado del Teorema 4.2.2. 
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4.3. 
Ecuaci´on caracter´ıstica 


En esta secci´on introducimos, entre otros, los conceptos de ecuaci´on caracter´ısti- 
ca y orden de un valor propio, estableciendo una relaci´on importante entre la dimen- 
si´on de un subespacio propio asociado a un valor propio y el orden del valor propio. 


Dada A ∈ Mn(R), un escalar λ ∈ R es un valor propio de A si existe un vector 
no nulo v = (v1, v2, . . . , vn) ∈ Rn tal que 


A · 


 


 
 
 
 


v1 
v2 
... 
vn 


 


 
 
 
 = λ 


 


 
 
 
 


v1 
v2 
... 
vn 


 


 
 
 
 


es decir 


(A − λI) · 


 


 
 
 
 


v1 
v2 
... 
vn 


 


 
 
 
 = 


 


 
 
 
 


0 
0 
... 
0 


 


 
 
 
 


Por tanto, para que un vector no nulo (v1, v2, . . . , vn) ∈ Rn sea soluci´on del sistema 
anterior, ese sistema debe ser compatible indeterminado (es un sistema homog´eneo). 
Entonces, aplicando el teorema de Rouch´e-Fr¨obenius, obtenemos 


λ es un valor propio de A 
⇐⇒ 
|A − λI| = 0 


Esto da lugar a la siguiente 


Deﬁnici´on 4.3.1 (Ecuaci´on caracter´ıstica) 
Sea A ∈ Mn(R). La ecuaci´on |A −λI| = 0 cuyas ra´ıces son los valores propios de 
A, se llama ecuaci´on caracter´ıstica de la matriz A. 


Ejemplo 4.3.1 
Veamos que aplicando la deﬁnici´on anterior es mucho m´as sencillo demostrar que 
la matriz del Ejemplo 4.2.1 no tiene valores propios reales. 


Calculamos la matriz A − λI: 


A − λI = 
0 
1 
−1 
0 


− λ 
1 
0 
0 
1 


= 
−λ 
1 
−1 
−λ 



Por tanto, la ecuaci´on caracter´ıstica de la matriz A es 
−λ 
1 
−1 
−λ 


= 0 


es decir 
1 + λ2 = 0 


ecuaci´on que no tiene ra´ıces reales, por tanto la matriz A no tiene valores propios 
reales. 
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Ejemplo 4.3.2 
Vamos a calcular los valores propios de la matriz 


A = 


 


 
1 
3 
2 
1 
−1 
−2 
−1 
3 
4 


 


 


La ecuaci´on caracter´ıstica de la matriz A es 


1 − λ 
3 
2 
1 
−1 − λ 
−2 
−1 
3 
4 − λ 


= 0 


es decir 


(1 − λ)(−1 − λ)(4 − λ) + 6 + 6 − 
2(1 + λ) − 6(1 − λ) + 3(4 − λ) 
= 0 ; 


(−1 + λ2)(4 − λ) + 12 − 2 − 2λ + 6 − 6λ − 12 + 3λ = 0 ; 


−4 + 4λ2 + λ − λ3 + 4 − 5λ = 0 ; 


−λ3 + 4λ2 − 4λ = 0 ; 


Resolvemos ahora la ecuaci´on caracter´ıstica de la matriz A: 


−λ(λ2 − 4λ + 4) = 0 


λ = 0 


λ2 − 4λ + 4 = 0 ; λ = 4 ± √16 − 16 


2 
= 2 


Por tanto 
σ(A) = {0, 2} 


Hay que se˜nalar aqu´ı que la ra´ız λ = 2 de la ecuaci´on caracter´ıstica de la matriz A 
es una ra´ız doble. 


Ejercicio 4.3.1 
Calcula los valores propios de la matriz 


A = 


 


 
2 
0 
0 
5 
2 
5 
−5 
0 
−3 


 


 



Deﬁnici´on 4.3.2 (Polinomio caracter´ıstico) 
Sea A ∈ Mn(R). El polinomio p(λ) = |A−λI| se llama polinomio caracter´ıstico 
de la matriz A. 
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Ejemplo 4.3.3 
El polinomio caracter´ıstico de la matriz del Ejemplo 4.3.2 es 


p(λ) = −λ3 + 4λ2 − 4λ 



Ejercicio 4.3.2 
Calcula el polinomio caracter´ıstico de la matriz del Ejercicio 4.3.1. 


Deﬁnici´on 4.3.3 (Orden de un valor propio) 
Sea A ∈ Mn(R). Si λ es una ra´ız de la ecuaci´on caracter´ıstica de A, con orden de 
multiplicidad α, se dice que λ es un valor propio de orden α de la matriz A. 


Ejemplo 4.3.4 
Los valores propios de la matriz del Ejemplo 4.3.2 son λ1 = 0 de orden α1 = 1 y 
λ2 = 2 de orden α2 = 2. 


Ejercicio 4.3.3 
Calcula el orden de los valores propios de la matriz del Ejercicio 4.3.1. 


Teorema 4.3.1 (Relaci´on entre dimensi´on y orden) 
Sean A ∈ Mn(R) y λ un valor propio de A de orden α. Entonces 


1 dim (Eλ) α 



Ejemplo 4.3.5 
Vamos a calcular los subespacios propios asociados a los valores propios de la ma- 
triz del Ejemplo 4.3.2. 


En primer lugar, calculamos el subespacio propio asociado al valor propio λ1 = 0: 


Eλ1 = {(x, y, z) ∈ R3 : (A − λ1I) · 


 


 
x 
y 
z 


 


 = 


 


 
0 
0 
0 


 


 


 
 


 


= {(x, y, z) ∈ R3 : 


 


 
1 
3 
2 
1 
−1 
−2 
−1 
3 
4 


 


 · 


 


 
x 
y 
z 


 


 = 


 


 
0 
0 
0 


 


 


 
 


 


Escalonando la matriz de coeﬁcientes del sistema que deﬁne el subespacio propio 
asociado al valor propio λ1 = 0 
 


 
1 
3 
2 
1 
−1 
−2 
−1 
3 
4 


 


 ∼ 


 


 
1 
3 
2 
0 
4 
4 
0 
6 
6 


 


 ∼ 


 


 
1 
3 
2 
0 
4 
4 
0 
0 
0 


 


 ∼ 


 


 
1 
3 
2 
0 
1 
1 
0 
0 
0 


 


 
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obtenemos 


Eλ1 = {(x, y, z) ∈ R3 : 


 


 
1 
3 
2 
0 
1 
1 
0 
0 
0 


 


 · 


 


 
x 
y 
z 


 


 = 


 


 
0 
0 
0 


 


 


 
 


 


= {(x, y, z) ∈ R3 : x + 3y + 2z = 0 , y + z = 0} 


Resolvemos ahora el sistema de ecuaciones lineal: 
x + 3y + 2z = 0 
y + z = 0 


Para resolverlo despejamos x en la primera ecuaci´on del sistema: 


x = −3y − 2z 


y despejamos y en la segunda ecuaci´on del sistema: 


y = −z 


Sustituyendo y = −z en la ecuaci´on x = −3y − 2z, obtenemos: 


x = z 


Por tanto, las inﬁnitas soluciones del sistema son: (z, −z, z), con z ∈ R. De aqu´ı: 


Eλ1 = {(x, y, z) ∈ R3 : x = z , y = −z} = {(z, −z, z) : z ∈ R} 
= {z(1, −1, 1) : z ∈ R} = (1, −1, 1) 


Entonces 
Bλ1 = {(1, −1, 1)} 


es una base del subespacio propio Eλ1, por tanto 


d1 = dim (Eλ1) = 1 = α1 


En segundo lugar, calculamos el subespacio propio asociado al valor propio λ2 = 2: 


Eλ2 = {(x, y, z) ∈ R3 : (A − λ2I) · 


 


 
x 
y 
z 


 


 = 


 


 
0 
0 
0 


 


 


 
 


 


= {(x, y, z) ∈ R3 : 


 


 
−1 
3 
2 
1 
−3 
−2 
−1 
3 
2 


 


 · 


 


 
x 
y 
z 


 


 = 


 


 
0 
0 
0 


 


 


 
 


 


Escalonando la matriz de coeﬁcientes del sistema que deﬁne el subespacio propio 
asociado al valor propio λ2 = 2 
 


 
−1 
3 
2 
1 
−3 
−2 
−1 
3 
2 


 


 ∼ 


 


 
−1 
3 
2 
0 
0 
0 
0 
0 
0 


 


 


Fundamentos Matem´aticos de la Ingenier´ıa – 503 – 1r. semestre 
2009/2010 
107 
c 
UJI 


5 
M. Barreda Rochera / J. A. López Ortí - ISBN: 978-84-692-9833-6 
Fundamentos Matemáticos de la Ingeniería. Parte I: Álgebra Lineal - UJI 


obtenemos 


Eλ2 = {(x, y, z) ∈ R3 : 


 


 
−1 
3 
2 
0 
0 
0 
0 
0 
0 


 


 · 


 


 
x 
y 
z 


 


 = 


 


 
0 
0 
0 


 


 


 
 


 


= {(x, y, z) ∈ R3 : −x + 3y + 2z = 0} 
= {(x, y, z) ∈ R3 : x = 3y + 2z} 
= {(3y + 2z, y, z) : y, z ∈ R} 
= {y(3, 1, 0) + z(2, 0, 1) : y, z ∈ R} 
= (3, 1, 0), (2, 0, 1) 


Entonces, como el sistema de vectores 


Bλ2 = {(3, 1, 0), (2, 0, 1)} 


es libre, obtenemos que Bλ2 es una base del subespacio propio Eλ2 y, as´ı 


d2 = dim (Eλ2) = 2 = α2 



Ejercicio 4.3.4 
Calcula los subespacios propios de la matriz del Ejercicio 4.3.1. 


4.4. 
Matrices diagonalizables 


En el Tema 1 introdujimos las matrices diagonales. Tales matrices tienen, entre 
otras, una propiedad importante: si D ∈ Mn(R) es una matriz diagonal 


D = diag [λ1, λ2, . . . , λn] 


entonces 
σ(D) = {λ1, λ2, . . . , λn} 


es decir, la matriz D tiene exactamente n valores propios, aunque puede que no sean 
todos distintos. 
Por otra parte, si k ∈ N, se tiene que 


Dk = 


k veces 
D · D · · · · · D 


es tambi´en una matriz diagonal y 


Dk = diag 
(λ1)k, (λ2)k, . . . , (λn)k 


En esta secci´on vamos a introducir un tipo de matrices, las matrices diagonali- 
zables, cuya potencia es tambi´en f´acil de calcular. Como aplicaci´on obtendremos el 
t´ermino general de la sucesi´on de Fibonacci. 
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Deﬁnici´on 4.4.1 (Matriz diagonalizable) 
Una matriz A ∈ Mn(R) se dice que es diagonalizable, si existe una matriz regular 
P ∈ Mn(R) tal que D = P −1 · A · P es una matriz diagonal. La matriz D se llama 
forma diagonal y la matriz regular P se llama matriz de paso. 


Evidentemente, toda matriz diagonal es diagonalizable. 


Deﬁnici´on 4.4.2 (Matrices semejantes) 
Dos matrices A, B ∈ Mn(R) se dice que son semejantes si existe una matriz 
regular P ∈ Mn(R) tal que B = P −1 · A · P. 


Por tanto, una matriz es diagonalizable si es semejante a una matriz diagonal. 


Veamos ahora que dos matrices semejantes tienen los mismos valores propios. 
Sean A, B ∈ Mn(R) dos matrices semejantes. Entonces, existe una matriz regular 
P ∈ Mn(R) tal que 
B = P −1 · A · P 


Calculamos ahora el polinomio caracter´ıstico de la matriz B: 


|B − λI| = |P −1 · A · P − λI| = |P −1 · A · P − λP −1 · I · P| 
= |P −1 · (A − λI) · P| = |P −1| |A − λI| |P| = |A − λI| 


Por tanto, el polinomio caracter´ıstico de la matriz B es igual al polinomio carac- 
ter´ıstico de la matriz A; es decir, las matrices A y B tienen los mismos valores 
propios. 
Tal y como hemos se˜nalado anteriormente, una matriz diagonal de tama˜no n 
tiene exactamente n valores propios (no necesariamente todos distintos). Entonces, 
para que una matriz A ∈ Mn(R) sea diagonalizable, al ser ´esta semejante a una 
matriz diagonal, deber´a tener n valores propios (no necesariamente todos distintos); 
adem´as, su forma diagonal estar´a formada por los valores propios de A (cada uno 
tantas veces como indique su orden). 
Este resultado se expresa m´as formalmente de la siguiente forma. Supongamos 
que A ∈ Mn(R) es una matriz diagonalizable tal que λ1, λ2, . . . , λr son sus valores 
propios con ´ordenes respectivos α1, α2, . . . , αr. Entonces 


α1 + α2 + · · · + αr = n 


Pero esta condici´on no es suﬁciente; para que lo sea, adem´as, las dimensiones de 
los subespacios propios tienen que ser iguales a los ´ordenes de los valores propios 
respectivos. 


Teorema 4.4.1 (Caracterizaci´on de matrices diagonalizables) 
Sea A ∈ Mn(R) tal que λ1, λ2, . . . , λr son sus valores propios con ´ordenes respec- 
tivos α1, α2, . . . , αr y sea di = dim (Eλi), para i = 1, 2, . . . , r. Entonces 


A es diagonalizable ⇐⇒ 
α1 + α2 + · · · + αr = n 
αi = di , i = 1, 2, . . . , r 


⇐⇒ d1 + d2 + · · · + dr = n 
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La forma diagonal de la matriz A es 


D = diag 


 


λ1, . . . , λ1 
α1 


, λ2, . . . , λ2 
α2 


, . . . , λr, . . . , λr 
αr 


 


 


y si 


Eλ1 = v1,v2, . . . ,vα1 
Eλ2 = vα1+1,vα1+2, . . . ,vα1+α2 
· · · · · · · · ·· · · · · · · · · · · · · · · · · ·· · · 


Eλr = 
vα1+α2+···+αr−1+1,vα1+α2+···+αr−1+2, . . . ,vn 


entonces la matriz de paso es 


P = 
V1, V2, . . . , Vα1, . . . , Vα1+···+αr−1+1, . . . , Vn 


donde Vi es la matriz columna cuyos elementos son las componentes del vector vi, 
para i = 1, 2, . . . , n. 


Ejemplo 4.4.1 
La matriz del Ejemplo 4.3.2 es diagonalizable, ya que como hemos visto en el Ejem- 
plo 4.3.5 
d1 = 1 = α1 , d2 = 2 = α2 


o lo que es lo mismo 


d1 + d2 = 1 + 2 = 3 = tama˜no de la matriz A 


Adem´as, la forma diagonal es 


D = 


 


 
0 
0 
0 
0 
2 
0 
0 
0 
2 


 


 


y la matriz de paso es 


P = 


 


 
1 
3 
2 
−1 
1 
0 
1 
0 
1 


 


 



Ejercicio 4.4.1 
Demuestra que la matriz 


A = 


 


 
4 
0 
−20 
2 
0 
−10 
1 
−1 
−2 


 


 


no es diagonalizable. 
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Ejercicio 4.4.2 
Determina si la matriz 


A = 


 


 
4 
−2 
−2 
3 
−3 
−6 
2 
−4 
−1 


 


 


es diagonalizable y, en caso aﬁrmativo, calcula su forma diagonal y la matriz de 
paso. 


Como consecuencia del teorema de caracterizaci´on de matrices diagonaliza- 
bles, Teorema 4.4.1, si una matriz A ∈ Mn(R) tiene n valores propios distintos 
λ1, λ2, . . . , λn, es decir, sus ´ordenes respectivos son α1 = 1, α2 = 1, . . . , αn = 1, 
como 
1 di αi , i = 1, 2, . . . n 


se tiene que 
di = 1 , i = 1, 2, . . . n 


y entonces 
d1 + d2 + · · · + dn = n 


Por tanto, la matriz A es diagonalizable. 


Ejercicio 4.4.3 
Determina, en cada caso, si la matriz es diagonalizable y en caso aﬁrmativo, calcula 
su forma diagonal, as´ı como la matriz de paso: 


(a) A = 
1 
0 
1 
2 


(b) A = 
7 
5 
−10 
−8 



(c) A = 
−1 
−2 
4 
5 


(d) A = 
−7 
−5 
16 
17 



(e) A = 
0 
−1 
1 
0 


(f ) A = 
2 
4 
0 
2 



(g) A = 


 


 
−1 
0 
1 
−7 
2 
5 
3 
0 
1 


 


 
(h) A = 


 


 
7 
−10 
0 
3 
−4 
0 
1 
−2 
2 


 


 


(i) A = 


 


 
−2 
0 
0 
−5 
−2 
−5 
5 
0 
3 


 


 
(j) A = 


 


 
2 
0 
3 
0 
2 
−1 
0 
0 
2 


 


 


(k) A = 


 


 
−4 
6 
−12 
3 
−1 
6 
3 
−3 
8 


 


 
(l) A = 


 


 
−3 
5 
−20 
2 
0 
8 
2 
1 
7 


 


 
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(m) A = 


 


 
−2 
0 
−1 
0 
2 
0 
3 
0 
2 


 


 
(n) A = 


 


 
1 
1 
3 
3 
−1 
3 
0 
0 
−2 


 


 



Una aplicaci´on de las matrices diagonalizables es el c´alculo de una potencia de 
una matriz (diagonalizable). Supongamos que A ∈ Mn(R) es una matriz diagona- 
lizable. Entonces, existe una matriz regular P ∈ Mn(R) tal que 


D = P −1 · A · P 


es una matriz diagonal. Multiplicando a la izquierda por la matriz P en ambas partes 
de la igualdad, obtenemos 


P · D = P · P −1 · A · P = I · A · P = A · P 


Multiplicando ahora a la derecha por la matriz P −1 en ambas partes de la igualdad, 
tenemos 
P · D · P −1 = A · P · P −1 = A · I = A 


es decir 
A = P · D · P −1 


De aqu´ı 


A2 = P · D · P −1 · P · D · P −1 = P · D · I · D · P −1 


= P · D · D · P −1 = P · D2 · P −1 


Igualmente se puede demostrar que 


A3 = P · D3 · P −1 


y en general 
An = P · Dn · P −1 


para cualquier n ∈ N. 


Ejercicio 4.4.4 
Calcula A30 en los siguientes casos: 


(a) A = 
−1 
−2 
4 
5 



(b) A = 


 


 
−3 
5 
−20 
2 
0 
8 
2 
1 
7 


 


 


(c) A = 


 


 
2 
0 
0 
−1 
1 
1 
1 
1 
1 


 


 
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En ocasiones debemos determinar si una familia de matrices que dependen de 
un par´ametro son o no diagonalizables en funci´on de dicho par´ametro. 


Ejemplo 4.4.2 
Dada la matriz 


A = 


 


 
5 
0 
0 
0 
−1 
0 
3 
0 
a 


 


 


vamos a estudiar si es diagonalizable en funci´on de los valores del par´ametro a ∈ R. 


En primer lugar calculamos los valores propios de la matriz A: 


5 − λ 
0 
0 
0 
−1 − λ 
0 
3 
0 
a − λ 


= 0 


es decir 
(5 − λ)(−1 − λ)(a − λ) = 0 


Por tanto, los valores propios de A son λ1 = 5, λ2 = −1 y λ3 = a. 
Esto da lugar a los tres siguientes casos: 


1 a = −1, 5. 


Los valores propios de la matriz A, en este caso, son 


λ1 = 5 , α1 = 1 ; 
1 d1 α1 = 1 
=⇒ 
d1 = 1 
λ2 = −1 , α2 = 1 ; 
1 d2 α2 = 1 
=⇒ 
d2 = 1 
λ3 = a , α3 = 1 ; 
1 d3 α3 = 1 
=⇒ 
d3 = 1 


De aqu´ı 
d1 + d2 + d3 = 3 = tama˜no de la matriz A 


y entonces, la matriz A es diagonalizable. 


2 a = −1. 


Los valores propios de la matriz A, en este caso, son 


λ1 = 5 , α1 = 1 ; 
1 d1 α1 = 1 
=⇒ 
d1 = 1 
λ2 = −1 , α2 = 2 ; 
1 d1 α2 = 2 
=⇒ 
d2 = 1 ´o d2 = 2 


Pero, de acuerdo con lo expuesto despu´es del Teorema 4.2.1: 


d2 = dim(Eλ2) = 3 − rang(A − λ2I) 


y como 


A − λ2I = 


 


 
6 
0 
0 
0 
0 
0 
3 
0 
0 


 


 ∼ 


 


 
6 
0 
0 
0 
0 
0 
0 
0 
0 


 


 


se tiene que 
rang(A − λ2I) = 1 
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y as´ı, d2 = 2. Por tanto 


d1 + d2 = 3 = tama˜no de la matriz A 


y entonces, la matriz A es diagonalizable. 


3 a = 5. 


Los valores propios de la matriz A, en este caso, son 


λ1 = 5 , α1 = 2 ; 
1 d1 α1 = 2 
=⇒ 
d1 = 1 ´o d1 = 2 
λ2 = −1 , α2 = 1 ; 
1 d1 α2 = 1 
=⇒ 
d2 = 1 


Pero, de acuerdo con lo expuesto despu´es del Teorema 4.2.1: 


d1 = dim(Eλ1) = 3 − rang(A − λ1I) 


y como 


A − λ1I = 


 


 
0 
0 
0 
0 
−6 
0 
3 
0 
0 


 


 ∼ 


 


 
1 
0 
0 
0 
1 
0 
0 
0 
0 


 


 


se tiene que 
rang(A − λ2I) = 2 


y as´ı, d1 = 1. Por tanto 


d1 + d2 = 2 < 3 = tama˜no de la matriz A 


y entonces, la matriz A no es diagonalizable. 



Ejercicio 4.4.5 
Estudia, en funci´on del par´ametro a, si la matriz 


A = 


 


 
5 
0 
0 
0 
−1 
a 
3 
0 
5 


 


 


es diagonalizable. 


La Sucesi´on de Fibonacci 


La sucesi´on num´erica 


1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, . . . 


se llama sucesi´on de Fibonacci. Los dos primeros t´erminos son 1 y 1, y el siguiente 
se obtiene al sumar los dos anteriores. As´ı, en general, si u0 = 1 y u1 = 1, entonces 
para n 2: 
un = un−1 + un−2 
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La sucesi´on de Fibonacci fue introducida por el matem´atico italiano Leonardo 
Pisano (1175-1226), conocido como Fibonacci, que a su vez es contracci´on de ﬁlius 
Bonaci (literalmente “hijo de Bonacci”). En su “libro del ´abaco” plantea el siguiente 
problema: 
Sea una pareja de conejos j´ovenes. Los conejos tardan un mes en madurar y a 
partir del mes siguiente, cada pareja madura de conejos da a luz a una pareja de 
conejos j´ovenes. ¿Cu´al es el n´umero de parejas de conejos que habr´a transcurridos 
n meses? 


La siguiente tabla muestra el n´umero de parejas de conejos que hay en cada uno 
de los seis primeros meses, tal y como se desprende de la ﬁgura anterior: 


Mes 
Primero 
Segundo 
Trecero 
Cuarto 
Quinto 
Sexto 


Parejas de conejos 
1 
1 
2 
3 
5 
8 


La sucesi´on de Fibonacci se puede escribir en forma matricial, como: 
un 
un−1 


= 
1 
1 
1 
0 


un−1 
un−2 


, 
n 2 


con u0 = 1 y u1 = 1. 
El problema consiste, pues, en obtener el t´ermino general de la sucesi´on de 
Fibonacci, un, ya que un corresponde con el n´umero de parejas de conejos que hay 
en el n-´esimo mes. 
Haciendo 


A = 
1 
1 
1 
0 
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w0 = 
u1 
u0 


= 
1 
1 


, w1 = 
u2 
u1 


= 
2 
1 


, . . . , 


wn−2 = 
un−1 
un−2 


, wn−1 = 
un 
un−1 



obtenemos que la sucesi´on de Fibonacci se puede escribir en la forma: 


wn−1 = A · wn−2 , 
n 2 


y as´ı 


w1 = A · w0 
w2 = A · w1 = A · (A · w0) = A2 · w0 
... 
wn−1 = A · wn−2 = · · · = An−1 · w0 


El problema consiste, pues, en calcular An−1, lo cual es bastante engorroso cuando n 
es grande. En cambio, si la matriz A es diagonalizable el problema es m´as sencillo. 
Veamos si ´este es el caso. Empezamos calculando los valores propios de A: 


|A − λI| = 0 ; 
1 − λ 
1 
1 
−λ 


= 0 ; 
(1 − λ)(−λ) − 1 = 0 ; 
λ2 − λ − 1 = 0; 


λ = 1 ± √1 + 4 


2 
= 1 ± 
√ 


5 


2 
es decir, los valores propios (simples) de A son: 


λ1 = 1 + 
√ 


5 


2 
, 
λ2 = 1 − 
√ 


5 


2 


Al ser la matriz de tama˜no 2 y tener dos valores propios distintos, se tiene que la 
matriz A es diagonalizable. La forma diagonal, es: 


D = 


 


 
 


1 + 
√ 


5 


2 
0 


0 
1 − 
√ 


5 


2 


 


 
 = 1 


2 


1 + 
√ 


5 
0 
0 
1 − 
√ 


5 



Calculamos ahora la matriz de paso. Para ello, obtenemos los subespacios pro- 
pios correspondientes a los valores propios λ1 y λ2. 
Comenzamos con Eλ1: 


Eλ1 = {(x, y) ∈ R2 : (A − λ1I) 
x 
y 


= 
0 
0 



Pero 


A − λ1I = 


 


 
 
1 − 1 + 
√ 


5 


2 
1 


1 
−1 + 
√ 


5 


2 


 


 
 = 


 


 
 


1 − 
√ 


5 


2 
1 


1 
−1 + 
√ 


5 


2 


 


 
 
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y as´ı, la expresi´on matricial del sistema 


(A − λ1I) 
x 
y 


= 
0 
0 



es 
 


 
 


1 − 
√ 


5 


2 
1 


1 
−1 + 
√ 


5 


2 


 


 
 
x 
y 


= 
0 
0 



es decir 
 
 
 


 
 


1 − 
√ 


5 


2 
x 
+ 
y 
= 
0 


x 
− 
1 + 
√ 


5 


2 
y 
= 
0 


Multiplicando la segunda ecuaci´on por (1 − 
√ 


5)/2, obtenemos 


1 − 
√ 


5 


2 
x − 1 − 
√ 


5 


2 
1 + 
√ 


5 


2 
y = 0 =⇒ 1 − 
√ 


5 


2 
x − 1 − 5 


4 
y = 0 


=⇒ 1 − 
√ 


5 


2 
x + y = 0 


es decir, la primera ecuaci´on. Por tanto, las dos ecuaciones no son independientes. 
Entonces 


Eλ1 = 
(x, y) ∈ R2 : y = 
√ 


5−1 
2 
x 
= 
(x, 
√ 


5−1 
2 
x) : x ∈ R 


= 
x (1, 
√ 


5−1 
2 
) : x ∈ R 
= 
(1, 
√ 


5−1 
2 
) 
= 
(2, 
√ 


5 − 1) 


Ahora hacemos lo mismo con Eλ2: 


Eλ2 = {(x, y) ∈ R2 : (A − λ2I) 
x 
y 


= 
0 
0 



Pero 


A − λ2I = 


 


 
 
1 − 1 − 
√ 


5 


2 
1 


1 
−1 − 
√ 


5 


2 


 


 
 = 


 


 
 


1 + 
√ 


5 


2 
1 


1 
√ 


5 − 1 
2 


 


 
 


y as´ı, la expresi´on matricial del sistema 


(A − λ2I) 
x 
y 


= 
0 
0 



es 
 


 
 


1 + 
√ 


5 


2 
1 


1 


√ 


5 − 1 
2 


 


 
 


x 
y 


= 
0 
0 
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es decir 
 
 
 


 
 


1 + 
√ 


5 


2 
x 
+ 
y 
= 
0 


x 
+ 


√ 


5 − 1 
2 
y 
= 
0 


Multiplicando la segunda ecuaci´on por (1 + 
√ 


5)/2, obtenemos 


1 + 
√ 


5 


2 
x + 1 + 
√ 


5 


2 


√ 


5 − 1 
2 
y = 0 =⇒ 1 + 
√ 


5 


2 
x + 5 − 1 


4 
y = 0 


=⇒ 1 + 
√ 


5 


2 
x + y = 0 


es decir, la primera ecuaci´on. Por tanto, las dos ecuaciones no son independientes. 
As´ı 


Eλ2 = 
(x, y) ∈ R2 : y = −1+ 
√ 


5 


2 
x 
= 
(x, −1+ 
√ 


5 


2 
x) : x ∈ R 


= 
x (1, −1+ 
√ 


5 


2 
) : x ∈ R 
= 
(1, −1+ 
√ 


5 


2 
) 
= 
(2, −1 − 
√ 


5) 


Entonces 


P = 
2 
2 
√ 


5 − 1 
−1 − 
√ 


5 



es la matriz de paso, y de aqu´ı 


D = P −1 · A · P 
A = P · D · P −1 


De este modo 
An−1 = P · Dn−1 · P −1 


Por tanto 
un 
un−1 


= wn−1 = An−1 · w0 = An−1 · 
1 
1 


= P · Dn−1 · P −1 · 
1 
1 



Calculamos ahora la matriz inversa de P. Omitiendo el proceso, obtenemos: 


P −1 = 


 


 
 
 


1 + 
√ 


5 


4 
√ 


5 
1 


2 
√ 


5 
√ 


5 − 1 
4 
√ 


5 
− 1 


2 
√ 


5 


 


 
 
 


Entonces 


P −1 · 
1 
1 


= 


 


 
 
 


1 + 
√ 


5 


4 
√ 


5 
1 


2 
√ 


5 
√ 


5 − 1 
4 
√ 


5 
− 1 


2 
√ 


5 


 


 
 
 


1 
1 


= 


 


 
 
 


1 + 
√ 


5 


4 
√ 


5 
+ 
1 


2 
√ 


5 
√ 


5 − 1 
4 
√ 


5 
− 
1 


2 
√ 


5 


 


 
 
 = 


 


 
 
 


3 + 
√ 


5 


4 
√ 


5 
√ 


5 − 3 
4 
√ 


5 


 


 
 
 
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Pero 


Dn−1 = 


 


 
 
 
 
 


1 + 
√ 


5 


2 


n−1 
0 


0 


1 − 
√ 


5 


2 


n−1 


 


 
 
 
 
 


y as´ı 


Dn−1 · P −1 
1 
1 


= 


 


 
 
 
 
 


1 + 
√ 


5 


2 


n−1 
0 


0 


1 − 
√ 


5 


2 


n−1 


 


 
 
 
 
 


 


 
 
 


3 + 
√ 


5 


4 
√ 


5 
√ 


5 − 3 
4 
√ 


5 


 


 
 
 


= 


 


 
 
 
 
 


3 + 
√ 


5 


4 
√ 


5 


1 + 
√ 


5 


2 


n−1 


√ 


5 − 3 
4 
√ 


5 


1 − 
√ 


5 


2 


n−1 


 


 
 
 
 
 


de donde 


P · Dn−1 · P −1 · 
1 
1 


= 
2 
2 
√ 


5 − 1 
−1 − 
√ 


5 



 


 
 
 
 
 


3 + 
√ 


5 


4 
√ 


5 


1 + 
√ 


5 


2 


n−1 


√ 


5 − 3 
4 
√ 


5 


1 − 
√ 


5 


2 


n−1 


 


 
 
 
 
 


Entonces 


un = 2 


1 + 
√ 


5 


2 


n−1 
3 + 
√ 


5 


4 
√ 


5 
+ 2 


1 − 
√ 


5 


2 


n−1 √ 


5 − 3 
4 
√ 


5 


= 


1 + 
√ 


5 


2 


n−1 
3 + 
√ 


5 


2 
√ 


5 
+ 


1 − 
√ 


5 


2 


n−1 √ 


5 − 3 
2 
√ 


5 


= 1 
√ 


5 


 


 


1 + 
√ 


5 


2 


n−1 
3 + 
√ 


5 


2 
− 


1 − 
√ 


5 


2 


n−1 
3 − 
√ 


5 


2 


 


 


= 1 
√ 


5 


 


 
1 + 
√ 


5 


2 


n−1 
1 + 
√ 


5 


2 
+ 1 



− 


1 − 
√ 


5 


2 


n−1 
1 − 
√ 


5 


2 
+ 1 


 


 


= 1 
√ 


5 


 


 


1 + 
√ 


5 


2 


n 
+ 


1 + 
√ 


5 


2 


n−1 
− 


1 − 
√ 


5 


2 


n 
− 


1 − 
√ 


5 


2 


n−1 


 
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Se puede demostrar que 
1 + 
√ 


5 


2 


n+1 
= 


1 + 
√ 


5 


2 


n 
+ 


1 + 
√ 


5 


2 


n−1 


e igualmente 
1 − 
√ 


5 


2 


n+1 
= 


1 − 
√ 


5 


2 


n 
+ 


1 − 
√ 


5 


2 


n−1 


Por tanto 


un = 
1 
√ 


5 


 


 
1 + 
√ 


5 


2 


n+1 
− 


1 − 
√ 


5 


2 


n+1 


 


As´ı, por ejemplo, transcurridos 2 a˜nos tendremos: 


u24 = 75025 pares de conejos 


y transcurridos 10 a˜nos: 


u120 = 8670007398507948658051921 pares de conejos 


es decir, del orden de 1025 conejos. 


4.5. 
Diagonalizaci´on de matrices reales sim´etricas 


Tal y como hemos visto en la secci´on anterior, no todas las matrices son dia- 
gonalizables. En cambio, cuando se trata de una matriz real y sim´etrica, caso ´este 
muy frecuente en F´ısica e Ingenier´ıa, la simetr´ıa de la matriz conduce a que ´esta sea 
diagonalizable. 
Una aplicaci´on importante de este hecho lo abordaremos en la segunda parte 
de este material, cuando consideremos problemas de extremos libres y de extre- 
mos condicionados (multiplicadores de Lagrange). En el estudio de los extremos 
(m´aximos y m´ınimos) interviene la matriz hessiana, la cual es una matriz sim´etrica. 


Teorema 4.5.1 
Toda matriz real y sim´etrica es diagonalizable. 


Ejercicio 4.5.1 
Estudia, en funci´on del par´ametro real a, si las siguientes matrices son diagonaliza- 
bles: 


(a) A = 


 


 
5 
0 
0 
0 
−1 
0 
0 
0 
a 


 


 


(b) A = 


 


 
5 
1 
0 
1 
a 
4 
0 
4 
3 


 


 
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